mathanalysis101.m

来自「数值分析的实验程序」· M 代码 · 共 36 行

format long;
e=0.5e-10;
a=1;
b=2;
%b=1.5;
while(1)
    c=(a+b)/2;
    fc=sin(c)-c.^2/2;
    fa=sin(a)-a.^2/2;
%  fc=c.^3-c-1;
%  fa=a.^3-a-1;
   if fc==0   ;
    else
        if fa.*fc>0
            a=c;
        else
            b=c;
        end
        if b-a<=e
            c=(a+b)/2;
            break;
        else
            continue;
        end
    end
end
printa=['The a is ',num2str(a,15),'.'];
printb=['The b is ',num2str(b,15),'.'];
printc=['The c is ',num2str(c,15),'.'];
printx=['So the x is ',num2str(c,15),'.'];
disp(printa);
disp(printb);
disp(printc);
disp(printx);
format short;