用MATLAB编出牛顿迭代法的程序

用MATLAB编出牛顿迭代法的程序

function newton(x0,e,N) %输入xo为估计的迭代初值,e为规定的误差,N为最大迭代次数.%输出x,y为最后迭代的两个近似根,k为迭代次数.clcformat long;disp('迭代次数 近似根')k=0;x1=0;x2=x0;while (abs(x2-x1))>ex1=x2;x2=x1-f(x1)./df(x1);k=k+1;if k>Nreturn;end%%%%%%记录并输出%%%%%%%%%o1=sprintf('%3d',k);o2=sprintf('%3.8f',x2);OL=[o1,' ' o2];disp(OL);y(k)=x2;end%%%%画图%%%%%%%i=1:k;figure(2)plot(i,y,'rD-')grid onxlabel('迭代次数')ylabel('近似根')title(['牛顿法求出的该方程的近似根 x^*=', num2str(x2,9)])function y=f(x)y=x^2/2-sin(x)-1;function y=df(x)y=x-cos(x);