fixedpoint.m

来自「国外计算机科学教材系列 numerial method using matla」· M 代码 · 共 15 行

function [k,p,err,P] = fixedpoint(g,p0,tol,max1)

P(1)= p0;

for k=2:max1
	P(k)=feval(g,P(k-1));
	err=abs(P(k)-P(k-1));
	relerr=err/(abs(P(k))+eps);
	p=P(k);
	if (err<tol) | (relerr<tol),break;end
end

P=P';