迭代求解的困惑
我遇到个问题,求x^5-3*x^3-2*x^2+2=x的解程序为
function [k,p,err,P]=fixpt(g,p0,tol,max1)
p0=1;max1=10000;tol=10^-6;
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
if k==max
disp('maximum number of iterations exceeded')
end
p=p';
可就是运行不过去 ,哪位大侠行行好,帮忙下