xjf.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function xjf(x,y,e,f)
x1=x;
y1=y;
syms t
f1=inline(f);
xx=linspace(x,y,1000);
for i=2:1000;
    yy(i)=(f1(xx(i))-f1(xx(1)));
end
yy(1)=1;
if ismember(yy,0)==zeros(1,1000);%判断是否两点的斜率为零，注意这也是局部的
    syms t2
    f3=inline(t2-f1(t2)*(t2-x1)/(f1(t2)-f1(x1)));%对迭代函数求导，判断是否满足收敛的条件
    syms t3
    y2=diff(f3(t3));
    f2=inline(y2);
    xx1=linspace(x,y,100);
    s=abs(f2(xx1));
    s1=zeros(1,100);
    s2=ones(1,100);
    s(1)=0.5;
    if s>s1 & s<s2;
        i=0;
        while abs(x-y)>e;
            w=x;
            x=y;
            y=x-f1(x)*(x-x1)/(f1(x)-f1(x1));
            i=i+1;
            fprintf('%d%13.8f\n',i,y);
        end
    else 
        disp('此多项式不能进行迭代');
    end
else
    disp('此多项式不能进行迭代(是由于迭代的函数的导数为零无意义)');
end