文档介绍：精选优质文档-----倾情为你奉上
精选优质文档-----倾情为你奉上
专心---专注---专业
专心---专注---专业
精选优质文档-----倾情为你奉上
专心---专注---专业
《数值分析》报告
运用M0))/(feval('f',x1)-feval('f',x0));
err=abs(x2-1.);
x0=x1;
x1=x2;
y=feval('f',x1);
fprintf('k=%.0f x%d=%.8f e%d=%.6e y%d=%.8e\n',k,k,x1,k,err,k,y)
if (err<e)|(y==0)|(k==max)
break;
end
end
精选优质文档-----倾情为你奉上
精选优质文档-----倾情为你奉上
专心---专注---专业
专心---专注---专业
精选优质文档-----倾情为你奉上
专心---专注---专业
function pwxmethod(f,x0,x1,x2,e,max)
%f 是要求根的方程(f(x)=0);
%x0,x1，x2是所给初值,位于x*附近；
%e是给定允许误差；
%max是迭代的最大次数；
%x3是弦截法求得的方程的近似解；
%err是误差估计；
%k是迭代次数；
%y是f(x)值
fprintf('k=%.0f x%d=%.8f y%d=%.8e\n',0,0,x0,0,feval('f',x0))
fprintf('k=%.0f x%d=%.8f y%d=%.8e\n',1,1,x1,1,feval('f',x1))
fprintf('k=%.0f x%d=%.8f y%d=%.8e\n',2,2,x2,2,feval('f',x2))
for k=3:max
f0=feval('f',x0);
f1=feval('f',x1);
f2=feval('f',x2);
a=(f0-f2)/(x0-x2);
b=(f1-f2)/(x1-x2);
c=(a-b)/(x0-x1);
w=b+c*(x2-x1);
if w<0
x3=x2-(2*f2/(w-sqrt(w^2-4*c*f2)));
elseif w>0
x3=x2-(2*f2/(w+sqrt(w^2-4*c*f2)));
end
err=abs(x3-1.);
x0=x1;
精选优质文档-----倾情为你奉上
精选优质文档-----倾情为你奉上
专心---专注---专业
专心---专注---专业
精选优质文档-----倾情为你奉上
专心---专注---专业
x1=x2;
x2=x3;
y=feval('f',x2);
fprintf('k=%.0f x%d=%.8f e%d=%.6e y%d=%.8e\n',k,k,x2,k