Advanced Search

Question: dummy index summation in maple?

Question:dummy index summation in maple?

mehdibgh 240 Maple 2016
linearalgebra indexing summation
September 12 2017
1



I did dummy indices implication using add command as below. is it OK or there are mistakes?

Is there another way to imply dummy index summation in maple instead of using add command?

any suggestion???

> 

restart
> 

II := 9:
> 

JJ := 9:
> 

with(LinearAlgebra):

 

 
> 

F := add(add(add(add(R[i, m]*R[k, m]*(U[i, j]*U[k, j]+U[i, j]*V[k, j]+V[i, j]*V[k, j]+W[i, j]*W[k, j])/((2*m+1)*(2*j+1)), i = 0 .. II), k = 0 .. II), m = 0 .. II), j = 0 .. JJ):
> 

EqU := seq(seq(diff(F, U[i, j]), j = 0 .. JJ), i = 0 .. II):

> 

EqV := seq(seq(diff(F, V[i, j]), j = 0 .. JJ), i = 0 .. II):

> 

EqW := seq(seq(diff(F, W[i, j]), j = 0 .. JJ), i = 0 .. II):

> 

var := [seq(seq(U[i, j], j = 0 .. JJ), i = 0 .. II), seq(seq(V[i, j], j = 0 .. JJ), i = 0 .. II), seq(seq(W[i, j], j = 0 .. JJ), i = 0 .. II)]:

> 

NULL

> 

sys := [EqU, EqV, EqW]:

> 

Aa, bb := GenerateMatrix(sys, var):

> 

Aa

RTABLE(18446744074191847422, anything, Matrix, rectangular, Fortran_order, [], 2, 1 .. 300, 1 .. 300)

(1.1)

``


 

Download Mesal.mw

 

 

 
Please Wait...