erfenfa.m

来自「介绍了二分法求解非线性方程的算法」· M 代码 · 共 28 行

function x3=erfenfa(x1,x2)
% function varargout=erfenfa(x1,x2)
% 用二分法求解非线性方程，
% 输入变量为x1,x2,输出变量为x3
clc
tic
nargout=3;
x1=0;x2=1;%非线性方程的区间
k=0;%迭代步数
while k<=4
% while abs(x2-x1)/2>1e-10
    x3=(x2+x1)/2;
    fx1=3*x1+sin(x1)-exp(x1);
    fx3=3*x3+sin(x3)-exp(x3);
%     fx1=x1^3+x1^2-3*x1-3;
%     fx3=x3^3+x3^2-3*x3-3;

    if fx1*fx3<0
        x2=x3;
    else
        x1=x3;
    end
    k=k+1;
end
x3
k
fx3
toc