高等数学中非线性方程求解的一种方法——牛顿法的matlab实现。

matlab中有fzero可用来求根，这里作为教学例子供学生理解牛顿法原理时参考。

 

本程序需要使用者提供至少三个输入：

fun 非线性方程

dfun 非线性方程的导函数

x0   猜想的初始根

下面是牛顿法的源代码。

function [sol,fval] = newtonalgorithm0(fun,dfun,x0,tolerance)
if  nargin < 3
    disp('Newton algorithm requires at least three input arguments.');
    return;
end
if  nargin < 4
    tolerance = 1e-6;
end
f0 = feval(fun,x0);
if abs(f0) < tolerance
    sol = x0; fval = f0; k = 0;
    return;
end
k = 0; error=1; x1=x0+1;
while error > tolerance
    f0 = feval(fun,x0);
    df0 = feval(dfun,x0);
    x1 = x0 - f0/df0;
    error = abs(x1 - x0);
    x0 = x1;
    k = k+1
end
sol = x1;
fval = feval(fun,sol);
%iters = k;