已知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编辑过]