bisection.m

来自「数值计算中的二分法和牛顿迭代法的演示」· M 代码 · 共 18 行

clear all
clc
format('long');
M = 2;
a = 1;
b = 2;
k = 0;
while b-a > eps
    x = (a + b)/2;
    if x^2 > M
        b = x;
    else
        a = x;
    end
    k = k + 1;
    disp(['a=',num2str(a),', b=',num2str(b),', x=',num2str(x),', k=',num2str(k)]);    
end
format('short');