二分法

a=1.0
b=1.5
f=sym('x^4-x-2')
e=eps(realmin).*10
while (b-a)>0.0001
	c=(a+b)/2
	if sign(subs(f,'x',c))==sign(subs(f,'x',a))
		a=c
	else
		b=c
end
end

a =

     1


b =

    1.5000

 
f =
 
x^4 - x - 2
 

e =

  4.9407e-323


c =

    1.2500


a =

    1.2500


c =

    1.3750


b =

    1.3750


c =

    1.3125


a =

    1.3125


c =

    1.3438


a =

    1.3438


c =

    1.3594


b =

    1.3594


c =

    1.3516


a =

    1.3516


c =

    1.3555


b =

    1.3555


c =

    1.3535


b =

    1.3535


c =

    1.3525


a =

    1.3525


c =

    1.3530


a =

    1.3530


c =

    1.3533


b =

    1.3533


c =

    1.3531


a =

    1.3531


c =

    1.3532


b =

    1.3532