N=70 epsilon=10**(-15) a=1.0 b=2.0 i=0 m=(a+b)/2 print i,m abserr=(b-a)/2 while i <= N and abserr > epsilon*abs(m): if m**2-2==0: a=b=m if (a**2-2)*(m**2-2)<0: b=m else: a=m i=i+1 m=(a+b)/2 print i,m abserr=(b-a)/2