Kirby, you are working too hard! > >>> simpson(g,0,3,10000) > 8.9999999999999893 Why ten thousand intervals when one interval gives the exact result? >>> simpson(g,0,3,1) 9.0 Henrik