近似值1.jpeg
近似值2.jpeg
function Newton :
function [] = Newton(x0, n)
disp(strcat('x0 ：', num2str(x0))); % 打印初值
m = 1;
while m <= n
p = [1, 0, 0, -6]; % f(x) = x^3 - 6
p1 = polyder(p); % f'(x)
x0 = x0 - polyval(p,x0)/polyval(p1,x0); % x1 = x0 - f(x0)/f'(x0)
disp(strcat('x', num2str(m), ' ：', num2str(x0, '%.9f'))); % 打印
m = m + 1;
end
end
runNewton.m
format long;
Newton(2, 5); % 初值 x0=2 ，迭代计算5次
计算结果：
>> runNewton
x0 ：2
x1 ：1.833333333
x2 ：1.817263545
x3 ：1.817120604
x4 ：1.817120593
x5 ：1.817120593
近似值3.jpeg
近似值4.jpeg
近似值5.jpeg
近似值6.jpeg
