mathanalysis102.m

数值分析的实验程序

format long
%x0=0.5;
%x0=1;
x0=0.55;
e1=0.5e-10;
e2=0.5e-10;
N=100;
count=0;
while(1)
    %fx0=abs(x0.*exp(x0)-1);
    %difffx0=exp(x0)+x0.*exp(x0);
    %%fx0=abs(x0.^3-x0-1);
    %%difffx0=3*x0.^2-1;
    fx0=abs((x0-1).^2.*(2.*x0-1));
    difffx0=2.*(x0-1).*(2.*x0-1)+(x0-1).^2.*2;
    if fx0<e1
        answer=x0;
        break;
    else 
        x1=x0-2*fx0./difffx0;
        count=count+1;
        if abs(x1-x0)<e2
            answer=x1;
        else
            if count>N
                answer=x1;
                break;
            end
        end        
        x0=x1;
    end    
end
printanswer=['The answer is ',num2str(answer,15),'.'];
disp(printanswer);