xia329.m

来自「为配套王能超版本数值分析简明教程编写,方便学生教师演示!」· M 代码 · 共 28 行

function xia329(a,b,j,f3)
a1=a;
b1=b;
f=inline(f3);
syms x
y=diff(f(x));
f2=inline(y);
xx=linspace(a,b,100);
    s=abs(f2(xx));
    s1=zeros(1,100);
    s2=ones(1,100);
    if s>s1 & s<s2;
        i=0;
        while abs(a-b)>j;
           if b<=b1 & b>=a1;
            a=b;
            b=f(a);
            i=i+1;
            fprintf('%d%13.8f\n',i,b);
        else
           disp('此多项式不能进行迭代,有能区取法有误');
           break;
       end
        end
    else 
        disp('此多项式不能进行迭代');
    end