interpolation_xiuzheng.m

给出插值算子的算法

% function x3=interpolation_xiuzheng(x1,x2)
% function varargout=interpolation(x1,x2)
clc
tic


x1=0;x2=1;
fx1=3*x1+sin(x1)-exp(x1);
fx2=3*x2+sin(x2)-exp(x2);
% fx1=x1^3+x1^2-3*x1-3;
% fx2=x2^3+x2^2-3*x2-3;

k=0;
save=fx1;F1=fx1;F2=fx2;
while k<=4 | (abs(fx3)>1e-5)
% while (abs(x2-x1)>1e-6) | (abs(fx3)>1e-10)
    x3=x2-fx2*(x2-x1)/(fx2-fx1);
    fx3=3*x3+sin(x3)-exp(x3);
    if fx1*fx3<0
        x2=x3;
        F2=fx3;
        if fx3*save>0
            F1=F1/2;
        end
    else
        x1=x3;
        F1=fx3;
        if fx3*save>0
            F2=F2/2;
        end
    end
    save=fx3;
    k=k+1;
end
x3
k
fx3
toc