newtonmweizhi.m

来自「这是非线性方程（组）的数值解法」· M 代码 · 共 18 行

function [x,n]=NewtonMweizhi(f,x0,eps)
% f 有m重根
% m未知，Newton迭代法
if nargin ==2
    eps=1e-5;
end
g=inline(diff(sym(f))); % f 的导函数
gg=inline(diff(sym(f),2)); % f 的二阶导数
x1=x0-f(x0)*g(x0)/((g(x0))^2-f(x0)*gg(x0));
n=1;
while abs(x1-x0)>eps && abs(f(x1))>eps
    x0=x1;
    x1=x0-f(x0)*g(x0)/((g(x0))^2-f(x0)*gg(x0));
    n=n+1;
end
x=x1;
end