c4_4.m

这是我从一本基于MATLAB实现的数值分析算法的书中自带的光盘程序

%C4_4.m
%用牛顿法求f(x)=0的一个根
function C4_4
ep=input('Tolerance:ep?');
x=input('Initial guess?');
a=input('多项式的系数(升幂排列):');
m=length(a);
n=0;
xb=x;
xp=0;
fprintf('It.No,    x(n-1),    y(n-1),    x(n) \n');
while abs(xb-xp)>ep
    n=n+1;
    xp=xb;
    y=0;
    for i=1:m
        y=y+a(i)*xb^(i-1);
    end
    y_d=0;
    for i=1:m-1
        y_d=y_d+i*a(i+1)*xb^(i-1);
    end
    xb=xp-y_d\y;
    fprintf('%1d,  %6.5e,  %6.5e,  %6.5e \n',n,xp,y,xb);
end
fprintf('Final solution=%2.0e \n',xb);