from math import sin,pi N=60 i=0 a=-1.0 b=20.0 eps=10**(-10) def f(x): return sin(x) m=(a+b)/2.0 abserr=(b-a)/2.0 print i,m,abserr/abs(m) while i<=N and abserr>eps*abs(m): if f(m)==0: a=b=m if f(a)*f(m)<0: b=m else: a=m i=i+1 m=(a+b)/2.0 abserr=(b-a)/2.0 print i,m,(m-4*pi)/(4*pi)