bisection1.m

来自「a file for determining a polynomial degr」· M 代码 · 共 15 行

function x3=bisection1(x1,x2)
y1=x1.^2+5*x1+6;
y2=x2.^2+5*x2+6;
while y1*y2<0&&abs(x1-x2)<.0025
    x3=(x1+x2)/2;
    y3=x3.^2+5*x3+6;
if y1*y3<0
    x2=x3;
else
    x1=x3;
end
end
disp(x3)