xia328.m

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

function xia328(a,b,j,f1)
i=0;
f=inline(f1);
while abs(a-b)>j;
    xm=a+(b-a)/2;
    if sign(f(a))==sign(f(xm));
        a=xm;
    else
        b=xm;
    end
    i=i+1;
    fprintf('%d%13.8f\n',i,a);
end
if sign(f(a))==sign(f(b))
   fprintf('error');
end