newton.m

来自「用matlab实现Romberg积分」· M 代码 · 共 8 行

function[x,k]=Newton(f,df,x0,e)
x=x0-feval('f',x0)/feval('df',x0);
k=1;
while abs(x-x0)>e & k<2000
    x0=x;
    x=x0-feval('f',x0)/feval('df',x0);
    k=k+1;
end