gexiantwo4.m

来自「数值分析方法,主要里面涵盖了计算方法方面matlab程序！」· M 代码 · 共 33 行

%GeXian
clear
clc
eps=1e-6;
delta=1e-8;
N=100;
k=0;
x1=1.0;
x0=1.9;
while(1)
    x2=x1-(func2_4_1(x1)*(x1-x0))/(func2_4_1(x1)-func2_4_1(x0));
    k=k+1;
    if(k>N|abs(x2)<eps)
        disp('GeXian method failed');
        break;
    end
    if(func2_4_1(x0)*func2_4_1(x2)<=0)
        x0=x0;
        x1=x2;
    elseif(func2_4_1(x1)*func2_4_1(x2)<0)
        x0=x1;
        x1=x2;
    end
    if abs(x1)<1
        d=x1-x0;
    else
        d=(x1-x0)/x1;
    end
    if(abs(d)<eps|abs(func2_4_1(x1))<delta)
        break;
    end
end
fprintf('%f\n',x1);