newton.m

用matlab实现Romberg积分

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