牛顿迭代法求解根号a的Matlab代码

主程序：

function [y,k]=newton(fname,dfname,x0,tol,N,m)

%输入：初值x0，最大迭代步数N，误差限tol,m=1为牛顿迭代法，m>1为修改的牛顿迭代法

%输出：近似根y，迭代次数k

y=x0;x0=y+2*tol;k=0;

while abs(x0-y)>tol&k<N

    k=k+1

    x0=y

    y=x0-m*fname(x0)/dfname(x0);

    disp(y)

end

if k==N

    disp(warning)

end

然后在命令行输入以下语句，我这里算的是根号4，你们可以自行更改

>> f=inline('x^2-4');df=inline('2*x');

>> [y,k]=newton(f,df,0.5,0.5*10^(-4),500,1)