计算平方根的另一种办法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求一个数的平方根，各种程序都会有sqrt函数，SAS当然也有。本人不才，今天在SAS文档中查其他资料时，偶然发现这样一个求解平方根的算法。
proc iml;                   /* begin IML session */

start MySqrt(x);            /* begin module */
   y = 1;                   /* initialize y */
   do until(w<1e-3);        /* begin DO loop */
      z = y;                /* set z=y */
      y = 0.5#(z+x/z);      /* estimate square root */
      w = abs(y-z);         /* compute change in estimate */
      print x y z ;
   end;                     /* end DO loop */
   return(y);               /* return approximation */
finish;

t = MySqrt({3,4,7,9});      /* call function MySqrt  */
s = sqrt({3,4,7,9});        /* compare with true values */
diff = t - s;               /* compute differences */
print t s diff;             /* print matrices */

部分输出结果：

xyz
31.73205081.7320508
422.0000001
72.64575132.645767
933.0000916

tsdiff
1.73205081.73205080
222.22E-15
2.64575132.64575134.678E-11
331.397E-9

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：平方根 Initial session module begin 平方根 start 程序 资料

欢迎大家积极共享自己的心得体会！