newtoneqs.m

来自「Gauss-Legendre 采用五点 Gauss-Legendre 求积公式计」· M 代码 · 共 11 行

function y=newtoneqs(x0)
x1=x0-fc(x0)*inv(df2(x0));
n=1;
while(norm(x1-x0)>=1.0e-6)&(n<=1000000000)
    x0=x1;
    x1=x0-fc(x0)*inv(df2(x0));
    n=n=1;
end 
n
x1