[SciPy-user] Numerical integration of a matrix

[Nils Wagner]

Hi all,

How do I integrate the following matrix elements in an efficient manner
with scipy

c_{j,k} = \int\limits_0^\pi \varepsilon a(x) \sin (k x) \sin (j x) dx   
(Latex notation)

j, k = 1, ..., 9

\varepsilon = 0.1

a(x) = x^2 (\pi-x)^2 - \delta, \delta = 201.0

Thanks in advance.


                   Nils