另一个牛顿法求解非线性方程的程序,优点是速度特别快,缺点是需要建立两个子程序,分别是待求函数子程序和待求函数的导数子程序。
主函数如下:
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