matlab利用牛顿法求解非线性方程02

 另一个牛顿法求解非线性方程的程序，优点是速度特别快，缺点是需要建立两个子程序，分别是待求函数子程序和待求函数的导数子程序。

主函数如下：

function gen=fnewton(func,dfunc,x,tol)
% func是函数表达式
% dfunc是这个函数的导数
% tol为截断误差
% x是初始值
x0=x;
d=feval(func,x0)/feval(dfunc,x0);将x0代入，利用feval函数求值
while abs(d)>tol
x1=x0-d;
x0=x1;
d=feval(func,x0)/feval(dfunc,x0);
end
gen=x0;

 待求函数子程序如下：

function fx1=func(x)
fx1=x^2-4*x-5;

 待求函数的导数子程序如下：

function fx2=dfunc(x)
fx2=2*x-4;

输入：

fnewton(@func,@dfunc,10.5,1e-6)

运行结果如下：

ans =

   5.000000000336842