sroot.m

来自「关于数值计算的matlab程序代码,包含了很多的实用的算法.程序全部可以直接用的」· M 代码 · 共 12 行

function r = sroot(a)

p0=1;

for k=1:50
   p1=(p0+a/p0)/2;
   disp(p1);
   if abs(p1-p0)/p1 < eps, break, end;
   p0=p1;
end