c4_5.m

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

%C4_5.m
%用牛顿法求f(x)=0的一个复根
function C4_5
format long e
ep=input('Tolerance:ep?');
x0=input('Initial guess x0?');
a=input('多项式的系数(升幂排列):');
m=length(a);
n=0;
xb=x0;
xp=0;
y=zeros(1,2);y_d=zeros(1,2);
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
    f=y_d\y;
    xb=xp-f;
end
x=xb