[SciPy-User] Facing a problem with integrations

[Ralf Gommers]

On Wed, Feb 29, 2012 at 3:47 AM, Marcel Caraciolo <caraciol at gmail.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 [1].
>
> 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 [1].
>
> 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.

Ralf
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20120303/c8ad11fb/attachment.html>