多元高次方程的求解 搓一把

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

已知SAS 的iml提供一元多次方程函数，多元一次方程函数

没有讲道多元高次方程函数，现在利用矩阵计算基本规律编程在SAS实现多元高次方程函数的求解

警告：由于本人金刚不坏，坏蛋拍一火车砖头我都不怕！！!

proc iml;
a={2 -3 5 4,6 2 3 -2,3 -5 -4 1,-2 5 3 -2};/*高次方程左边的系数矩阵*/
c={18881, 61261 ,27119,-17203};/*多元高次方程右边常数部分 列向量*/
x7=inv(a)*c;
x8=x7[1]**(1/4)//x7[2]**(1/3)//x7[3]**(1/2)//x7[4];
print "/*需要注意偶次方的负根*/";
print  x8 ;

/*一个假设的多元高次方程
c1=2*a**4-3*b**3+5*c**2+4*d;
c2=6*a**4+2*b**3+3*c**2-2*d;
c3=3*a**4-5*b**3-4*c**2+1*d;
c4=-2*a**4+5*b**3+3*c**2-2*d;*/

                                                             SAS 系统                   2008年03月15日 星期六 下午12时57分49秒   6

                                                     /*需要注意偶次方的负根(x8)*/

                                                 X8

                                                 10

                                                 8

                                                 9

                                                 3

[此贴子已经被作者于2008-3-15 14:29:28编辑过]

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：高次方 print 基本规律 一次方程 int 方程 求解

不错,顶一个.