请问这个程序如何理解？

shiyiming

注释里注明问号的语句，不理解是什么含义？请高手不吝赐教！
非常感谢！
proc iml;
reset deflib=resdat;
start approx(x);   /* ?  */
y=1;
do until(w<1e-3);
z=y;
y=0.5#(z+x/z);    /* ?  */
w=abs(y-z);               /* ?  */
end;
return(y);
finish approx;
t=approx({3,5,7,9});
print t;
quit;
run;

shiyiming

我想你问的并不是这个code的意思而是为什么这种算法可以找出平方根。这种算法也叫Babylonian method。是一种循环求解的方式。其初始值y可以是任意的，并不局限于某个值，比如这里的1.
京剧

shiyiming

之前没有接触过，非常感谢jingju11！