%用牛頓法求解非線性方程的根
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