from jhplot import *
import time

f1=F1D('cos(x*x)+sin(1.0/x)*x^2')
methods=['trapezium','gauss4', 'gauss8',  'richardson', 'simpson'] 
for m in methods:
   start = time.clock()
   d=f1.integral(m,5000,10,100)
   t = time.clock()-start 
   print m,' =',d,' time (ms)=',t*1000