[Edu-sig] RE: Integration correction

Kirby Urner urnerk at qwest.net
Tue Mar 29 20:49:09 CEST 2005

> Kirby got the trapezoidal integration rule wrong.

Right, I was just doing a simple average, not any trapezoid.  My rectangles
took the mean between f(x-h) and f(x+h), nothing more.  Not the best
approximation, I agree, but simple to think about, and some text books show

> This is the corrected version.
> def integrate(f,a,b,n=1000):
>     sum = 0
>     h = (b-a)/float(n)
>     for i in range(1,n):
>            sum += f(a+i*h)
>     return h*(0.5*f(a)+sum+0.5*f(b))

So now we should come up with the @integral decorator to wrap any function
in its integral.  How would that look? 


More information about the Edu-sig mailing list