|
function root= fixedpoint(func,xr,es,maxit)
if nargin<5, maxit=50;end
if nargin<4, es=0.001;end
iter=0;
while(1)
xrold=xr;
xr=func(xr)+xr; <<<
iter=iter+1;
if xr~=0, ea = abs((xr-xrold)/xr)*100; end
if ea<=es | iter>=maxit, break, end
end
disp
|
|
메뉴펼치기
