[SciPy-User] Facing a problem with integrations
Sat Mar 3 08:12:59 CST 2012
On Wed, Feb 29, 2012 at 3:47 AM, Marcel Caraciolo <email@example.com>wrote:
> Hi all,
> My name is Marcel and I am lecturing scientific computing with Python here
> at Brazil. One of my students came to me with a problem that he is
> currently solving it with matlab but he decided to change his code to
> Python (thanks to the course!)
> The problem is calculate numerically the coefficients aim that are defined
> by the following integral .
> It must be calculated using integrals. In the example showed above he
> wants to use the trapezoid rule adapted for 2-D arrays or if there is any
> another solutions easily with scipy it would be match perfectly also.
> Here is the matrix input (U) and the corresponding coefficients
> (solution). The goal is to calculate the corresponding coefficients
> by the formula (integral) shown at .
> Could anyone give some a solution using scipy.integrate ? I tried several
> proposals but it didn't worked.
It's a double integral, so your first try should be integrate.dblquad. I
suggest that you try that, then if you get stuck show us where exactly and
we'll try to help you.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-User