newton.m

来自「实现matlab的牛顿法。牛顿法是非常有高效的」· M 代码 · 共 29 行

% Created by colent 12/7
% Copyright (C) 2007-2008 XMU 数值分析 Experiment. All rights reserved.\nAuthors: 廖绮绮(ID:04120049).
% Thanks 鞠颖老师
% 牛顿迭代法

function  [] = Newton()

format long     % 显示多位有效数字
count = 0       % 记录迭代次数
x0 = 1.5        % 用户给的初值
pause;
f = x0^3-x0-1;  % 用户给的函数
ff = 3*x0^2-1;  % 用户给的函数的导数

x = x0 - f/ff   % 计算第一次迭代
xprev = x0;     % xprev表示上次的迭代值
pause;
while abs(x - xprev) > eps*abs(x)   % eps近似为2.2e-16 若两次的迭代相差很小，则说明已经解决问题
    xprev = x; 
    f = xprev^3 -xprev -1;
    ff = 3*xprev^2-1;
    x = xprev - f/ff        
    count = count +1
    pause;
end