牛顿迭代算法的matlab程序

% 用牛顿下山法求解方程function [x,k]=myfun_newton(f,x0,emg)% f表示非线形方程% x0迭代初值,此种方法是局部收敛,初值要选择恰当% emg是精度指标% k,u分别表示迭代次数和下山因子% d1表示非线形方程f在x0处的导数值[f1,d1]=feval(f,x0);k=1;x(1)=x(0);x(2)=x(1)-f1/d1;while abs(f1)>emg u=1; k=k+1; [f1,d1]=feval(f,x(k)); x(k+1)=x(k)-u*f1/d1; while abs(f2)>abs(f1) u=u/2; x(k+1)=x(k)-u*f1/d1; [f2,d2]=feval(f,x(k+1)); endend这个收敛速度快,建议给你用下