File: makeQ0-example.mac

package info (click to toggle)
maxima 5.47.0-9
  • links: PTS
  • area: main
  • in suites: forky, sid
  • size: 193,104 kB
  • sloc: lisp: 434,678; fortran: 14,665; tcl: 10,990; sh: 4,577; makefile: 2,763; ansic: 447; java: 328; python: 262; perl: 201; xml: 60; awk: 28; sed: 15; javascript: 2
file content (31 lines) | stat: -rw-r--r-- 1,532 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
 
load("makeQ0.mac")$
makeQ0(3,3,g,T)$
components(T([],[c]),g([a,b],[])*T([],[a,b,c]))$
EQ:T0([],[a,b,c])$
EQ:EQ,T$
EQ:subst(T0[a][b,c],'T0[a,b,c],ic_convert('T0([],[a,b,c])=EQ))$
lg:ident(3)$
ug:ident(3)$
dim:3$
T0:[0,0,0]$
for i thru dim do T0[i]:zeromatrix(dim,dim)$
for i thru 3 do for j thru 3 do for k thru 3 do T[i,j,k]:0$
T[3,3,3]:m*z^3$
ev(EQ)$
T0:factor(T0);
/* output shown in paper:
                                                      [      3                 ]
        [                 3 ]  [ 0    0       0    ]  [   m z                  ]
        [              m z  ]  [                   ]  [ - ----    0       0    ]
        [   0     0  - ---- ]  [                 3 ]  [    5                   ]
        [               5   ]  [              m z  ]  [                        ]
        [                   ]  [ 0    0     - ---- ]  [              3         ]
(%o15) [[   0     0    0    ], [               5   ], [           m z          ]]
        [                   ]  [                   ]  [   0     - ----    0    ]
        [      3            ]  [         3         ]  [            5           ]
        [   m z             ]  [      m z          ]  [                        ]
        [ - ----  0    0    ]  [ 0  - ----    0    ]  [                      3 ]
        [    5              ]  [       5           ]  [                 2 m z  ]
                                                      [   0       0     ------ ]
                                                      [                   5    ]
 */