newton.m

来自「自己写的matlab程序 数值分析课程里的程序」· M 代码 · 共 18 行

%用牛顿法求解非线性方程的根
function x1=newton(x0)
e=10^-7;
%求迭代函数在x0附近的最大范数L
L1=abs(fminsearch(@fung,x0))
L2=abs(fminsearch(@fungg,x0))
L=max(L1,L2)
%求函数f(x)在x0附近的根
if L<1
    %迭代函数为-fung
    x1=fung(x0);
    while norm(x1-x0)>L*e/(1-L)
        x0=x1;
        x1=fung(x0);
    end
else
    disp('diedai error')
end