[SciPy-user] Double integration

[Gary Ruben]

Hi Colin,

Try this:

In [8]: import scipy.integrate.quadpack as siq

In [9]: siq.dblquad(lambda y,x:y, 0, 1, lambda x:x, lambda x:1)
Out[9]: (0.33333333333333331, 3.7007434154171879e-015)


Gary R.

Colin Gillespie wrote:
> Dear all,
> 
> What would be the best way of solving this double integration using scipy:
> 
> int_0^1 int_x^1 y dy dx
> 
> Many thanks
> 
> Colin