본문내용
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('iter')
disp(iter)
disp('ea')
disp(ea)
disp('root')
disp(xr)
-------------------------------------------
function root= gawichi(func,xl,xu,es,maxit)
if func(xl)*func(xu)>0, error('no bracket'), return, end
if nargin<5, maxit=50;end
if nargin<4, es=0.001;end
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('iter')
disp(iter)
disp('ea')
disp(ea)
disp('root')
disp(xr)
-------------------------------------------
function root= gawichi(func,xl,xu,es,maxit)
if func(xl)*func(xu)>0, error('no bracket'), return, end
if nargin<5, maxit=50;end
if nargin<4, es=0.001;end