bisector.m

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

function d=bisector(x1,x2,n)
%y1=input('enter y1 :');
%y1=input('enter y2 :');
y1=x1.^2+5*x1+6;x=x1;
y2=x2.^2+5*x2+6;x=x2;
for i=1:n
    x3=(x1+x2)/2;y3=x3.^2+5*x3+6;x=x3;
    disp(x)
    if y1*y3<0
        x2=x3;
        if y2*y2<0
            x1=x3;
            
        end
    end
end