数值分析matlab代码(共25页).doc

1、%用牛顿法求f(x)=x-sin x 的零点，e=10(-6)disp(牛顿法);i=1;n0=180;p0=pi/3;tol=10(-6);for i=1:n0 p=p0-(p0-sin(p0)/(1-cos(p0); if abs(p-p0)=10(-6) disp(用牛顿法求得方程的根为) disp(p); disp(迭代次数为：) disp(i) break; end p0=p;endif i=n0&(abs(p-p0)=10(-6) disp(n0) disp(次牛顿迭代后无法求出方程的解)end2、disp(Steffensen加速);p0=pi/3;for i=1:n0 p1=0.5*p0+0.5*cos(p0); p2=0.5*p1+0.5*cos(p1); p=p0-(p1-p0).2)./(