Matlab 牛顿方法求解方程

function [jie ,n]= newton(f,g,x0,tol)
%输出：jie为牛顿法函数解；n为迭代次数；
%输入：f为求解函数，g为其导数，x0为初始估计，tol为控制精度。
if g(x0)==0
    error('初始值不合适');
end
%避免出现导致牛顿法导数分母为零情况。
i=1;x(1)=x0;
x(i+1)=x(i)-f(x(i))/g(x(i));
while (abs(x(i+1)-x(i))>tol)
    i=i+1;
    x(i+1)=x(i)-f(x(i))/g(x(i));
end
jie=x(i+1);n=i;
end

例：求解函数在区间【0，3】上的根
>> f=@(x)14*x.*exp(x-2)-12*exp(x-2)-7*x.^3+20*x.^2-26*x+12;
>> g=@(x)40*x + 2*exp(x - 2) + 14*x*exp(x - 2) - 21*x^2 - 26;
>> [jie,n]=newton(f,g,0.9,0.5*10^-6)

jie =

   0.857142857142857

n =

     4