用SAS直接解方程

求解方程,在matlab中简直是小菜一碟,并且还提供数值解和解析解,功能很强大。当然,R也不错的功能强大的求解备选软件之一。

以前幻想着SAS里面也可以做得到,曾经在SASor发帖咨询,未果!【当然IML和用类似于C语言一样的算法求解除外,但不具备第四代语言的风格。】

作为号称第四代语言的SAS,现在9.2版本可以做到了:

代码示例:解方程:y=x+x**2,y=25,求x值。

proc fcmp;
      /* define the function */
   function inversesqrt(x);
      return(x+x**2);
   endsub;

   y=25;
   x=solve("inversesqrt", {.}, y, .);
   put x;
run;

最后得到输出值:

4.5249378106

注:1.不可解方程组;2.不可解虚根,一律missing。


相关博文

《用SAS直接解方程》上的一个想法

发表评论

电子邮件地址不会被公开。 必填项已用*标注