Question: Any comment, idea or innovation to calculate this parametric integral?

Any comment, idea or innovation to calculate this parametric integral?

Note M, II, JJ are arbitrary positive integers (0<M, II, JJ<11).

F must be a function of Pm at the final!

``

restart

M := 3:

JJ := 5

II := 5

with(ArrayTools):

W := RandomArray(II, JJ, M):

V := ArrayTools:-RandomArray(II, JJ, M):

w := add(add(add(W[i, j, m]*LegendreP(i-1, x)*LegendreP(j-1, y)*p[m], i = 1 .. II), j = 1 .. JJ), m = 1 .. M):

v := add(add(add(V[i, j, m]*LegendreP(i-1, x)*LegendreP(j-1, y)*p[m], i = 1 .. II), j = 1 .. JJ), m = 1 .. M):

L := add(add((LegendreP(i-1, x)*LegendreP(j-1, y))^2, i = 1 .. II), j = 1 .. JJ):

H := 1+tanh(w-v)

F := int(H*L, [y = -1 .. 1, x = -1 .. 1])

Warning,  computation interrupted

 

``

``

Download integralproblem.mw

Please Wait...