Solution 4 to problem over
Expressions |
Parameters |
Inequalities |
Relevance |
Back to problem over
Expressions
The solution is given through the following expressions:
r10=0
r11=0
r12=0
r13=0
r14=0
r15=0
2
b11 *r29
r20=----------
2
4*a33
r21=0
2 2
a11*b11 *r29 - a33*b11 *r29
r22=-----------------------------
2 2
4*a11*a33 - 4*a22*a33
r23=0
r24=0
r26=0
r27=0
r28=0
r210=0
r212=0
r213=0
a11*r29 - a22*r29
r214=-------------------
a11 - a33
r215=0
r216=0
r218=0
r219=0
r30=0
r31=0
r32=0
r33=0
r34=0
r35=0
r36=0
r37=0
r38=0
r39=0
r310=0
r311=0
r312=0
r313=0
r314=0
r315=0
r316=0
r317=0
r318=0
r319=0
r320=0
r321=0
r322=0
r323=0
r324=0
r325=0
r326=0
r327=0
r328=0
r329=0
r330=0
r331=0
r332=0
r333=0
r334=0
r335=0
r336=0
r337=0
r338=0
r339=0
r340=0
r341=0
r342=0
r343=0
r344=0
r345=0
r346=0
r347=0
r348=0
r349=0
r350=0
r351=0
r352=0
r353=0
r354=0
r355=0
5 2 2 5 4 4 2 2
r40=( - 4*a11 *a33 *b11 *r446 - a11 *b11 *r464 + 20*a11 *a22*a33 *b11 *r446
4 3 4 4 4 4
+ 8*a11 *a22*a33*b11 *r4102 - 3*a11 *a22*b11 *r4119 - 2*a11 *a22*b11 *r464
4 2 3 4 4 4 4
- 8*a11 *a33 *b11 *r4102 + 3*a11 *a33*b11 *r4119 + 7*a11 *a33*b11 *r464
3 2 2 2 3 2 3
- 36*a11 *a22 *a33 *b11 *r446 - 16*a11 *a22 *a33*b11 *r4102
3 2 4 3 2 4
+ 3*a11 *a22 *b11 *r4119 + 8*a11 *a22 *b11 *r464
3 3 2 3 4
- 8*a11 *a22*a33 *b11 *r446 + 6*a11 *a22*a33*b11 *r4119
3 4 3 4 2
- 8*a11 *a22*a33*b11 *r464 + 4*a11 *a33 *b11 *r446
3 3 3 3 2 4
+ 16*a11 *a33 *b11 *r4102 - 9*a11 *a33 *b11 *r4119
3 2 4 2 3 2 2
- 10*a11 *a33 *b11 *r464 + 28*a11 *a22 *a33 *b11 *r446
2 3 3 2 3 4
+ 8*a11 *a22 *a33*b11 *r4102 - a11 *a22 *b11 *r4119
2 3 4 2 2 3 2
- 6*a11 *a22 *b11 *r464 + 24*a11 *a22 *a33 *b11 *r446
2 2 2 3 2 2 4
+ 24*a11 *a22 *a33 *b11 *r4102 - 6*a11 *a22 *a33*b11 *r4119
2 2 4 2 4 2
- 6*a11 *a22 *a33*b11 *r464 - 12*a11 *a22*a33 *b11 *r446
2 3 3 2 2 4
- 24*a11 *a22*a33 *b11 *r4102 - 3*a11 *a22*a33 *b11 *r4119
2 2 4 2 4 3
+ 18*a11 *a22*a33 *b11 *r464 - 8*a11 *a33 *b11 *r4102
2 3 4 2 3 4
+ 10*a11 *a33 *b11 *r4119 + 4*a11 *a33 *b11 *r464
4 2 2 4 4 4 4
- 8*a11*a22 *a33 *b11 *r446 + a11*a22 *b11 *r4119 + a11*a22 *b11 *r464
3 3 2 3 2 3
- 24*a11*a22 *a33 *b11 *r446 - 16*a11*a22 *a33 *b11 *r4102
3 4 3 4
- 2*a11*a22 *a33*b11 *r4119 + 8*a11*a22 *a33*b11 *r464
2 4 2 2 2 4
+ 12*a11*a22 *a33 *b11 *r446 + 9*a11*a22 *a33 *b11 *r4119
2 2 4 4 3
- 6*a11*a22 *a33 *b11 *r464 + 16*a11*a22*a33 *b11 *r4102
3 4 3 4
- 4*a11*a22*a33 *b11 *r4119 - 8*a11*a22*a33 *b11 *r464
4 4 4 3 2 4 4
- 4*a11*a33 *b11 *r4119 + 8*a22 *a33 *b11 *r446 - a22 *a33*b11 *r4119
4 4 3 4 2 3 3 3
- a22 *a33*b11 *r464 - 4*a22 *a33 *b11 *r446 + 8*a22 *a33 *b11 *r4102
3 2 4 3 2 4 2 4 3
+ 3*a22 *a33 *b11 *r4119 - 2*a22 *a33 *b11 *r464 - 8*a22 *a33 *b11 *r4102
2 3 4 2 3 4 4 4
- 6*a22 *a33 *b11 *r4119 + 4*a22 *a33 *b11 *r464 + 4*a22*a33 *b11 *r4119)/
5 4 4 4 3 2 4 2 3 4
(32*a11 *a33 - 160*a11 *a22*a33 + 320*a11 *a22 *a33 - 320*a11 *a22 *a33
4 4 5 4
+ 160*a11*a22 *a33 - 32*a22 *a33 )
r41=0
5 2 2 5 4 4 2 2
r42=( - 4*a11 *a33 *b11 *r446 - a11 *b11 *r464 + 16*a11 *a22*a33 *b11 *r446
4 3 4 4
+ 4*a11 *a22*a33*b11 *r4102 - 2*a11 *a22*b11 *r4119
4 3 2 4 2 3 4 4
+ 4*a11 *a33 *b11 *r446 - 4*a11 *a33 *b11 *r4102 + 2*a11 *a33*b11 *r4119
4 4 3 2 2 2
+ 5*a11 *a33*b11 *r464 - 22*a11 *a22 *a33 *b11 *r446
3 2 3 3 2 4 3 2 4
- 4*a11 *a22 *a33*b11 *r4102 + a11 *a22 *b11 *r4119 + a11 *a22 *b11 *r464
3 3 2 3 2 3
- 20*a11 *a22*a33 *b11 *r446 - 8*a11 *a22*a33 *b11 *r4102
3 4 3 4
+ 6*a11 *a22*a33*b11 *r4119 - 2*a11 *a22*a33*b11 *r464
3 4 2 3 3 3 3 2 4
+ 2*a11 *a33 *b11 *r446 + 12*a11 *a33 *b11 *r4102 - 7*a11 *a33 *b11 *r4119
3 2 4 2 3 2 2
- 9*a11 *a33 *b11 *r464 + 10*a11 *a22 *a33 *b11 *r446
2 3 3 2 3 4
- 4*a11 *a22 *a33*b11 *r4102 - a11 *a22 *b11 *r4119
2 3 4 2 2 3 2
+ 2*a11 *a22 *b11 *r464 + 36*a11 *a22 *a33 *b11 *r446
2 2 2 3 2 2 4
+ 24*a11 *a22 *a33 *b11 *r4102 - 9*a11 *a22 *a33*b11 *r464
2 4 2 2 3 3
- 6*a11 *a22*a33 *b11 *r446 - 12*a11 *a22*a33 *b11 *r4102
2 2 4 2 2 4
- 9*a11 *a22*a33 *b11 *r4119 + 12*a11 *a22*a33 *b11 *r464
2 4 3 2 3 4 2 3 4
- 8*a11 *a33 *b11 *r4102 + 10*a11 *a33 *b11 *r4119 + 5*a11 *a33 *b11 *r464
4 2 2 4 3
+ 2*a11*a22 *a33 *b11 *r446 + 4*a11*a22 *a33*b11 *r4102
4 4 4 4
+ 2*a11*a22 *b11 *r4119 - 2*a11*a22 *b11 *r464
3 3 2 3 2 3
- 28*a11*a22 *a33 *b11 *r446 - 8*a11*a22 *a33 *b11 *r4102
3 4 3 4
- 6*a11*a22 *a33*b11 *r4119 + 4*a11*a22 *a33*b11 *r464
2 4 2 2 3 3
+ 6*a11*a22 *a33 *b11 *r446 - 12*a11*a22 *a33 *b11 *r4102
2 2 4 2 2 4
+ 9*a11*a22 *a33 *b11 *r4119 + 3*a11*a22 *a33 *b11 *r464
4 3 3 4
+ 16*a11*a22*a33 *b11 *r4102 - 10*a11*a22*a33 *b11 *r464
4 4 5 2 2 4 3 2
- 5*a11*a33 *b11 *r4119 - 2*a22 *a33 *b11 *r446 + 8*a22 *a33 *b11 *r446
4 2 3 4 4 4 4
- 4*a22 *a33 *b11 *r4102 - 2*a22 *a33*b11 *r4119 + 2*a22 *a33*b11 *r464
3 4 2 3 3 3 3 2 4
- 2*a22 *a33 *b11 *r446 + 12*a22 *a33 *b11 *r4102 + 7*a22 *a33 *b11 *r4119
3 2 4 2 4 3 2 3 4
- 6*a22 *a33 *b11 *r464 - 8*a22 *a33 *b11 *r4102 - 10*a22 *a33 *b11 *r4119
2 3 4 4 4 5 4
+ 5*a22 *a33 *b11 *r464 + 5*a22*a33 *b11 *r4119)/(16*a11 *a33
4 4 3 2 4 2 3 4
- 80*a11 *a22*a33 + 160*a11 *a22 *a33 - 160*a11 *a22 *a33
4 4 5 4
+ 80*a11*a22 *a33 - 16*a22 *a33 )
r43=0
5 2 2 5 4 4 2 2
r44=( - 4*a11 *a33 *b11 *r446 - a11 *b11 *r464 + 12*a11 *a22*a33 *b11 *r446
4 4 4 4 4 3 2
- a11 *a22*b11 *r4119 + 2*a11 *a22*b11 *r464 + 8*a11 *a33 *b11 *r446
4 4 4 4 3 2 2 2
+ a11 *a33*b11 *r4119 + 3*a11 *a33*b11 *r464 - 8*a11 *a22 *a33 *b11 *r446
3 2 3 3 2 4
+ 8*a11 *a22 *a33*b11 *r4102 + a11 *a22 *b11 *r4119
3 2 4 3 3 2
- 5*a11 *a22 *b11 *r464 - 32*a11 *a22*a33 *b11 *r446
3 2 3 3 4
- 16*a11 *a22*a33 *b11 *r4102 + 2*a11 *a22*a33*b11 *r4119
3 4 3 3 3
+ 2*a11 *a22*a33*b11 *r464 + 8*a11 *a33 *b11 *r4102
3 2 4 3 2 4
- 3*a11 *a33 *b11 *r4119 - 7*a11 *a33 *b11 *r464
2 3 2 2 2 3 3
- 8*a11 *a22 *a33 *b11 *r446 - 16*a11 *a22 *a33*b11 *r4102
2 3 4 2 3 4
- 4*a11 *a22 *b11 *r4119 + 8*a11 *a22 *b11 *r464
2 2 3 2 2 2 2 3
+ 48*a11 *a22 *a33 *b11 *r446 + 24*a11 *a22 *a33 *b11 *r4102
2 2 4 2 2 4
+ 9*a11 *a22 *a33*b11 *r4119 - 9*a11 *a22 *a33*b11 *r464
2 2 4 2 2 4
- 12*a11 *a22*a33 *b11 *r4119 + 6*a11 *a22*a33 *b11 *r464
2 4 3 2 3 4 2 3 4
- 8*a11 *a33 *b11 *r4102 + 7*a11 *a33 *b11 *r4119 + 5*a11 *a33 *b11 *r464
4 2 2 4 3
+ 12*a11*a22 *a33 *b11 *r446 + 8*a11*a22 *a33*b11 *r4102
4 4 4 4
+ 4*a11*a22 *b11 *r4119 - 4*a11*a22 *b11 *r464
3 3 2 3 4
- 32*a11*a22 *a33 *b11 *r446 - 8*a11*a22 *a33*b11 *r4119
2 3 3 2 2 4
- 24*a11*a22 *a33 *b11 *r4102 + 3*a11*a22 *a33 *b11 *r4119
2 2 4 4 3
+ 9*a11*a22 *a33 *b11 *r464 + 16*a11*a22*a33 *b11 *r4102
3 4 3 4
+ 6*a11*a22*a33 *b11 *r4119 - 10*a11*a22*a33 *b11 *r464
4 4 5 2 2 4 3 2
- 5*a11*a33 *b11 *r4119 - 4*a22 *a33 *b11 *r446 + 8*a22 *a33 *b11 *r446
4 2 3 4 4 4 4
- 8*a22 *a33 *b11 *r4102 - 4*a22 *a33*b11 *r4119 + 4*a22 *a33*b11 *r464
3 3 3 3 2 4
+ 16*a22 *a33 *b11 *r4102 + 12*a22 *a33 *b11 *r4119
3 2 4 2 4 3 2 3 4
- 8*a22 *a33 *b11 *r464 - 8*a22 *a33 *b11 *r4102 - 13*a22 *a33 *b11 *r4119
2 3 4 4 4 5 4
+ 5*a22 *a33 *b11 *r464 + 5*a22*a33 *b11 *r4119)/(32*a11 *a33
4 4 3 2 4 2 3 4
- 160*a11 *a22*a33 + 320*a11 *a22 *a33 - 320*a11 *a22 *a33
4 4 5 4
+ 160*a11*a22 *a33 - 32*a22 *a33 )
r45=0
r46=0
r47=0
r48=0
5 2 2 5 4 4 2 2
r49=( - 2*a11 *a33 *b11 *r446 - a11 *b11 *r464 + 10*a11 *a22*a33 *b11 *r446
4 3 4 4 4 4
+ 4*a11 *a22*a33*b11 *r4102 - a11 *a22*b11 *r4119 + a11 *a22*b11 *r464
4 2 3 4 4 4 4
- 4*a11 *a33 *b11 *r4102 + a11 *a33*b11 *r4119 + 4*a11 *a33*b11 *r464
3 2 2 2 3 2 3
- 18*a11 *a22 *a33 *b11 *r446 - 8*a11 *a22 *a33*b11 *r4102
3 2 4 3 2 4 3 3 2
+ a11 *a22 *b11 *r4119 + a11 *a22 *b11 *r464 - 4*a11 *a22*a33 *b11 *r446
3 4 3 4
+ 2*a11 *a22*a33*b11 *r4119 - 6*a11 *a22*a33*b11 *r464
3 4 2 3 3 3 3 2 4
+ 2*a11 *a33 *b11 *r446 + 8*a11 *a33 *b11 *r4102 - 3*a11 *a33 *b11 *r4119
3 2 4 2 3 2 2
- 5*a11 *a33 *b11 *r464 + 14*a11 *a22 *a33 *b11 *r446
2 3 3 2 3 4 2 3 4
+ 4*a11 *a22 *a33*b11 *r4102 - a11 *a22 *b11 *r4119 - a11 *a22 *b11 *r464
2 2 3 2 2 2 2 3
+ 12*a11 *a22 *a33 *b11 *r446 + 12*a11 *a22 *a33 *b11 *r4102
2 4 2 2 3 3
- 6*a11 *a22*a33 *b11 *r446 - 12*a11 *a22*a33 *b11 *r4102
2 2 4 2 2 4
- 3*a11 *a22*a33 *b11 *r4119 + 9*a11 *a22*a33 *b11 *r464
2 4 3 2 3 4 2 3 4
- 4*a11 *a33 *b11 *r4102 + 4*a11 *a33 *b11 *r4119 + 2*a11 *a33 *b11 *r464
4 2 2 4 4
- 4*a11*a22 *a33 *b11 *r446 + a11*a22 *b11 *r4119
3 3 2 3 2 3
- 12*a11*a22 *a33 *b11 *r446 - 8*a11*a22 *a33 *b11 *r4102
3 4 3 4
- 2*a11*a22 *a33*b11 *r4119 + 2*a11*a22 *a33*b11 *r464
2 4 2 2 2 4
+ 6*a11*a22 *a33 *b11 *r446 + 3*a11*a22 *a33 *b11 *r4119
2 2 4 4 3
- 3*a11*a22 *a33 *b11 *r464 + 8*a11*a22*a33 *b11 *r4102
3 4 4 4 4 3 2
- 4*a11*a22*a33 *b11 *r464 - 2*a11*a33 *b11 *r4119 + 4*a22 *a33 *b11 *r446
4 4 3 4 2 3 3 3
- a22 *a33*b11 *r4119 - 2*a22 *a33 *b11 *r446 + 4*a22 *a33 *b11 *r4102
3 2 4 3 2 4 2 4 3
+ 3*a22 *a33 *b11 *r4119 - a22 *a33 *b11 *r464 - 4*a22 *a33 *b11 *r4102
2 3 4 2 3 4 4 4
- 4*a22 *a33 *b11 *r4119 + 2*a22 *a33 *b11 *r464 + 2*a22*a33 *b11 *r4119)/
5 4 4 4 3 2 4 2 3 4
(16*a11 *a33 - 80*a11 *a22*a33 + 160*a11 *a22 *a33 - 160*a11 *a22 *a33
4 4 5 4
+ 80*a11*a22 *a33 - 16*a22 *a33 )
r410=0
5 2 2 5 4 4 2 2
r411=( - 2*a11 *a33 *b11 *r446 - a11 *b11 *r464 + 6*a11 *a22*a33 *b11 *r446
4 4 4 4 4 3 2
- a11 *a22*b11 *r4119 + 2*a11 *a22*b11 *r464 + 4*a11 *a33 *b11 *r446
4 4 4 4 3 2 2 2
+ a11 *a33*b11 *r4119 + 3*a11 *a33*b11 *r464 - 4*a11 *a22 *a33 *b11 *r446
3 2 3 3 2 4
+ 4*a11 *a22 *a33*b11 *r4102 + 2*a11 *a22 *b11 *r4119
3 2 4 3 3 2
- 3*a11 *a22 *b11 *r464 - 16*a11 *a22*a33 *b11 *r446
3 2 3 3 4
- 8*a11 *a22*a33 *b11 *r4102 - 2*a11 *a22*a33*b11 *r464
3 3 3 3 2 4 3 2 4
+ 4*a11 *a33 *b11 *r4102 - 2*a11 *a33 *b11 *r4119 - 5*a11 *a33 *b11 *r464
2 3 2 2 2 3 3
- 4*a11 *a22 *a33 *b11 *r446 - 8*a11 *a22 *a33*b11 *r4102
2 3 4 2 3 4
- 3*a11 *a22 *b11 *r4119 + 4*a11 *a22 *b11 *r464
2 2 3 2 2 2 2 3
+ 24*a11 *a22 *a33 *b11 *r446 + 12*a11 *a22 *a33 *b11 *r4102
2 2 4 2 2 4
+ 3*a11 *a22 *a33*b11 *r4119 - 3*a11 *a22 *a33*b11 *r464
2 2 4 2 2 4
- 3*a11 *a22*a33 *b11 *r4119 + 6*a11 *a22*a33 *b11 *r464
2 4 3 2 3 4 2 3 4
- 4*a11 *a33 *b11 *r4102 + 3*a11 *a33 *b11 *r4119 + 3*a11 *a33 *b11 *r464
4 2 2 4 3
+ 6*a11*a22 *a33 *b11 *r446 + 4*a11*a22 *a33*b11 *r4102
4 4 4 4
+ 2*a11*a22 *b11 *r4119 - 2*a11*a22 *b11 *r464
3 3 2 3 4
- 16*a11*a22 *a33 *b11 *r446 - 2*a11*a22 *a33*b11 *r4119
2 3 3 2 2 4
- 12*a11*a22 *a33 *b11 *r4102 + 3*a11*a22 *a33 *b11 *r464
4 3 3 4
+ 8*a11*a22*a33 *b11 *r4102 + 2*a11*a22*a33 *b11 *r4119
3 4 4 4
- 6*a11*a22*a33 *b11 *r464 - 2*a11*a33 *b11 *r4119
5 2 2 4 3 2 4 2 3
- 2*a22 *a33 *b11 *r446 + 4*a22 *a33 *b11 *r446 - 4*a22 *a33 *b11 *r4102
4 4 4 4 3 3 3
- 2*a22 *a33*b11 *r4119 + 2*a22 *a33*b11 *r464 + 8*a22 *a33 *b11 *r4102
3 2 4 3 2 4 2 4 3
+ 5*a22 *a33 *b11 *r4119 - 4*a22 *a33 *b11 *r464 - 4*a22 *a33 *b11 *r4102
2 3 4 2 3 4 4 4
- 5*a22 *a33 *b11 *r4119 + 3*a22 *a33 *b11 *r464 + 2*a22*a33 *b11 *r4119)
5 4 4 4 3 2 4 2 3 4
/(16*a11 *a33 - 80*a11 *a22*a33 + 160*a11 *a22 *a33 - 160*a11 *a22 *a33
4 4 5 4
+ 80*a11*a22 *a33 - 16*a22 *a33 )
r412=0
r413=0
4 2 4 3 3 2
r415=(a11 *a33*b11 *r4102 - a11 *b11 *r464 + 2*a11 *a22*a33 *b11*r446
3 2 3 3 3 3
+ a11 *a22*a33*b11 *r4102 - a11 *a22*b11 *r4119 + a11 *a22*b11 *r464
3 3 3 2 2 3 3
- 2*a11 *a33 *b11*r446 - 5*a11 *a33 *b11 *r4102 + a11 *a33*b11 *r4119
3 3 2 2 2
+ 3*a11 *a33*b11 *r464 - 6*a11 *a22 *a33 *b11*r446
2 2 2 2 2 3
- 5*a11 *a22 *a33*b11 *r4102 + a11 *a22 *b11 *r464
2 3 2 2 2
+ 6*a11 *a22*a33 *b11*r446 + 7*a11 *a22*a33 *b11 *r4102
2 3 2 3
+ 3*a11 *a22*a33*b11 *r4119 - 5*a11 *a22*a33*b11 *r464
2 3 2 2 2 3 2 2 3
+ 4*a11 *a33 *b11 *r4102 - 3*a11 *a33 *b11 *r4119 - 2*a11 *a33 *b11 *r464
3 2 3 2
+ 6*a11*a22 *a33 *b11*r446 + 3*a11*a22 *a33*b11 *r4102
3 3 3 3 2 3
+ a11*a22 *b11 *r4119 - a11*a22 *b11 *r464 - 6*a11*a22 *a33 *b11*r446
2 2 2 2 3
+ a11*a22 *a33 *b11 *r4102 - 3*a11*a22 *a33*b11 *r4119
2 3 3 2
+ a11*a22 *a33*b11 *r464 - 8*a11*a22*a33 *b11 *r4102
2 3 3 3 4 2
+ 4*a11*a22*a33 *b11 *r464 + 2*a11*a33 *b11 *r4119 - 2*a22 *a33 *b11*r446
3 3 3 2 2 3 3
+ 2*a22 *a33 *b11*r446 - 3*a22 *a33 *b11 *r4102 - a22 *a33*b11 *r4119
3 3 2 3 2 2 2 3
+ a22 *a33*b11 *r464 + 4*a22 *a33 *b11 *r4102 + 3*a22 *a33 *b11 *r4119
2 2 3 3 3 4 3
- 2*a22 *a33 *b11 *r464 - 2*a22*a33 *b11 *r4119)/(4*a11 *a33
3 3 2 2 3 3 3 4 3
- 16*a11 *a22*a33 + 24*a11 *a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 )
r416=0
4 2 4 3 3 2
r417=(a11 *a33*b11 *r4102 - a11 *b11 *r464 + 2*a11 *a22*a33 *b11*r446
3 2 3 3 3 3
+ 2*a11 *a22*a33*b11 *r4102 - a11 *a22*b11 *r4119 - 2*a11 *a33 *b11*r446
3 2 2 3 3 3 3
- 6*a11 *a33 *b11 *r4102 + a11 *a33*b11 *r4119 + 4*a11 *a33*b11 *r464
2 2 2 2 2 2
- 6*a11 *a22 *a33 *b11*r446 - 7*a11 *a22 *a33*b11 *r4102
2 2 3 2 2 3 2 3
- a11 *a22 *b11 *r4119 + 3*a11 *a22 *b11 *r464 + 6*a11 *a22*a33 *b11*r446
2 2 2 2 3
+ 8*a11 *a22*a33 *b11 *r4102 + 5*a11 *a22*a33*b11 *r4119
2 3 2 3 2
- 6*a11 *a22*a33*b11 *r464 + 5*a11 *a33 *b11 *r4102
2 2 3 2 2 3
- 4*a11 *a33 *b11 *r4119 - 3*a11 *a33 *b11 *r464
3 2 3 2
+ 6*a11*a22 *a33 *b11*r446 + 4*a11*a22 *a33*b11 *r4102
3 3 3 3 2 3
+ 2*a11*a22 *b11 *r4119 - 2*a11*a22 *b11 *r464 - 6*a11*a22 *a33 *b11*r446
2 2 2 2 3
+ 2*a11*a22 *a33 *b11 *r4102 - 4*a11*a22 *a33*b11 *r4119
3 2 2 3
- 10*a11*a22*a33 *b11 *r4102 - a11*a22*a33 *b11 *r4119
2 3 3 3 4 2
+ 6*a11*a22*a33 *b11 *r464 + 3*a11*a33 *b11 *r4119 - 2*a22 *a33 *b11*r446
3 3 3 2 2 3 3
+ 2*a22 *a33 *b11*r446 - 4*a22 *a33 *b11 *r4102 - 2*a22 *a33*b11 *r4119
3 3 2 3 2 2 2 3
+ 2*a22 *a33*b11 *r464 + 5*a22 *a33 *b11 *r4102 + 5*a22 *a33 *b11 *r4119
2 2 3 3 3 4 3
- 3*a22 *a33 *b11 *r464 - 3*a22*a33 *b11 *r4119)/(4*a11 *a33
3 3 2 2 3 3 3 4 3
- 16*a11 *a22*a33 + 24*a11 *a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 )
r418=0
r419=0
r420=0
r421=0
2 2 2 2 2 3
r422=(a11 *a22*a33 *b11*r446 + 2*a11 *a22*a33*b11 *r4102 - a11 *a22*b11 *r464
2 3 2 2 2 2 3
- a11 *a33 *b11*r446 - 2*a11 *a33 *b11 *r4102 + a11 *a33*b11 *r464
2 2 2 2
- 2*a11*a22 *a33 *b11*r446 - 2*a11*a22 *a33*b11 *r4102
2 3 2 3 3
- a11*a22 *b11 *r4119 + a11*a22 *b11 *r464 + 2*a11*a22*a33 *b11*r446
3 3 2 2 3
+ 2*a11*a22*a33*b11 *r4119 + 2*a11*a33 *b11 *r4102 - a11*a33 *b11 *r4119
2 3 3 2 2 3
- a11*a33 *b11 *r464 + a22 *a33 *b11*r446 - a22 *a33 *b11*r446
2 2 2 2 3 2 3
+ 2*a22 *a33 *b11 *r4102 + a22 *a33*b11 *r4119 - a22 *a33*b11 *r464
3 2 2 3 2 3
- 2*a22*a33 *b11 *r4102 - 2*a22*a33 *b11 *r4119 + a22*a33 *b11 *r464
3 3 3 3 2 3 2 3
+ a33 *b11 *r4119)/(2*a11 *a33 - 6*a11 *a22*a33 + 6*a11*a22 *a33
3 3
- 2*a22 *a33 )
r423=0
r424=0
3 3 2 3 2
r425=(4*a11 *a22*a33*b11*r4102 - 2*a11 *a22*b11 *r4119 - 3*a11 *a22*b11 *r464
3 2 3 2 3 2
- 4*a11 *a33 *b11*r4102 + 2*a11 *a33*b11 *r4119 + 3*a11 *a33*b11 *r464
2 2 2 2 2
- 8*a11 *a22 *a33*b11*r4102 + a11 *a22 *b11 *r4119
2 2 2 2 2
+ 6*a11 *a22 *b11 *r464 + 4*a11 *a22*a33 *b11*r4102
2 2 2 2
+ 4*a11 *a22*a33*b11 *r4119 - 3*a11 *a22*a33*b11 *r464
2 3 2 2 2 2 2 2
+ 4*a11 *a33 *b11*r4102 - 5*a11 *a33 *b11 *r4119 - 3*a11 *a33 *b11 *r464
3 3 2 3 2
+ 4*a11*a22 *a33*b11*r4102 + a11*a22 *b11 *r4119 - 3*a11*a22 *b11 *r464
2 2 2 2
+ 4*a11*a22 *a33 *b11*r4102 - 5*a11*a22 *a33*b11 *r4119
2 2 3
- 3*a11*a22 *a33*b11 *r464 - 8*a11*a22*a33 *b11*r4102
2 2 2 2
+ a11*a22*a33 *b11 *r4119 + 6*a11*a22*a33 *b11 *r464
3 2 3 2 3 2
+ 3*a11*a33 *b11 *r4119 - 4*a22 *a33 *b11*r4102 - a22 *a33*b11 *r4119
3 2 2 3 2 2 2
+ 3*a22 *a33*b11 *r464 + 4*a22 *a33 *b11*r4102 + 4*a22 *a33 *b11 *r4119
2 2 2 3 2 4 2
- 3*a22 *a33 *b11 *r464 - 3*a22*a33 *b11 *r4119)/(4*a11 *a33
3 2 2 2 2 3 2 4 2
- 16*a11 *a22*a33 + 24*a11 *a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 )
r426=0
4 2 3 2 3
r427=( - 2*a11 *a33 *r446 + 8*a11 *a22*a33 *r446 + 4*a11 *a22*a33*b11*r4102
3 2 3 2 3 2
- a11 *a22*b11 *r4119 - 2*a11 *a22*b11 *r464 - 4*a11 *a33 *b11*r4102
3 2 3 2 2 2 2
+ a11 *a33*b11 *r4119 + 2*a11 *a33*b11 *r464 - 12*a11 *a22 *a33 *r446
2 2 2 2 2
- 8*a11 *a22 *a33*b11*r4102 - a11 *a22 *b11 *r4119
2 2 2 2 2
+ 4*a11 *a22 *b11 *r464 + 4*a11 *a22*a33 *b11*r4102
2 2 2 2
+ 5*a11 *a22*a33*b11 *r4119 - 2*a11 *a22*a33*b11 *r464
2 3 2 2 2 2 2 2
+ 4*a11 *a33 *b11*r4102 - 4*a11 *a33 *b11 *r4119 - 2*a11 *a33 *b11 *r464
3 2 3 3 2
+ 8*a11*a22 *a33 *r446 + 4*a11*a22 *a33*b11*r4102 + 2*a11*a22 *b11 *r4119
3 2 2 2
- 2*a11*a22 *b11 *r464 + 4*a11*a22 *a33 *b11*r4102
2 2 2 2
- 4*a11*a22 *a33*b11 *r4119 - 2*a11*a22 *a33*b11 *r464
3 2 2
- 8*a11*a22*a33 *b11*r4102 - a11*a22*a33 *b11 *r4119
2 2 3 2 4 2
+ 4*a11*a22*a33 *b11 *r464 + 3*a11*a33 *b11 *r4119 - 2*a22 *a33 *r446
3 2 3 2 3 2
- 4*a22 *a33 *b11*r4102 - 2*a22 *a33*b11 *r4119 + 2*a22 *a33*b11 *r464
2 3 2 2 2 2 2 2
+ 4*a22 *a33 *b11*r4102 + 5*a22 *a33 *b11 *r4119 - 2*a22 *a33 *b11 *r464
3 2 4 2 3 2
- 3*a22*a33 *b11 *r4119)/(4*a11 *a33 - 16*a11 *a22*a33
2 2 2 3 2 4 2
+ 24*a11 *a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 )
r428=0
r429=0
4 2 4 2 3 2
r430=( - 2*a11 *a33 *r446 - a11 *b11 *r464 + 8*a11 *a22*a33 *r446
3 3 2 3 2
+ 4*a11 *a22*a33*b11*r4102 - a11 *a22*b11 *r4119 - 4*a11 *a33 *b11*r4102
3 2 3 2 2 2 2
+ a11 *a33*b11 *r4119 + 4*a11 *a33*b11 *r464 - 12*a11 *a22 *a33 *r446
2 2 2 2 2
- 8*a11 *a22 *a33*b11*r4102 + 3*a11 *a22 *b11 *r464
2 2 2 2
+ 4*a11 *a22*a33 *b11*r4102 + 3*a11 *a22*a33*b11 *r4119
2 2 2 3
- 6*a11 *a22*a33*b11 *r464 + 4*a11 *a33 *b11*r4102
2 2 2 2 2 2 3 2
- 3*a11 *a33 *b11 *r4119 - 3*a11 *a33 *b11 *r464 + 8*a11*a22 *a33 *r446
3 3 2 3 2
+ 4*a11*a22 *a33*b11*r4102 + a11*a22 *b11 *r4119 - 2*a11*a22 *b11 *r464
2 2 2 2
+ 4*a11*a22 *a33 *b11*r4102 - 3*a11*a22 *a33*b11 *r4119
3 2 2
- 8*a11*a22*a33 *b11*r4102 + 6*a11*a22*a33 *b11 *r464
3 2 4 2 3 2
+ 2*a11*a33 *b11 *r4119 - 2*a22 *a33 *r446 - 4*a22 *a33 *b11*r4102
3 2 3 2 2 3
- a22 *a33*b11 *r4119 + 2*a22 *a33*b11 *r464 + 4*a22 *a33 *b11*r4102
2 2 2 2 2 2 3 2
+ 3*a22 *a33 *b11 *r4119 - 3*a22 *a33 *b11 *r464 - 2*a22*a33 *b11 *r4119)
4 2 3 2 2 2 2 3 2
/(4*a11 *a33 - 16*a11 *a22*a33 + 24*a11 *a22 *a33 - 16*a11*a22 *a33
4 2
+ 4*a22 *a33 )
2 2
r431=(a11 *a33*r4102 - a11 *b11*r464 - a11*a22*a33*r4102 - a11*a22*b11*r4119
2
+ a11*a22*b11*r464 - a11*a33 *r4102 + a11*a33*b11*r4119
2
+ a11*a33*b11*r464 + a22*a33 *r4102 + a22*a33*b11*r4119
2 2 2
- a22*a33*b11*r464 - a33 *b11*r4119)/(a11 *a33 - 2*a11*a22*a33 + a22 *a33
)
r432=0
r433=0
2
r434=(a11 *r464 - a11*a22*r4119 - a11*a22*r464 + a11*a33*r4119 - a11*a33*r464
2 2 2
+ a22*a33*r4119 + a22*a33*r464 - a33 *r4119)/(2*a11 - 4*a11*a22 + 2*a22
)
r435=0
2 3 2 2 3
a11*a33*b11 *r4102 - a11*b11 *r464 - a33 *b11 *r4102 + a33*b11 *r464
r436=----------------------------------------------------------------------
3 3
4*a11*a33 - 4*a22*a33
r437=0
2 2 2 3 2 2
r438=(a11 *a33*b11 *r4102 - a11 *b11 *r464 - 2*a11*a33 *b11 *r4102
3 3 2 2 3 2 3
+ 2*a11*a33*b11 *r464 + a33 *b11 *r4102 - a33 *b11 *r464)/(4*a11 *a33
3 2 3
- 8*a11*a22*a33 + 4*a22 *a33 )
r439=0
r440=0
r441=0
r442=0
r444=0
r445=0
r447=0
r448=0
r449=0
r450=0
r451=0
2
a11*a33*r4102 - a11*b11*r464 - a33 *r4102 + a33*b11*r464
r452=----------------------------------------------------------
a11*a33 - a22*a33
r453=0
r454=0
2 2 2 2
r455=( - 2*a11 *a33 *r446 + 4*a11*a22*a33 *r446 - a11*a22*b11 *r4119
2 2 2
- a11*a22*b11 *r464 + a11*a33*b11 *r4119 + a11*a33*b11 *r464
2 2 2 2 2
- 2*a22 *a33 *r446 + a22 *b11 *r464 + a22*a33*b11 *r4119
2 2 2 2 2 2
- a22*a33*b11 *r464 - a33 *b11 *r4119)/(4*a11 *a33 - 8*a11*a22*a33
2 2
+ 4*a22 *a33 )
r456=0
r458=0
r459=0
2 2 2 2 2
r460=( - 2*a11 *a33 *r446 - a11 *b11 *r464 + 2*a11*a22*a33 *r446
2 3 2
- a11*a22*b11 *r4119 + 2*a11*a33 *r446 + a11*a33*b11 *r4119
2 3 2
+ 2*a11*a33*b11 *r464 - 2*a22*a33 *r446 + a22*a33*b11 *r4119
2 2 2 2 2 2 2
- a33 *b11 *r4119 - a33 *b11 *r464)/(4*a11 *a33 - 4*a11*a22*a33
3 3
- 4*a11*a33 + 4*a22*a33 )
r461=(a11*a33*r4102 - a11*b11*r464 - a22*a33*r4102 - a22*b11*r4119
+ a22*b11*r464 + a33*b11*r4119)/(a11*a33 - a22*a33)
r462=0
r463=0
r465=0
a33*r4102 - b11*r464
r466=----------------------
a33
r467=0
r468=0
2
r469=(a11 *r464 + a11*a22*r4119 - a11*a22*r464 - a11*a33*r4119 - a11*a33*r464
2 2 2
- a22*a33*r4119 + a22*a33*r464 + a33 *r4119)/(2*a11 - 4*a11*a33 + 2*a33
)
r470=0
r471=0
r472=0
r473=0
4 2 4 2 4 3
r474=( - 2*a11 *a33 *b11*r446 + a11 *a33*b11 *r4102 - a11 *b11 *r464
3 2 3 2
+ 8*a11 *a22*a33 *b11*r446 + a11 *a22*a33*b11 *r4102
3 3 3 2 2 3 3
- 2*a11 *a22*b11 *r4119 - 5*a11 *a33 *b11 *r4102 + 2*a11 *a33*b11 *r4119
3 3 2 2 2
+ 4*a11 *a33*b11 *r464 - 12*a11 *a22 *a33 *b11*r446
2 2 2 2 2 3
- 5*a11 *a22 *a33*b11 *r4102 + a11 *a22 *b11 *r4119
2 2 3 2 2 2
+ 3*a11 *a22 *b11 *r464 + 7*a11 *a22*a33 *b11 *r4102
2 3 2 3
+ 4*a11 *a22*a33*b11 *r4119 - 6*a11 *a22*a33*b11 *r464
2 3 2 2 2 3 2 2 3
+ 4*a11 *a33 *b11 *r4102 - 5*a11 *a33 *b11 *r4119 - 3*a11 *a33 *b11 *r464
3 2 3 2
+ 8*a11*a22 *a33 *b11*r446 + 3*a11*a22 *a33*b11 *r4102
3 3 3 3 2 2 2
+ a11*a22 *b11 *r4119 - 2*a11*a22 *b11 *r464 + a11*a22 *a33 *b11 *r4102
2 3 3 2
- 5*a11*a22 *a33*b11 *r4119 - 8*a11*a22*a33 *b11 *r4102
2 3 2 3
+ a11*a22*a33 *b11 *r4119 + 6*a11*a22*a33 *b11 *r464
3 3 4 2 3 2 2
+ 3*a11*a33 *b11 *r4119 - 2*a22 *a33 *b11*r446 - 3*a22 *a33 *b11 *r4102
3 3 3 3 2 3 2
- a22 *a33*b11 *r4119 + 2*a22 *a33*b11 *r464 + 4*a22 *a33 *b11 *r4102
2 2 3 2 2 3 3 3
+ 4*a22 *a33 *b11 *r4119 - 3*a22 *a33 *b11 *r464 - 3*a22*a33 *b11 *r4119)
4 3 3 3 2 2 3 3 3
/(4*a11 *a33 - 16*a11 *a22*a33 + 24*a11 *a22 *a33 - 16*a11*a22 *a33
4 3
+ 4*a22 *a33 )
r475=0
4 2 4 2 4 3
r476=( - 2*a11 *a33 *b11*r446 + a11 *a33*b11 *r4102 - a11 *b11 *r464
3 2 3 2
+ 8*a11 *a22*a33 *b11*r446 + 2*a11 *a22*a33*b11 *r4102
3 3 3 2 2 3 3
- a11 *a22*b11 *r4119 - 6*a11 *a33 *b11 *r4102 + a11 *a33*b11 *r4119
3 3 2 2 2
+ 4*a11 *a33*b11 *r464 - 12*a11 *a22 *a33 *b11*r446
2 2 2 2 2 3
- 7*a11 *a22 *a33*b11 *r4102 - a11 *a22 *b11 *r4119
2 2 3 2 2 2
+ 3*a11 *a22 *b11 *r464 + 8*a11 *a22*a33 *b11 *r4102
2 3 2 3
+ 5*a11 *a22*a33*b11 *r4119 - 6*a11 *a22*a33*b11 *r464
2 3 2 2 2 3 2 2 3
+ 5*a11 *a33 *b11 *r4102 - 4*a11 *a33 *b11 *r4119 - 3*a11 *a33 *b11 *r464
3 2 3 2
+ 8*a11*a22 *a33 *b11*r446 + 4*a11*a22 *a33*b11 *r4102
3 3 3 3
+ 2*a11*a22 *b11 *r4119 - 2*a11*a22 *b11 *r464
2 2 2 2 3
+ 2*a11*a22 *a33 *b11 *r4102 - 4*a11*a22 *a33*b11 *r4119
3 2 2 3
- 10*a11*a22*a33 *b11 *r4102 - a11*a22*a33 *b11 *r4119
2 3 3 3 4 2
+ 6*a11*a22*a33 *b11 *r464 + 3*a11*a33 *b11 *r4119 - 2*a22 *a33 *b11*r446
3 2 2 3 3 3 3
- 4*a22 *a33 *b11 *r4102 - 2*a22 *a33*b11 *r4119 + 2*a22 *a33*b11 *r464
2 3 2 2 2 3 2 2 3
+ 5*a22 *a33 *b11 *r4102 + 5*a22 *a33 *b11 *r4119 - 3*a22 *a33 *b11 *r464
3 3 4 3 3 3
- 3*a22*a33 *b11 *r4119)/(4*a11 *a33 - 16*a11 *a22*a33
2 2 3 3 3 4 3
+ 24*a11 *a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 )
r477=0
r478=0
4 2 4 3 3 2
r479=( - 2*a11 *a33 *b11*r446 - a11 *b11 *r464 + 8*a11 *a22*a33 *b11*r446
3 2 3 3
+ 4*a11 *a22*a33*b11 *r4102 - a11 *a22*b11 *r4119
3 2 2 3 3 3 3
- 4*a11 *a33 *b11 *r4102 + a11 *a33*b11 *r4119 + 4*a11 *a33*b11 *r464
2 2 2 2 2 2
- 12*a11 *a22 *a33 *b11*r446 - 8*a11 *a22 *a33*b11 *r4102
2 2 3 2 2 2
+ 3*a11 *a22 *b11 *r464 + 4*a11 *a22*a33 *b11 *r4102
2 3 2 3
+ 3*a11 *a22*a33*b11 *r4119 - 6*a11 *a22*a33*b11 *r464
2 3 2 2 2 3 2 2 3
+ 4*a11 *a33 *b11 *r4102 - 3*a11 *a33 *b11 *r4119 - 3*a11 *a33 *b11 *r464
3 2 3 2
+ 8*a11*a22 *a33 *b11*r446 + 4*a11*a22 *a33*b11 *r4102
3 3 3 3 2 2 2
+ a11*a22 *b11 *r4119 - 2*a11*a22 *b11 *r464 + 4*a11*a22 *a33 *b11 *r4102
2 3 3 2
- 3*a11*a22 *a33*b11 *r4119 - 8*a11*a22*a33 *b11 *r4102
2 3 3 3 4 2
+ 6*a11*a22*a33 *b11 *r464 + 2*a11*a33 *b11 *r4119 - 2*a22 *a33 *b11*r446
3 2 2 3 3 3 3
- 4*a22 *a33 *b11 *r4102 - a22 *a33*b11 *r4119 + 2*a22 *a33*b11 *r464
2 3 2 2 2 3 2 2 3
+ 4*a22 *a33 *b11 *r4102 + 3*a22 *a33 *b11 *r4119 - 3*a22 *a33 *b11 *r464
3 3 4 3 3 3
- 2*a22*a33 *b11 *r4119)/(4*a11 *a33 - 16*a11 *a22*a33
2 2 3 3 3 4 3
+ 24*a11 *a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 )
r480=0
r481=0
r482=0
r484=0
r485=0
r486=0
r487=0
2 2
r488=(a11 *a33*r4102 - a11*a22*a33*r4102 - a11*a22*b11*r4119 - a11*a33 *r4102
2 2
+ a11*a33*b11*r4119 + a22*a33 *r4102 + a22*a33*b11*r4119 - a33 *b11*r4119
2 2
)/(a11 *a33 - 2*a11*a22*a33 + a22 *a33)
r489=0
r490=0
r491=0
r492=0
r493=0
r495=0
r496=0
r497=0
r498=0
r499=0
r4100=0
r4101=0
r4103=0
r4104=0
2 2 2 2 2
r4105=( - 2*a11 *a33 *r446 - 4*a11 *a33*b11*r4102 + a11 *b11 *r4119
2 2 2
+ 2*a11 *b11 *r464 + 4*a11*a22*a33 *r446 + 4*a11*a22*a33*b11*r4102
2 2 2
+ a11*a22*b11 *r4119 - 2*a11*a22*b11 *r464 + 4*a11*a33 *b11*r4102
2 2 2 2
- 3*a11*a33*b11 *r4119 - 2*a11*a33*b11 *r464 - 2*a22 *a33 *r446
2 2 2
- 4*a22*a33 *b11*r4102 - a22*a33*b11 *r4119 + 2*a22*a33*b11 *r464
2 2 2 2 2 2 2
+ 2*a33 *b11 *r4119)/(4*a11 *a33 - 8*a11*a22*a33 + 4*a22 *a33 )
r4106=0
2 2 2 2 2
r4107=( - 2*a11 *a33 *r446 - 4*a11 *a33*b11*r4102 + a11 *b11 *r4119
2 2 2
+ 2*a11 *b11 *r464 + 4*a11*a22*a33 *r446 + 4*a11*a22*a33*b11*r4102
2 2 2
+ 2*a11*a22*b11 *r4119 - 2*a11*a22*b11 *r464 + 4*a11*a33 *b11*r4102
2 2 2 2
- 4*a11*a33*b11 *r4119 - 2*a11*a33*b11 *r464 - 2*a22 *a33 *r446
2 2 2
- 4*a22*a33 *b11*r4102 - 2*a22*a33*b11 *r4119 + 2*a22*a33*b11 *r464
2 2 2 2 2 2 2
+ 3*a33 *b11 *r4119)/(4*a11 *a33 - 8*a11*a22*a33 + 4*a22 *a33 )
r4108=0
r4109=0
r4111=0
r4112=0
r4113=0
a11*r4119 - a33*r4119
r4114=-----------------------
a11 - a22
r4115=0
r4117=0
r4118=0
r4120=0
r4121=0
a11*b11*r4119 - a33*b11*r4119
r4122=-------------------------------
a11*a33 - a22*a33
r4123=0
r4124=0
m3=0
m2=0
m1=0
n3=0
n2=0
n1=0
2 2 2 2 2 2 2 2
a11 *a33*b11 - a11*a22 *b11 - a11*a33 *b11 + a22 *a33*b11
c33=---------------------------------------------------------------
2 2 2 2 2
4*a11 *a33 - 8*a11*a22*a33 + 4*a22 *a33
c23=0
2 2 2
a11*a22*b11 - a33 *b11
c22=--------------------------
2 2
4*a11*a33 - 4*a22*a33
c13=0
c12=0
c11=0
- a11*a22*b11 + a11*a33*b11
b33=------------------------------
a11*a33 - a22*a33
b32=0
b31=0
b23=0
b22=0
b21=0
b13=0
b12=0
a23=0
a13=0
a12=0
2 2 2
r494=(a11*a33 *r446 + a11*a33*b11*r4102 - a11*b11 *r464 - a22*a33 *r446
2 2 2 2
- a33 *b11*r4102 + a33*b11 *r464)/(a11*a33 - a22*a33 )
2 2 2 2 2 2
r483=(a11 *a33 *r446 + a11 *a33*b11*r4102 - a11 *b11 *r464 - 2*a11*a22*a33 *r446
2 2
- a11*a22*a33*b11*r4102 - a11*a22*b11 *r4119 + a11*a22*b11 *r464
2 2 2
- a11*a33 *b11*r4102 + a11*a33*b11 *r4119 + a11*a33*b11 *r464
2 2 2 2
+ a22 *a33 *r446 + a22*a33 *b11*r4102 + a22*a33*b11 *r4119
2 2 2 2 2 2
- a22*a33*b11 *r464 - a33 *b11 *r4119)/(a11 *a33 - 2*a11*a22*a33
2 2
+ a22 *a33 )
b11*r29
r217=---------
a33
Parameters
Apart from the condition that they must not vanish to give
a non-trivial solution and a non-singular solution with
non-vanishing denominators, the following parameters are free:
r29, r4119, r4102, r446, r464, b11, a33, a11, a22
Inequalities
In the following not identically vanishing expressions are shown.
Any auxiliary variables g00?? are used to express that at least
one of their coefficients must not vanish, e.g. g0019*p4 + g0020*p3
means that either p4 or p3 or both are non-vanishing.
{a11*a22 - a11*a33 - 3*a22*a33,
a11*a22 + a11*a33 - a22*a33,
2
a11 - ---*a22,
3
a11,
a22,
a33,
a11 - a22,
a22 - a33,
a11 - a33,
a11*a22 - a11*a33 + a22*a33,
2
a11*a22 - a33 }
Relevance for the application:
The new Hamiltonian in form of a list of vanishing expressions:
{a12,
a13,
a23,
b12,
b13,
b21,
b22,
b23,
b31,
b32,
a11*a22*b11 - a11*a33*b11 + a11*a33*b33 - a22*a33*b33,
c11,
c12,
c13,
- a11*a22*b11**2 + 4*a11*a33**2*c22 - 4*a22*a33**2*c22 + a33**2*b11**2,
c23,
4*a11**2*a33**2*c33 - a11**2*a33*b11**2 + a11*a22**2*b11**2 - 8*a11*a22*a33**2*
c33 + a11*a33**2*b11**2 + 4*a22**2*a33**2*c33 - a22**2*a33*b11**2,
n1,
n2,
n3,
m1,
m2,
m3}$
The system of equations related to the Hamiltonian HAM:
2 2 2
HAM=u1 *a11 + u1*v1*b11 + u2 *a22 + u3 *a33
2 2 2
- a11*a22*b11 + a11*a33*b11 2 a11*a22*b11 - a33 *b11
+ u3*v3*------------------------------ + v2 *--------------------------
a11*a33 - a22*a33 2 2
4*a11*a33 - 4*a22*a33
2 2 2 2 2 2 2 2
2 a11 *a33*b11 - a11*a22 *b11 - a11*a33 *b11 + a22 *a33*b11
+ v3 *---------------------------------------------------------------
2 2 2 2 2
4*a11 *a33 - 8*a11*a22*a33 + 4*a22 *a33
has apart from the Hamiltonian and Casimirs the following 5 first integrals:
2 2 5 2 2 4 2 2 4 3 2
FI=u1 *v2 *(16*a11 *a33 *b11 - 48*a11 *a22*a33 *b11 - 32*a11 *a33 *b11
3 2 2 2 3 3 2
+ 48*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11
3 4 2 2 3 2 2
+ 16*a11 *a33 *b11 - 16*a11 *a22 *a33 *b11
2 2 3 2 2 4 2
- 96*a11 *a22 *a33 *b11 - 48*a11 *a22*a33 *b11
3 3 2 2 4 2 3 4 2
+ 32*a11*a22 *a33 *b11 + 48*a11*a22 *a33 *b11 - 16*a22 *a33 *b11
2 2 5 2 2 4 2 2
) + u1 *v3 *(16*a11 *a33 *b11 - 48*a11 *a22*a33 *b11
4 3 2 3 2 2 2 3 3 2
- 32*a11 *a33 *b11 + 48*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11
3 4 2 2 3 2 2 2 2 3 2
+ 16*a11 *a33 *b11 - 16*a11 *a22 *a33 *b11 - 96*a11 *a22 *a33 *b11
2 4 2 3 3 2 2 4 2
- 48*a11 *a22*a33 *b11 + 32*a11*a22 *a33 *b11 + 48*a11*a22 *a33 *b11
3 4 2 5 2 2
- 16*a22 *a33 *b11 ) + u1*u2*v1*v2*( - 32*a11 *a33 *b11
4 2 2 4 3 2 3 2 2 2
+ 96*a11 *a22*a33 *b11 + 64*a11 *a33 *b11 - 96*a11 *a22 *a33 *b11
3 3 2 3 4 2 2 3 2 2
- 192*a11 *a22*a33 *b11 - 32*a11 *a33 *b11 + 32*a11 *a22 *a33 *b11
2 2 3 2 2 4 2 3 3 2
+ 192*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11 - 64*a11*a22 *a33 *b11
2 4 2 3 4 2
- 96*a11*a22 *a33 *b11 + 32*a22 *a33 *b11 ) + u1*u3*v1*v3*(
5 2 2 4 2 2 4 3 2
- 32*a11 *a33 *b11 + 96*a11 *a22*a33 *b11 + 64*a11 *a33 *b11
3 2 2 2 3 3 2 3 4 2
- 96*a11 *a22 *a33 *b11 - 192*a11 *a22*a33 *b11 - 32*a11 *a33 *b11
2 3 2 2 2 2 3 2
+ 32*a11 *a22 *a33 *b11 + 192*a11 *a22 *a33 *b11
2 4 2 3 3 2 2 4 2
+ 96*a11 *a22*a33 *b11 - 64*a11*a22 *a33 *b11 - 96*a11*a22 *a33 *b11
3 4 2 3 5 3 4 2 3
+ 32*a22 *a33 *b11 ) + u1*v1 *( - 8*a11 *a33*b11 + 40*a11 *a33 *b11
3 2 3 3 2 3 3 3 3
+ 24*a11 *a22 *a33*b11 - 48*a11 *a22*a33 *b11 - 56*a11 *a33 *b11
2 3 3 2 2 2 3 2 3 3
- 16*a11 *a22 *a33*b11 - 24*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11
2 4 3 3 2 3 2 3 3
+ 24*a11 *a33 *b11 + 32*a11*a22 *a33 *b11 - 24*a11*a22 *a33 *b11
4 3 3 3 3 2 4 3
- 48*a11*a22*a33 *b11 - 16*a22 *a33 *b11 + 24*a22 *a33 *b11 ) + u1*v1
2 5 3 4 2 3 3 2 3
*v2 *( - 8*a11 *a33*b11 + 40*a11 *a33 *b11 + 24*a11 *a22 *a33*b11
3 2 3 3 3 3 2 3 3
- 48*a11 *a22*a33 *b11 - 56*a11 *a33 *b11 - 16*a11 *a22 *a33*b11
2 2 2 3 2 3 3 2 4 3
- 24*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11 + 24*a11 *a33 *b11
3 2 3 2 3 3 4 3
+ 32*a11*a22 *a33 *b11 - 24*a11*a22 *a33 *b11 - 48*a11*a22*a33 *b11
3 3 3 2 4 3 2
- 16*a22 *a33 *b11 + 24*a22 *a33 *b11 ) + u1*v1*v3 *(
5 3 4 2 3 3 2 3
- 8*a11 *a33*b11 + 40*a11 *a33 *b11 + 24*a11 *a22 *a33*b11
3 2 3 3 3 3 2 3 3
- 48*a11 *a22*a33 *b11 - 56*a11 *a33 *b11 - 16*a11 *a22 *a33*b11
2 2 2 3 2 3 3 2 4 3
- 24*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11 + 24*a11 *a33 *b11
3 2 3 2 3 3 4 3
+ 32*a11*a22 *a33 *b11 - 24*a11*a22 *a33 *b11 - 48*a11*a22*a33 *b11
3 3 3 2 4 3 4 5 4
- 16*a22 *a33 *b11 + 24*a22 *a33 *b11 ) + u2 *(16*a11 *a33
4 4 3 2 4 2 3 4
- 80*a11 *a22*a33 + 160*a11 *a22 *a33 - 160*a11 *a22 *a33
4 4 5 4 3 5 3
+ 80*a11*a22 *a33 - 16*a22 *a33 ) + u2 *v2*( - 32*a11 *a33 *b11
4 3 4 4 3 2 3
+ 128*a11 *a22*a33 *b11 + 32*a11 *a33 *b11 - 192*a11 *a22 *a33 *b11
3 4 2 3 3 2 2 4
- 128*a11 *a22*a33 *b11 + 128*a11 *a22 *a33 *b11 + 192*a11 *a22 *a33 *b11
4 3 3 4 4 4 2
- 32*a11*a22 *a33 *b11 - 128*a11*a22 *a33 *b11 + 32*a22 *a33 *b11) + u2
2 5 4 4 4 4 5 3 2 4
*u3 *(32*a11 *a33 - 128*a11 *a22*a33 - 32*a11 *a33 + 192*a11 *a22 *a33
3 5 2 3 4 2 2 5
+ 128*a11 *a22*a33 - 128*a11 *a22 *a33 - 192*a11 *a22 *a33
4 4 3 5 4 5 2
+ 32*a11*a22 *a33 + 128*a11*a22 *a33 - 32*a22 *a33 ) + u2 *u3*v3*(
5 3 4 3 4 4
- 32*a11 *a33 *b11 + 128*a11 *a22*a33 *b11 + 32*a11 *a33 *b11
3 2 3 3 4 2 3 3
- 192*a11 *a22 *a33 *b11 - 128*a11 *a22*a33 *b11 + 128*a11 *a22 *a33 *b11
2 2 4 4 3 3 4
+ 192*a11 *a22 *a33 *b11 - 32*a11*a22 *a33 *b11 - 128*a11*a22 *a33 *b11
4 4 2 2 5 2 2
+ 32*a22 *a33 *b11) + u2 *v1 *( - 8*a11 *a33 *b11
4 2 2 4 3 2 3 2 2 2
+ 24*a11 *a22*a33 *b11 + 16*a11 *a33 *b11 - 24*a11 *a22 *a33 *b11
3 3 2 3 4 2 2 3 2 2
- 48*a11 *a22*a33 *b11 - 8*a11 *a33 *b11 + 8*a11 *a22 *a33 *b11
2 2 3 2 2 4 2 3 3 2
+ 48*a11 *a22 *a33 *b11 + 24*a11 *a22*a33 *b11 - 16*a11*a22 *a33 *b11
2 4 2 3 4 2 2 2
- 24*a11*a22 *a33 *b11 + 8*a22 *a33 *b11 ) + u2 *v3 *(
4 2 2 4 3 2 3 2 2 2
- 8*a11 *a22*a33 *b11 + 8*a11 *a33 *b11 + 24*a11 *a22 *a33 *b11
3 3 2 3 4 2 2 3 2 2
- 16*a11 *a22*a33 *b11 - 8*a11 *a33 *b11 - 24*a11 *a22 *a33 *b11
2 4 2 4 2 2 3 3 2
+ 24*a11 *a22*a33 *b11 + 8*a11*a22 *a33 *b11 + 16*a11*a22 *a33 *b11
2 4 2 4 3 2 3 4 2 2
- 24*a11*a22 *a33 *b11 - 8*a22 *a33 *b11 + 8*a22 *a33 *b11 ) + u2*u3
5 3 4 3 4 4
*v2*( - 32*a11 *a33 *b11 + 96*a11 *a22*a33 *b11 + 64*a11 *a33 *b11
3 2 3 3 4 3 5
- 96*a11 *a22 *a33 *b11 - 192*a11 *a22*a33 *b11 - 32*a11 *a33 *b11
2 3 3 2 2 4 2 5
+ 32*a11 *a22 *a33 *b11 + 192*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11
3 4 2 5 3 5
- 64*a11*a22 *a33 *b11 - 96*a11*a22 *a33 *b11 + 32*a22 *a33 *b11) + u2
3 5 3 4 3 4 2 3
*v2 *( - 8*a11 *a33*b11 + 16*a11 *a22*a33*b11 + 24*a11 *a33 *b11
3 2 3 3 2 3 3 3 3
- 8*a11 *a22 *a33*b11 - 48*a11 *a22*a33 *b11 - 24*a11 *a33 *b11
2 2 2 3 2 3 3 2 4 3
+ 24*a11 *a22 *a33 *b11 + 48*a11 *a22*a33 *b11 + 8*a11 *a33 *b11
2 3 3 4 3 2 4 3
- 24*a11*a22 *a33 *b11 - 16*a11*a22*a33 *b11 + 8*a22 *a33 *b11 ) +
2 5 3 4 3 4 2 3
u2*v2*v3 *( - 8*a11 *a33*b11 + 24*a11 *a22*a33*b11 + 16*a11 *a33 *b11
3 2 3 3 2 3
- 24*a11 *a22 *a33*b11 - 48*a11 *a22*a33 *b11
3 3 3 2 3 3
- 8*a11 *a33 *b11 + 8*a11 *a22 *a33*b11
2 2 2 3 2 3 3
+ 48*a11 *a22 *a33 *b11 + 24*a11 *a22*a33 *b11
3 2 3 2 3 3
- 16*a11*a22 *a33 *b11 - 24*a11*a22 *a33 *b11
3 3 3 4 5 4 4 4
+ 8*a22 *a33 *b11 ) + u3 *(16*a11 *a33 - 48*a11 *a22*a33
4 5 3 2 4 3 5 3 6
- 32*a11 *a33 + 48*a11 *a22 *a33 + 96*a11 *a22*a33 + 16*a11 *a33
2 3 4 2 2 5 2 6
- 16*a11 *a22 *a33 - 96*a11 *a22 *a33 - 48*a11 *a22*a33
3 5 2 6 3 6 3
+ 32*a11*a22 *a33 + 48*a11*a22 *a33 - 16*a22 *a33 ) + u3 *v3*(
5 3 4 3 4 4
- 32*a11 *a33 *b11 + 96*a11 *a22*a33 *b11 + 64*a11 *a33 *b11
3 2 3 3 4 3 5
- 96*a11 *a22 *a33 *b11 - 192*a11 *a22*a33 *b11 - 32*a11 *a33 *b11
2 3 3 2 2 4 2 5
+ 32*a11 *a22 *a33 *b11 + 192*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11
3 4 2 5 3 5 2
- 64*a11*a22 *a33 *b11 - 96*a11*a22 *a33 *b11 + 32*a22 *a33 *b11) + u3
2 5 2 2 4 3 2 3 2 2 2
*v1 *( - 8*a11 *a33 *b11 + 40*a11 *a33 *b11 + 24*a11 *a22 *a33 *b11
3 3 2 3 4 2 2 3 2 2
- 48*a11 *a22*a33 *b11 - 56*a11 *a33 *b11 - 16*a11 *a22 *a33 *b11
2 2 3 2 2 4 2 2 5 2
- 24*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11 + 24*a11 *a33 *b11
3 3 2 2 4 2 5 2
+ 32*a11*a22 *a33 *b11 - 24*a11*a22 *a33 *b11 - 48*a11*a22*a33 *b11
3 4 2 2 5 2 2 2
- 16*a22 *a33 *b11 + 24*a22 *a33 *b11 ) + u3 *v2 *(
4 2 2 4 3 2 3 2 2 2
- 16*a11 *a22*a33 *b11 + 16*a11 *a33 *b11 + 32*a11 *a22 *a33 *b11
3 4 2 2 3 2 2 2 2 3 2
- 32*a11 *a33 *b11 - 16*a11 *a22 *a33 *b11 - 48*a11 *a22 *a33 *b11
2 4 2 2 5 2 3 3 2
+ 48*a11 *a22*a33 *b11 + 16*a11 *a33 *b11 + 32*a11*a22 *a33 *b11
5 2 3 4 2 2 5 2 2 2
- 32*a11*a22*a33 *b11 - 16*a22 *a33 *b11 + 16*a22 *a33 *b11 ) + u3 *v3
4 2 2 4 3 2 3 2 2 2
*( - 24*a11 *a22*a33 *b11 + 24*a11 *a33 *b11 + 48*a11 *a22 *a33 *b11
3 4 2 2 3 2 2 2 2 3 2
- 48*a11 *a33 *b11 - 24*a11 *a22 *a33 *b11 - 72*a11 *a22 *a33 *b11
2 4 2 2 5 2 3 3 2
+ 72*a11 *a22*a33 *b11 + 24*a11 *a33 *b11 + 48*a11*a22 *a33 *b11
5 2 3 4 2 2 5 2 2
- 48*a11*a22*a33 *b11 - 24*a22 *a33 *b11 + 24*a22 *a33 *b11 ) + u3*v1
4 3 4 2 3 3 2 3
*v3*( - 16*a11 *a22*a33*b11 + 16*a11 *a33 *b11 + 32*a11 *a22 *a33*b11
3 3 3 2 3 3 2 2 2 3
- 32*a11 *a33 *b11 - 16*a11 *a22 *a33*b11 - 48*a11 *a22 *a33 *b11
2 3 3 2 4 3 3 2 3
+ 48*a11 *a22*a33 *b11 + 16*a11 *a33 *b11 + 32*a11*a22 *a33 *b11
4 3 3 3 3 2 4 3
- 32*a11*a22*a33 *b11 - 16*a22 *a33 *b11 + 16*a22 *a33 *b11 ) + u3
2 5 3 4 2 3 3 2 3
*v2 *v3*( - 8*a11 *a33*b11 + 40*a11 *a33 *b11 + 24*a11 *a22 *a33*b11
3 2 3 3 3 3 2 3 3
- 48*a11 *a22*a33 *b11 - 56*a11 *a33 *b11 - 16*a11 *a22 *a33*b11
2 2 2 3 2 3 3
- 24*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11
2 4 3 3 2 3 2 3 3
+ 24*a11 *a33 *b11 + 32*a11*a22 *a33 *b11 - 24*a11*a22 *a33 *b11
4 3 3 3 3 2 4 3
- 48*a11*a22*a33 *b11 - 16*a22 *a33 *b11 + 24*a22 *a33 *b11 ) +
3 5 3 4 3 4 2 3
u3*v3 *( - 8*a11 *a33*b11 + 8*a11 *a22*a33*b11 + 32*a11 *a33 *b11
3 2 3 3 2 3 3 3 3
+ 8*a11 *a22 *a33*b11 - 48*a11 *a22*a33 *b11 - 40*a11 *a33 *b11
2 3 3 2 3 3 2 4 3
- 8*a11 *a22 *a33*b11 + 72*a11 *a22*a33 *b11 + 16*a11 *a33 *b11
3 2 3 2 3 3
+ 16*a11*a22 *a33 *b11 - 24*a11*a22 *a33 *b11
4 3 3 3 3 2 4 3 2
- 32*a11*a22*a33 *b11 - 8*a22 *a33 *b11 + 16*a22 *a33 *b11 ) + v1
2 5 4 4 4 4 4 3 2 4
*v2 *( - 2*a11 *b11 + 2*a11 *a22*b11 + 8*a11 *a33*b11 - 4*a11 *a22 *b11
3 2 4 2 3 4 2 3 4
- 16*a11 *a33 *b11 + 4*a11 *a22 *b11 + 16*a11 *a33 *b11
3 4 2 2 4 3 4
- 8*a11*a22 *a33*b11 + 12*a11*a22 *a33 *b11 - 8*a11*a22*a33 *b11
4 4 3 2 4 2 3 4
- 6*a11*a33 *b11 + 4*a22 *a33 *b11 - 8*a22 *a33 *b11
4 4 2 2 5 4 4 4
+ 6*a22*a33 *b11 ) + v1 *v3 *( - 2*a11 *b11 + 10*a11 *a33*b11
3 2 4 3 4 3 2 4
+ 2*a11 *a22 *b11 - 4*a11 *a22*a33*b11 - 18*a11 *a33 *b11
2 2 4 2 2 4 2 3 4
- 6*a11 *a22 *a33*b11 + 12*a11 *a22*a33 *b11 + 14*a11 *a33 *b11
2 2 4 3 4 4 4
+ 6*a11*a22 *a33 *b11 - 12*a11*a22*a33 *b11 - 4*a11*a33 *b11
2 3 4 4 4 4 5 4 4 4
- 2*a22 *a33 *b11 + 4*a22*a33 *b11 ) + v2 *( - a11 *b11 + a11 *a22*b11
4 4 3 2 4 3 4
+ 4*a11 *a33*b11 - 4*a11 *a22 *b11 + 4*a11 *a22*a33*b11
3 2 4 2 3 4 2 2 4
- 10*a11 *a33 *b11 + 4*a11 *a22 *b11 - 6*a11 *a22*a33 *b11
2 3 4 3 4 2 2 4
+ 12*a11 *a33 *b11 - 8*a11*a22 *a33*b11 + 12*a11*a22 *a33 *b11
3 4 4 4 3 2 4
- 4*a11*a22*a33 *b11 - 5*a11*a33 *b11 + 4*a22 *a33 *b11
2 3 4 4 4 2 2 5 4
- 8*a22 *a33 *b11 + 5*a22*a33 *b11 ) + v2 *v3 *( - 2*a11 *b11
4 4 4 4 3 4
- 2*a11 *a22*b11 + 12*a11 *a33*b11 + 8*a11 *a22*a33*b11
3 2 4 2 3 4 2 2 4
- 28*a11 *a33 *b11 + 4*a11 *a22 *b11 - 12*a11 *a22 *a33*b11
2 3 4 3 4 2 2 4
+ 28*a11 *a33 *b11 - 8*a11*a22 *a33*b11 + 24*a11*a22 *a33 *b11
3 4 4 4 3 2 4
- 16*a11*a22*a33 *b11 - 10*a11*a33 *b11 + 4*a22 *a33 *b11
2 3 4 4 4 4 5 4
- 12*a22 *a33 *b11 + 10*a22*a33 *b11 ) + v3 *( - a11 *b11
4 4 4 4 3 2 4
- 3*a11 *a22*b11 + 8*a11 *a33*b11 + 5*a11 *a22 *b11
3 4 3 2 4 2 3 4
+ 2*a11 *a22*a33*b11 - 17*a11 *a33 *b11 - a11 *a22 *b11
2 2 4 2 2 4 2 3 4
- 12*a11 *a22 *a33*b11 + 9*a11 *a22*a33 *b11 + 14*a11 *a33 *b11
3 4 2 2 4 3 4
+ 2*a11*a22 *a33*b11 + 9*a11*a22 *a33 *b11 - 12*a11*a22*a33 *b11
4 4 3 2 4 2 3 4 4 4
- 4*a11*a33 *b11 - a22 *a33 *b11 - 2*a22 *a33 *b11 + 4*a22*a33 *b11 )
= a product of the elements of: {16,
a11 - a22,
2 2 4 2 2 3 2 2 3 3 2
u1 *v2 *(a11 *a33 *b11 - 2*a11 *a22*a33 *b11 - 2*a11 *a33 *b11
2 2 2 2 2 3 2 2 4 2
+ a11 *a22 *a33 *b11 + 4*a11 *a22*a33 *b11 + a11 *a33 *b11
2 3 2 4 2 2 4 2
- 2*a11*a22 *a33 *b11 - 2*a11*a22*a33 *b11 + a22 *a33 *b11 ) +
2 2 4 2 2 3 2 2 3 3 2
u1 *v3 *(a11 *a33 *b11 - 2*a11 *a22*a33 *b11 - 2*a11 *a33 *b11
2 2 2 2 2 3 2 2 4 2
+ a11 *a22 *a33 *b11 + 4*a11 *a22*a33 *b11 + a11 *a33 *b11
2 3 2 4 2 2 4 2
- 2*a11*a22 *a33 *b11 - 2*a11*a22*a33 *b11 + a22 *a33 *b11 ) + u1
4 2 2 3 2 2 3 3 2
*u2*v1*v2*( - 2*a11 *a33 *b11 + 4*a11 *a22*a33 *b11 + 4*a11 *a33 *b11
2 2 2 2 2 3 2 2 4 2
- 2*a11 *a22 *a33 *b11 - 8*a11 *a22*a33 *b11 - 2*a11 *a33 *b11
2 3 2 4 2 2 4 2
+ 4*a11*a22 *a33 *b11 + 4*a11*a22*a33 *b11 - 2*a22 *a33 *b11 )
4 2 2 3 2 2 3 3 2
+ u1*u3*v1*v3*( - 2*a11 *a33 *b11 + 4*a11 *a22*a33 *b11 + 4*a11 *a33 *b11
2 2 2 2 2 3 2 2 4 2
- 2*a11 *a22 *a33 *b11 - 8*a11 *a22*a33 *b11 - 2*a11 *a33 *b11
2 3 2 4 2 2 4 2 3
+ 4*a11*a22 *a33 *b11 + 4*a11*a22*a33 *b11 - 2*a22 *a33 *b11 ) + u1*v1
4 3 3 3 3 2 3
*( - a11 *a33*b11 - a11 *a22*a33*b11 + 5*a11 *a33 *b11
2 2 3 2 2 3 2 3 3
+ 2*a11 *a22 *a33*b11 - a11 *a22*a33 *b11 - 7*a11 *a33 *b11
2 2 3 3 3 4 3
- 4*a11*a22 *a33 *b11 + 5*a11*a22*a33 *b11 + 3*a11*a33 *b11
2 3 3 4 3 2 4 3
+ 2*a22 *a33 *b11 - 3*a22*a33 *b11 )/2 + u1*v1*v2 *( - a11 *a33*b11
3 3 3 2 3 2 2 3
- a11 *a22*a33*b11 + 5*a11 *a33 *b11 + 2*a11 *a22 *a33*b11
2 2 3 2 3 3 2 2 3
- a11 *a22*a33 *b11 - 7*a11 *a33 *b11 - 4*a11*a22 *a33 *b11
3 3 4 3 2 3 3
+ 5*a11*a22*a33 *b11 + 3*a11*a33 *b11 + 2*a22 *a33 *b11
4 3 2 4 3 3 3
- 3*a22*a33 *b11 )/2 + u1*v1*v3 *( - a11 *a33*b11 - a11 *a22*a33*b11
3 2 3 2 2 3 2 2 3
+ 5*a11 *a33 *b11 + 2*a11 *a22 *a33*b11 - a11 *a22*a33 *b11
2 3 3 2 2 3 3 3
- 7*a11 *a33 *b11 - 4*a11*a22 *a33 *b11 + 5*a11*a22*a33 *b11
4 3 2 3 3 4 3 4
+ 3*a11*a33 *b11 + 2*a22 *a33 *b11 - 3*a22*a33 *b11 )/2 + u2 *(
4 4 3 4 2 2 4 3 4
a11 *a33 - 4*a11 *a22*a33 + 6*a11 *a22 *a33 - 4*a11*a22 *a33
4 4 3 4 3 3 3
+ a22 *a33 ) + u2 *v2*( - 2*a11 *a33 *b11 + 6*a11 *a22*a33 *b11
3 4 2 2 3 2 4
+ 2*a11 *a33 *b11 - 6*a11 *a22 *a33 *b11 - 6*a11 *a22*a33 *b11
3 3 2 4 3 4 2 2
+ 2*a11*a22 *a33 *b11 + 6*a11*a22 *a33 *b11 - 2*a22 *a33 *b11) + u2 *u3 *
4 4 3 4 3 5 2 2 4
(2*a11 *a33 - 6*a11 *a22*a33 - 2*a11 *a33 + 6*a11 *a22 *a33
2 5 3 4 2 5 3 5 2
+ 6*a11 *a22*a33 - 2*a11*a22 *a33 - 6*a11*a22 *a33 + 2*a22 *a33 ) + u2
4 3 3 3 3 4
*u3*v3*( - 2*a11 *a33 *b11 + 6*a11 *a22*a33 *b11 + 2*a11 *a33 *b11
2 2 3 2 4 3 3
- 6*a11 *a22 *a33 *b11 - 6*a11 *a22*a33 *b11 + 2*a11*a22 *a33 *b11
2 4 3 4 2 2
+ 6*a11*a22 *a33 *b11 - 2*a22 *a33 *b11) + u2 *v1 *(
4 2 2 3 2 2 3 3 2
- a11 *a33 *b11 + 2*a11 *a22*a33 *b11 + 2*a11 *a33 *b11
2 2 2 2 2 3 2 2 4 2
- a11 *a22 *a33 *b11 - 4*a11 *a22*a33 *b11 - a11 *a33 *b11
2 3 2 4 2 2 4 2 2
+ 2*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 - a22 *a33 *b11 )/2 + u2
2 3 2 2 3 3 2 2 2 2 2
*v3 *( - a11 *a22*a33 *b11 + a11 *a33 *b11 + 2*a11 *a22 *a33 *b11
2 3 2 2 4 2 3 2 2
- a11 *a22*a33 *b11 - a11 *a33 *b11 - a11*a22 *a33 *b11
2 3 2 4 2 3 3 2
- a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 + a22 *a33 *b11
2 4 2 2 4 3
- a22 *a33 *b11 )/2 + u2*u3 *v2*( - 2*a11 *a33 *b11
3 3 3 4 2 2 3
+ 4*a11 *a22*a33 *b11 + 4*a11 *a33 *b11 - 2*a11 *a22 *a33 *b11
2 4 2 5 2 4
- 8*a11 *a22*a33 *b11 - 2*a11 *a33 *b11 + 4*a11*a22 *a33 *b11
5 2 5 3 4 3
+ 4*a11*a22*a33 *b11 - 2*a22 *a33 *b11) + u2*v2 *( - a11 *a33*b11
3 3 3 2 3 2 2 3
+ a11 *a22*a33*b11 + 3*a11 *a33 *b11 - 3*a11 *a22*a33 *b11
2 3 3 3 3 4 3 4 3
- 3*a11 *a33 *b11 + 3*a11*a22*a33 *b11 + a11*a33 *b11 - a22*a33 *b11 )
2 4 3 3 3 3 2 3
/2 + u2*v2*v3 *( - a11 *a33*b11 + 2*a11 *a22*a33*b11 + 2*a11 *a33 *b11
2 2 3 2 2 3 2 3 3
- a11 *a22 *a33*b11 - 4*a11 *a22*a33 *b11 - a11 *a33 *b11
2 2 3 3 3 2 3 3 4
+ 2*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 - a22 *a33 *b11 )/2 + u3 *(
4 4 3 4 3 5 2 2 4
a11 *a33 - 2*a11 *a22*a33 - 2*a11 *a33 + a11 *a22 *a33
2 5 2 6 2 5 6
+ 4*a11 *a22*a33 + a11 *a33 - 2*a11*a22 *a33 - 2*a11*a22*a33
2 6 3 4 3 3 3
+ a22 *a33 ) + u3 *v3*( - 2*a11 *a33 *b11 + 4*a11 *a22*a33 *b11
3 4 2 2 3 2 4
+ 4*a11 *a33 *b11 - 2*a11 *a22 *a33 *b11 - 8*a11 *a22*a33 *b11
2 5 2 4 5
- 2*a11 *a33 *b11 + 4*a11*a22 *a33 *b11 + 4*a11*a22*a33 *b11
2 5 2 2 4 2 2 3 2 2
- 2*a22 *a33 *b11) + u3 *v1 *( - a11 *a33 *b11 - a11 *a22*a33 *b11
3 3 2 2 2 2 2 2 3 2
+ 5*a11 *a33 *b11 + 2*a11 *a22 *a33 *b11 - a11 *a22*a33 *b11
2 4 2 2 3 2 4 2
- 7*a11 *a33 *b11 - 4*a11*a22 *a33 *b11 + 5*a11*a22*a33 *b11
5 2 2 4 2 5 2 2 2
+ 3*a11*a33 *b11 + 2*a22 *a33 *b11 - 3*a22*a33 *b11 )/2 + u3 *v2 *(
3 2 2 3 3 2 2 2 2 2
- a11 *a22*a33 *b11 + a11 *a33 *b11 + a11 *a22 *a33 *b11
2 3 2 2 4 2 2 3 2
+ a11 *a22*a33 *b11 - 2*a11 *a33 *b11 - 2*a11*a22 *a33 *b11
4 2 5 2 2 4 2 5 2
+ a11*a22*a33 *b11 + a11*a33 *b11 + a22 *a33 *b11 - a22*a33 *b11 ) +
2 2 3 2 2 3 3 2 2 2 2 2
u3 *v3 *( - 3*a11 *a22*a33 *b11 + 3*a11 *a33 *b11 + 3*a11 *a22 *a33 *b11
2 3 2 2 4 2 2 3 2
+ 3*a11 *a22*a33 *b11 - 6*a11 *a33 *b11 - 6*a11*a22 *a33 *b11
4 2 5 2 2 4 2
+ 3*a11*a22*a33 *b11 + 3*a11*a33 *b11 + 3*a22 *a33 *b11
5 2 2 3 3
- 3*a22*a33 *b11 )/2 + u3*v1 *v3*( - a11 *a22*a33*b11
3 2 3 2 2 3 2 2 3
+ a11 *a33 *b11 + a11 *a22 *a33*b11 + a11 *a22*a33 *b11
2 3 3 2 2 3 3 3
- 2*a11 *a33 *b11 - 2*a11*a22 *a33 *b11 + a11*a22*a33 *b11
4 3 2 3 3 4 3 2
+ a11*a33 *b11 + a22 *a33 *b11 - a22*a33 *b11 ) + u3*v2 *v3*(
4 3 3 3 3 2 3
- a11 *a33*b11 - a11 *a22*a33*b11 + 5*a11 *a33 *b11
2 2 3 2 2 3 2 3 3
+ 2*a11 *a22 *a33*b11 - a11 *a22*a33 *b11 - 7*a11 *a33 *b11
2 2 3 3 3 4 3
- 4*a11*a22 *a33 *b11 + 5*a11*a22*a33 *b11 + 3*a11*a33 *b11
2 3 3 4 3 3 4 3
+ 2*a22 *a33 *b11 - 3*a22*a33 *b11 )/2 + u3*v3 *( - a11 *a33*b11
3 2 3 2 2 3 2 2 3
+ 4*a11 *a33 *b11 + a11 *a22 *a33*b11 - 2*a11 *a22*a33 *b11
2 3 3 2 2 3 3 3
- 5*a11 *a33 *b11 - 2*a11*a22 *a33 *b11 + 4*a11*a22*a33 *b11
4 3 2 3 3 4 3 2 2
+ 2*a11*a33 *b11 + a22 *a33 *b11 - 2*a22*a33 *b11 )/2 + v1 *v2 *(
4 4 3 4 2 2 4 2 4
- a11 *b11 + 4*a11 *a33*b11 - 2*a11 *a22 *b11 + 4*a11 *a22*a33*b11
2 2 4 2 4 2 4
- 8*a11 *a33 *b11 + 4*a11*a22 *a33*b11 - 8*a11*a22*a33 *b11
3 4 2 2 4 3 4 4 4
+ 8*a11*a33 *b11 - 2*a22 *a33 *b11 + 4*a22*a33 *b11 - 3*a33 *b11 )/8 +
2 2 4 4 3 4 3 4
v1 *v3 *( - a11 *b11 - a11 *a22*b11 + 5*a11 *a33*b11
2 4 2 2 4 2 4
+ 3*a11 *a22*a33*b11 - 9*a11 *a33 *b11 - 3*a11*a22*a33 *b11
3 4 3 4 4 4 4
+ 7*a11*a33 *b11 + a22*a33 *b11 - 2*a33 *b11 )/8 + v2 *(
4 4 3 4 2 2 4 2 4
- a11 *b11 + 4*a11 *a33*b11 - 4*a11 *a22 *b11 + 8*a11 *a22*a33*b11
2 2 4 2 4 2 4
- 10*a11 *a33 *b11 + 8*a11*a22 *a33*b11 - 16*a11*a22*a33 *b11
3 4 2 2 4 3 4 4 4
+ 12*a11*a33 *b11 - 4*a22 *a33 *b11 + 8*a22*a33 *b11 - 5*a33 *b11 )/16
2 2 4 4 3 4 3 4
+ v2 *v3 *( - a11 *b11 - 2*a11 *a22*b11 + 6*a11 *a33*b11
2 2 4 2 4 2 2 4
- 2*a11 *a22 *b11 + 10*a11 *a22*a33*b11 - 14*a11 *a33 *b11
2 4 2 4 3 4
+ 4*a11*a22 *a33*b11 - 14*a11*a22*a33 *b11 + 14*a11*a33 *b11
2 2 4 3 4 4 4 4
- 2*a22 *a33 *b11 + 6*a22*a33 *b11 - 5*a33 *b11 )/8 + v3 *(
4 4 3 4 3 4 2 2 4
- a11 *b11 - 4*a11 *a22*b11 + 8*a11 *a33*b11 + a11 *a22 *b11
2 4 2 2 4 2 4
+ 10*a11 *a22*a33*b11 - 17*a11 *a33 *b11 - 2*a11*a22 *a33*b11
2 4 3 4 2 2 4
- 8*a11*a22*a33 *b11 + 14*a11*a33 *b11 + a22 *a33 *b11
3 4 4 4
+ 2*a22*a33 *b11 - 4*a33 *b11 )/16}
{HAM,FI} = 0
2 2 5 2 4 2 4 3 3 2 2
FI=u1 *v2 *( - 4*a11 *a33 + 16*a11 *a22*a33 + 4*a11 *a33 - 24*a11 *a22 *a33
3 3 2 3 2 2 2 3
- 16*a11 *a22*a33 + 16*a11 *a22 *a33 + 24*a11 *a22 *a33
4 2 3 3 4 3 2 2
- 4*a11*a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 ) + u1 *v3 *(
5 2 4 2 4 3 3 2 2
- 4*a11 *a33 + 16*a11 *a22*a33 + 4*a11 *a33 - 24*a11 *a22 *a33
3 3 2 3 2 2 2 3
- 16*a11 *a22*a33 + 16*a11 *a22 *a33 + 24*a11 *a22 *a33
4 2 3 3 4 3
- 4*a11*a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 ) + u1*u2*v1*v2*(
5 2 4 2 4 3 3 2 2
8*a11 *a33 - 32*a11 *a22*a33 - 8*a11 *a33 + 48*a11 *a22 *a33
3 3 2 3 2 2 2 3
+ 32*a11 *a22*a33 - 32*a11 *a22 *a33 - 48*a11 *a22 *a33
4 2 3 3 4 3
+ 8*a11*a22 *a33 + 32*a11*a22 *a33 - 8*a22 *a33 ) + u1*u3*v1*v3*(
5 2 4 2 4 3 3 2 2
8*a11 *a33 - 32*a11 *a22*a33 - 8*a11 *a33 + 48*a11 *a22 *a33
3 3 2 3 2 2 2 3
+ 32*a11 *a22*a33 - 32*a11 *a22 *a33 - 48*a11 *a22 *a33
4 2 3 3 4 3 3
+ 8*a11*a22 *a33 + 32*a11*a22 *a33 - 8*a22 *a33 ) + u1*v1 *(
5 4 4 2
- 4*a11 *a33*b11 + 16*a11 *a22*a33*b11 + 4*a11 *a33 *b11
3 2 3 2 2 3
- 24*a11 *a22 *a33*b11 - 16*a11 *a22*a33 *b11 + 16*a11 *a22 *a33*b11
2 2 2 4 3 2
+ 24*a11 *a22 *a33 *b11 - 4*a11*a22 *a33*b11 - 16*a11*a22 *a33 *b11
4 2 2 5 4
+ 4*a22 *a33 *b11) + u1*v1*v2 *( - 4*a11 *a33*b11 + 16*a11 *a22*a33*b11
4 2 3 2 3 2
+ 4*a11 *a33 *b11 - 24*a11 *a22 *a33*b11 - 16*a11 *a22*a33 *b11
2 3 2 2 2 4
+ 16*a11 *a22 *a33*b11 + 24*a11 *a22 *a33 *b11 - 4*a11*a22 *a33*b11
3 2 4 2 2 5
- 16*a11*a22 *a33 *b11 + 4*a22 *a33 *b11) + u1*v1*v3 *( - 4*a11 *a33*b11
4 4 2 3 2
+ 16*a11 *a22*a33*b11 + 4*a11 *a33 *b11 - 24*a11 *a22 *a33*b11
3 2 2 3 2 2 2
- 16*a11 *a22*a33 *b11 + 16*a11 *a22 *a33*b11 + 24*a11 *a22 *a33 *b11
4 3 2 4 2 2 2
- 4*a11*a22 *a33*b11 - 16*a11*a22 *a33 *b11 + 4*a22 *a33 *b11) + u2 *v1 *
5 2 4 2 4 3 3 2 2
( - 4*a11 *a33 + 16*a11 *a22*a33 + 4*a11 *a33 - 24*a11 *a22 *a33
3 3 2 3 2 2 2 3
- 16*a11 *a22*a33 + 16*a11 *a22 *a33 + 24*a11 *a22 *a33
4 2 3 3 4 3 2 2
- 4*a11*a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 ) + u2 *v3 *(
5 2 4 2 4 3 3 2 2
- 4*a11 *a33 + 16*a11 *a22*a33 + 4*a11 *a33 - 24*a11 *a22 *a33
3 3 2 3 2 2 2 3
- 16*a11 *a22*a33 + 16*a11 *a22 *a33 + 24*a11 *a22 *a33
4 2 3 3 4 3
- 4*a11*a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 ) + u2*u3*v2*v3*(
5 2 4 2 4 3 3 2 2
8*a11 *a33 - 32*a11 *a22*a33 - 8*a11 *a33 + 48*a11 *a22 *a33
3 3 2 3 2 2 2 3
+ 32*a11 *a22*a33 - 32*a11 *a22 *a33 - 48*a11 *a22 *a33
4 2 3 3 4 3 2 2
+ 8*a11*a22 *a33 + 32*a11*a22 *a33 - 8*a22 *a33 ) + u3 *v1 *(
5 2 4 2 4 3 3 2 2
- 4*a11 *a33 + 16*a11 *a22*a33 + 4*a11 *a33 - 24*a11 *a22 *a33
3 3 2 3 2 2 2 3
- 16*a11 *a22*a33 + 16*a11 *a22 *a33 + 24*a11 *a22 *a33
4 2 3 3 4 3 2 2
- 4*a11*a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 ) + u3 *v2 *(
5 2 4 2 4 3 3 2 2
- 4*a11 *a33 + 16*a11 *a22*a33 + 4*a11 *a33 - 24*a11 *a22 *a33
3 3 2 3 2 2 2 3
- 16*a11 *a22*a33 + 16*a11 *a22 *a33 + 24*a11 *a22 *a33
4 2 3 3 4 3 2
- 4*a11*a22 *a33 - 16*a11*a22 *a33 + 4*a22 *a33 ) + u3*v1 *v3*(
4 4 2 3 2
4*a11 *a22*a33*b11 - 4*a11 *a33 *b11 - 12*a11 *a22 *a33*b11
3 2 3 3 2 3
+ 8*a11 *a22*a33 *b11 + 4*a11 *a33 *b11 + 12*a11 *a22 *a33*b11
2 3 4 3 2
- 12*a11 *a22*a33 *b11 - 4*a11*a22 *a33*b11 - 8*a11*a22 *a33 *b11
2 3 4 2 3 3 2
+ 12*a11*a22 *a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11) + u3*v2 *v3*(
4 4 2 3 2
4*a11 *a22*a33*b11 - 4*a11 *a33 *b11 - 12*a11 *a22 *a33*b11
3 2 3 3 2 3
+ 8*a11 *a22*a33 *b11 + 4*a11 *a33 *b11 + 12*a11 *a22 *a33*b11
2 3 4 3 2
- 12*a11 *a22*a33 *b11 - 4*a11*a22 *a33*b11 - 8*a11*a22 *a33 *b11
2 3 4 2 3 3 3
+ 12*a11*a22 *a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11) + u3*v3 *(
4 4 2 3 2
4*a11 *a22*a33*b11 - 4*a11 *a33 *b11 - 12*a11 *a22 *a33*b11
3 2 3 3 2 3
+ 8*a11 *a22*a33 *b11 + 4*a11 *a33 *b11 + 12*a11 *a22 *a33*b11
2 3 4 3 2
- 12*a11 *a22*a33 *b11 - 4*a11*a22 *a33*b11 - 8*a11*a22 *a33 *b11
2 3 4 2 3 3 2 2
+ 12*a11*a22 *a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11) + v1 *v2 *(
5 2 4 2 4 2 3 2
- a11 *b11 + 2*a11 *a22*b11 + 3*a11 *a33*b11 - 8*a11 *a22*a33*b11
3 2 2 2 3 2 2 2 2
- 2*a11 *a33 *b11 - 2*a11 *a22 *b11 + 6*a11 *a22 *a33*b11
2 2 2 4 2 2 2 2
+ 6*a11 *a22*a33 *b11 + a11*a22 *b11 - 6*a11*a22 *a33 *b11
4 2 3 2 2 2 2 5 2
- a22 *a33*b11 + 2*a22 *a33 *b11 ) + v1 *v3 *( - a11 *b11
4 2 4 2 3 2 2
+ 4*a11 *a22*b11 + a11 *a33*b11 - 5*a11 *a22 *b11
3 2 3 2 2 2 3 2
- 6*a11 *a22*a33*b11 + a11 *a33 *b11 + 2*a11 *a22 *b11
2 2 2 2 3 2 3 2
+ 9*a11 *a22 *a33*b11 - a11 *a33 *b11 - 4*a11*a22 *a33*b11
2 2 2 3 2 3 2 2
- 3*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 + 2*a22 *a33 *b11
2 3 2 4 5 2 4 2 4 2
- a22 *a33 *b11 ) + v2 *( - a11 *b11 + 2*a11 *a22*b11 + 3*a11 *a33*b11
3 2 3 2 2 2 3 2
- 8*a11 *a22*a33*b11 - 2*a11 *a33 *b11 - 2*a11 *a22 *b11
2 2 2 2 2 2 4 2
+ 6*a11 *a22 *a33*b11 + 6*a11 *a22*a33 *b11 + a11*a22 *b11
2 2 2 4 2 3 2 2 2 2
- 6*a11*a22 *a33 *b11 - a22 *a33*b11 + 2*a22 *a33 *b11 ) + v2 *v3 *(
5 2 4 2 4 2 3 2 2
- 2*a11 *b11 + 6*a11 *a22*b11 + 4*a11 *a33*b11 - 5*a11 *a22 *b11
3 2 3 2 2 2 2 2
- 14*a11 *a22*a33*b11 - a11 *a33 *b11 + 15*a11 *a22 *a33*b11
2 2 2 2 3 2 4 2
+ 6*a11 *a22*a33 *b11 - a11 *a33 *b11 + a11*a22 *b11
3 2 2 2 2 3 2
- 4*a11*a22 *a33*b11 - 9*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11
4 2 3 2 2 2 3 2 4 5 2
- a22 *a33*b11 + 4*a22 *a33 *b11 - a22 *a33 *b11 ) + v3 *( - a11 *b11
4 2 4 2 3 2 2
+ 4*a11 *a22*b11 + a11 *a33*b11 - 5*a11 *a22 *b11
3 2 3 2 2 2 3 2
- 6*a11 *a22*a33*b11 + a11 *a33 *b11 + 2*a11 *a22 *b11
2 2 2 2 3 2 3 2
+ 9*a11 *a22 *a33*b11 - a11 *a33 *b11 - 4*a11*a22 *a33*b11
2 2 2 3 2 3 2 2
- 3*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 + 2*a22 *a33 *b11
2 3 2
- a22 *a33 *b11 )
= a product of the elements of: {4,
- a11 + a22,
a11 - a22,
a11 - a33,
2 2 2 2 2 2 2
u1 *v2 *(a11 *a33 - 2*a11*a22*a33 + a22 *a33 )
2 2 2 2 2 2 2
+ u1 *v3 *(a11 *a33 - 2*a11*a22*a33 + a22 *a33 )
2 2 2 2 2
+ u1*u2*v1*v2*( - 2*a11 *a33 + 4*a11*a22*a33 - 2*a22 *a33 )
2 2 2 2 2
+ u1*u3*v1*v3*( - 2*a11 *a33 + 4*a11*a22*a33 - 2*a22 *a33 )
3 2 2
+ u1*v1 *(a11 *a33*b11 - 2*a11*a22*a33*b11 + a22 *a33*b11)
2 2 2
+ u1*v1*v2 *(a11 *a33*b11 - 2*a11*a22*a33*b11 + a22 *a33*b11)
2 2 2
+ u1*v1*v3 *(a11 *a33*b11 - 2*a11*a22*a33*b11 + a22 *a33*b11)
2 2 2 2 2 2 2
+ u2 *v1 *(a11 *a33 - 2*a11*a22*a33 + a22 *a33 )
2 2 2 2 2 2 2
+ u2 *v3 *(a11 *a33 - 2*a11*a22*a33 + a22 *a33 )
2 2 2 2 2
+ u2*u3*v2*v3*( - 2*a11 *a33 + 4*a11*a22*a33 - 2*a22 *a33 )
2 2 2 2 2 2 2
+ u3 *v1 *(a11 *a33 - 2*a11*a22*a33 + a22 *a33 )
2 2 2 2 2 2 2 2
+ u3 *v2 *(a11 *a33 - 2*a11*a22*a33 + a22 *a33 ) + u3*v1 *v3
2 2 2 2
*( - a11*a22*a33*b11 + a11*a33 *b11 + a22 *a33*b11 - a22*a33 *b11) + u3*v2
2 2 2
*v3*( - a11*a22*a33*b11 + a11*a33 *b11 + a22 *a33*b11 - a22*a33 *b11)
3 2 2 2
+ u3*v3 *( - a11*a22*a33*b11 + a11*a33 *b11 + a22 *a33*b11 - a22*a33 *b11)
2 2 2 2 2 2
2 2 a11 *b11 - 2*a11*a33*b11 - a22 *b11 + 2*a22*a33*b11
+ v1 *v2 *---------------------------------------------------------
4
2 2 2 2 2 2
2 2 a11 *b11 - 2*a11*a22*b11 + 2*a22*a33*b11 - a33 *b11
+ v1 *v3 *---------------------------------------------------------
4
2 2 2 2 2 2
4 a11 *b11 - 2*a11*a33*b11 - a22 *b11 + 2*a22*a33*b11 2 2
+ v2 *--------------------------------------------------------- + v2 *v3 *(
4
2 2 2 2 2 2 2
2*a11 *b11 - 2*a11*a22*b11 - 2*a11*a33*b11 - a22 *b11 + 4*a22*a33*b11
2 2
- a33 *b11 )/4
2 2 2 2 2 2
4 a11 *b11 - 2*a11*a22*b11 + 2*a22*a33*b11 - a33 *b11
+ v3 *---------------------------------------------------------}
4
{HAM,FI} = 0
2 2 5 2 4 2 4 3
FI=u1 *v2 *( - 4*a11 *a33 *b11 + 12*a11 *a22*a33 *b11 + 8*a11 *a33 *b11
3 2 2 3 3 3 4
- 12*a11 *a22 *a33 *b11 - 24*a11 *a22*a33 *b11 - 4*a11 *a33 *b11
2 3 2 2 2 3
+ 4*a11 *a22 *a33 *b11 + 24*a11 *a22 *a33 *b11
2 4 3 3 2 4
+ 12*a11 *a22*a33 *b11 - 8*a11*a22 *a33 *b11 - 12*a11*a22 *a33 *b11
3 4 2 2 5 2
+ 4*a22 *a33 *b11) + u1 *v3 *( - 4*a11 *a33 *b11
4 2 4 3 3 2 2
+ 12*a11 *a22*a33 *b11 + 8*a11 *a33 *b11 - 12*a11 *a22 *a33 *b11
3 3 3 4 2 3 2
- 24*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 + 4*a11 *a22 *a33 *b11
2 2 3 2 4 3 3
+ 24*a11 *a22 *a33 *b11 + 12*a11 *a22*a33 *b11 - 8*a11*a22 *a33 *b11
2 4 3 4 2 5 3
- 12*a11*a22 *a33 *b11 + 4*a22 *a33 *b11) + u1*u2 *v1*(4*a11 *a33
4 3 4 4 3 2 3 3 4
- 16*a11 *a22*a33 - 4*a11 *a33 + 24*a11 *a22 *a33 + 16*a11 *a22*a33
2 3 3 2 2 4 4 3
- 16*a11 *a22 *a33 - 24*a11 *a22 *a33 + 4*a11*a22 *a33
3 4 4 4 5 2
+ 16*a11*a22 *a33 - 4*a22 *a33 ) + u1*u2*v1*v2*(4*a11 *a33 *b11
4 2 4 3 3 2 2
- 12*a11 *a22*a33 *b11 - 8*a11 *a33 *b11 + 12*a11 *a22 *a33 *b11
3 3 3 4 2 3 2
+ 24*a11 *a22*a33 *b11 + 4*a11 *a33 *b11 - 4*a11 *a22 *a33 *b11
2 2 3 2 4 3 3
- 24*a11 *a22 *a33 *b11 - 12*a11 *a22*a33 *b11 + 8*a11*a22 *a33 *b11
2 4 3 4 2 5 3
+ 12*a11*a22 *a33 *b11 - 4*a22 *a33 *b11) + u1*u3 *v1*(4*a11 *a33
4 3 4 4 3 2 3 3 4
- 12*a11 *a22*a33 - 8*a11 *a33 + 12*a11 *a22 *a33 + 24*a11 *a22*a33
3 5 2 3 3 2 2 4 2 5
+ 4*a11 *a33 - 4*a11 *a22 *a33 - 24*a11 *a22 *a33 - 12*a11 *a22*a33
3 4 2 5 3 5
+ 8*a11*a22 *a33 + 12*a11*a22 *a33 - 4*a22 *a33 ) + u1*u3*v1*v3*(
5 2 4 2 4 3
4*a11 *a33 *b11 - 12*a11 *a22*a33 *b11 - 8*a11 *a33 *b11
3 2 2 3 3 3 4
+ 12*a11 *a22 *a33 *b11 + 24*a11 *a22*a33 *b11 + 4*a11 *a33 *b11
2 3 2 2 2 3 2 4
- 4*a11 *a22 *a33 *b11 - 24*a11 *a22 *a33 *b11 - 12*a11 *a22*a33 *b11
3 3 2 4 3 4 3
+ 8*a11*a22 *a33 *b11 + 12*a11*a22 *a33 *b11 - 4*a22 *a33 *b11) + u1*v1 *
4 2 4 2 2 3 2 2
(4*a11 *a22*a33*b11 - 4*a11 *a33 *b11 - 8*a11 *a22 *a33*b11
3 3 2 2 3 2 2 2 2 2
+ 8*a11 *a33 *b11 + 4*a11 *a22 *a33*b11 + 12*a11 *a22 *a33 *b11
2 3 2 2 4 2 3 2 2
- 12*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11*a22 *a33 *b11
4 2 3 3 2 2 4 2 2
+ 8*a11*a22*a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11 ) + u1*v1*v2 *(
5 2 4 2 4 2 2
a11 *a33*b11 + 2*a11 *a22*a33*b11 - 7*a11 *a33 *b11
3 2 2 3 2 2 3 3 2
- 7*a11 *a22 *a33*b11 + 6*a11 *a22*a33 *b11 + 11*a11 *a33 *b11
2 3 2 2 2 2 2 2 3 2
+ 4*a11 *a22 *a33*b11 + 9*a11 *a22 *a33 *b11 - 18*a11 *a22*a33 *b11
2 4 2 3 2 2 2 3 2
- 5*a11 *a33 *b11 - 8*a11*a22 *a33 *b11 + 3*a11*a22 *a33 *b11
4 2 3 3 2 2 4 2 2
+ 10*a11*a22*a33 *b11 + 4*a22 *a33 *b11 - 5*a22 *a33 *b11 ) + u1*v1*v3
5 2 4 2 4 2 2 3 2 2
*(a11 *a33*b11 + a11 *a22*a33*b11 - 6*a11 *a33 *b11 - 5*a11 *a22 *a33*b11
3 2 2 3 3 2 2 3 2
+ 6*a11 *a22*a33 *b11 + 9*a11 *a33 *b11 + 3*a11 *a22 *a33*b11
2 2 2 2 2 3 2 2 4 2
+ 6*a11 *a22 *a33 *b11 - 15*a11 *a22*a33 *b11 - 4*a11 *a33 *b11
3 2 2 2 3 2 4 2
- 6*a11*a22 *a33 *b11 + 3*a11*a22 *a33 *b11 + 8*a11*a22*a33 *b11
3 3 2 2 4 2 3 5 3
+ 3*a22 *a33 *b11 - 4*a22 *a33 *b11 ) + u2 *v2*(4*a11 *a33
4 3 4 4 3 2 3 3 4
- 16*a11 *a22*a33 - 4*a11 *a33 + 24*a11 *a22 *a33 + 16*a11 *a22*a33
2 3 3 2 2 4 4 3
- 16*a11 *a22 *a33 - 24*a11 *a22 *a33 + 4*a11*a22 *a33
3 4 4 4 2 5 3
+ 16*a11*a22 *a33 - 4*a22 *a33 ) + u2 *u3*v3*(4*a11 *a33
4 3 4 4 3 2 3 3 4
- 16*a11 *a22*a33 - 4*a11 *a33 + 24*a11 *a22 *a33 + 16*a11 *a22*a33
2 3 3 2 2 4 4 3
- 16*a11 *a22 *a33 - 24*a11 *a22 *a33 + 4*a11*a22 *a33
3 4 4 4 2 5 3
+ 16*a11*a22 *a33 - 4*a22 *a33 ) + u2*u3 *v2*(4*a11 *a33
4 3 4 4 3 2 3 3 4
- 12*a11 *a22*a33 - 8*a11 *a33 + 12*a11 *a22 *a33 + 24*a11 *a22*a33
3 5 2 3 3 2 2 4 2 5
+ 4*a11 *a33 - 4*a11 *a22 *a33 - 24*a11 *a22 *a33 - 12*a11 *a22*a33
3 4 2 5 3 5 3
+ 8*a11*a22 *a33 + 12*a11*a22 *a33 - 4*a22 *a33 ) + u2*v2 *(
5 2 4 2 4 2 2
a11 *a33*b11 - 2*a11 *a22*a33*b11 - 3*a11 *a33 *b11
3 2 2 3 2 2 3 3 2
+ a11 *a22 *a33*b11 + 6*a11 *a22*a33 *b11 + 3*a11 *a33 *b11
2 2 2 2 2 3 2 2 4 2
- 3*a11 *a22 *a33 *b11 - 6*a11 *a22*a33 *b11 - a11 *a33 *b11
2 3 2 4 2 2 4 2
+ 3*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 - a22 *a33 *b11 ) + u2*v2
2 5 2 4 2 4 2 2
*v3 *(a11 *a33*b11 - 3*a11 *a22*a33*b11 - 2*a11 *a33 *b11
3 2 2 3 2 2 3 3 2
+ 3*a11 *a22 *a33*b11 + 6*a11 *a22*a33 *b11 + a11 *a33 *b11
2 3 2 2 2 2 2 2 3 2
- a11 *a22 *a33*b11 - 6*a11 *a22 *a33 *b11 - 3*a11 *a22*a33 *b11
3 2 2 2 3 2 3 3 2 3
+ 2*a11*a22 *a33 *b11 + 3*a11*a22 *a33 *b11 - a22 *a33 *b11 ) + u3
5 3 4 3 4 4 3 2 3
*v3*(4*a11 *a33 - 12*a11 *a22*a33 - 8*a11 *a33 + 12*a11 *a22 *a33
3 4 3 5 2 3 3 2 2 4
+ 24*a11 *a22*a33 + 4*a11 *a33 - 4*a11 *a22 *a33 - 24*a11 *a22 *a33
2 5 3 4 2 5 3 5
- 12*a11 *a22*a33 + 8*a11*a22 *a33 + 12*a11*a22 *a33 - 4*a22 *a33 )
2 2 4 2 4 3 3 2 2
+ u3 *v1 *(4*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11 *a22 *a33 *b11
3 4 2 3 2 2 2 3
+ 8*a11 *a33 *b11 + 4*a11 *a22 *a33 *b11 + 12*a11 *a22 *a33 *b11
2 4 2 5 3 3
- 12*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11*a22 *a33 *b11
5 3 4 2 5 2
+ 8*a11*a22*a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11) + u3
2 4 2 4 3 3 2 2
*v2 *(4*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11 *a22 *a33 *b11
3 4 2 3 2 2 2 3
+ 8*a11 *a33 *b11 + 4*a11 *a22 *a33 *b11 + 12*a11 *a22 *a33 *b11
2 4 2 5 3 3
- 12*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11*a22 *a33 *b11
5 3 4 2 5 2 2
+ 8*a11*a22*a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11) + u3 *v3 *(
4 2 4 3 3 2 2
4*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11 *a22 *a33 *b11
3 4 2 3 2 2 2 3
+ 8*a11 *a33 *b11 + 4*a11 *a22 *a33 *b11 + 12*a11 *a22 *a33 *b11
2 4 2 5 3 3
- 12*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11*a22 *a33 *b11
5 3 4 2 5 2
+ 8*a11*a22*a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11) + u3*v1 *v3*(
4 2 4 2 2 3 2 2
4*a11 *a22*a33*b11 - 4*a11 *a33 *b11 - 8*a11 *a22 *a33*b11
3 3 2 2 3 2 2 2 2 2
+ 8*a11 *a33 *b11 + 4*a11 *a22 *a33*b11 + 12*a11 *a22 *a33 *b11
2 3 2 2 4 2 3 2 2
- 12*a11 *a22*a33 *b11 - 4*a11 *a33 *b11 - 8*a11*a22 *a33 *b11
4 2 3 3 2 2 4 2 2
+ 8*a11*a22*a33 *b11 + 4*a22 *a33 *b11 - 4*a22 *a33 *b11 ) + u3*v2 *v3*
5 2 4 2 4 2 2
(a11 *a33*b11 + 2*a11 *a22*a33*b11 - 7*a11 *a33 *b11
3 2 2 3 2 2 3 3 2
- 7*a11 *a22 *a33*b11 + 6*a11 *a22*a33 *b11 + 11*a11 *a33 *b11
2 3 2 2 2 2 2 2 3 2
+ 4*a11 *a22 *a33*b11 + 9*a11 *a22 *a33 *b11 - 18*a11 *a22*a33 *b11
2 4 2 3 2 2 2 3 2
- 5*a11 *a33 *b11 - 8*a11*a22 *a33 *b11 + 3*a11*a22 *a33 *b11
4 2 3 3 2 2 4 2 3
+ 10*a11*a22*a33 *b11 + 4*a22 *a33 *b11 - 5*a22 *a33 *b11 ) + u3*v3 *(
5 2 4 2 4 2 2
a11 *a33*b11 + a11 *a22*a33*b11 - 6*a11 *a33 *b11
3 2 2 3 2 2 3 3 2
- 5*a11 *a22 *a33*b11 + 6*a11 *a22*a33 *b11 + 9*a11 *a33 *b11
2 3 2 2 2 2 2 2 3 2
+ 3*a11 *a22 *a33*b11 + 6*a11 *a22 *a33 *b11 - 15*a11 *a22*a33 *b11
2 4 2 3 2 2 2 3 2
- 4*a11 *a33 *b11 - 6*a11*a22 *a33 *b11 + 3*a11*a22 *a33 *b11
4 2 3 3 2 2 4 2 2 2
+ 8*a11*a22*a33 *b11 + 3*a22 *a33 *b11 - 4*a22 *a33 *b11 ) + v1 *v2 *(
3 2 3 3 3 3 2 3 2 3 3
a11 *a22 *b11 - 2*a11 *a22*a33*b11 + a11 *a33 *b11 - a11 *a22 *b11
2 2 3 2 3 3 3 3
+ 3*a11 *a22*a33 *b11 - 2*a11 *a33 *b11 + 2*a11*a22 *a33*b11
2 2 3 4 3 3 2 3
- 3*a11*a22 *a33 *b11 + a11*a33 *b11 - a22 *a33 *b11
2 3 3 4 3 2 2 4 3
+ 2*a22 *a33 *b11 - a22*a33 *b11 ) + v1 *v3 *(a11 *a22*b11
4 3 3 2 3 3 3 3 2 3
- a11 *a33*b11 - a11 *a22 *b11 - 2*a11 *a22*a33*b11 + 3*a11 *a33 *b11
2 2 3 2 3 3 2 2 3
+ 3*a11 *a22 *a33*b11 - 3*a11 *a33 *b11 - 3*a11*a22 *a33 *b11
3 3 4 3 2 3 3 4 3
+ 2*a11*a22*a33 *b11 + a11*a33 *b11 + a22 *a33 *b11 - a22*a33 *b11 ) +
4 3 2 3 3 3 3 2 3 2 3 3
v2 *(a11 *a22 *b11 - 2*a11 *a22*a33*b11 + a11 *a33 *b11 - a11 *a22 *b11
2 2 3 2 3 3 3 3
+ 3*a11 *a22*a33 *b11 - 2*a11 *a33 *b11 + 2*a11*a22 *a33*b11
2 2 3 4 3 3 2 3
- 3*a11*a22 *a33 *b11 + a11*a33 *b11 - a22 *a33 *b11
2 3 3 4 3 2 2 4 3
+ 2*a22 *a33 *b11 - a22*a33 *b11 ) + v2 *v3 *(a11 *a22*b11
4 3 3 3 3 2 3 2 3 3
- a11 *a33*b11 - 4*a11 *a22*a33*b11 + 4*a11 *a33 *b11 - a11 *a22 *b11
2 2 3 2 2 3 2 3 3
+ 3*a11 *a22 *a33*b11 + 3*a11 *a22*a33 *b11 - 5*a11 *a33 *b11
3 3 2 2 3 3 3
+ 2*a11*a22 *a33*b11 - 6*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11
4 3 3 2 3 2 3 3 4 3
+ 2*a11*a33 *b11 - a22 *a33 *b11 + 3*a22 *a33 *b11 - 2*a22*a33 *b11 )
4 4 3 4 3 3 2 3 3 3
+ v3 *(a11 *a22*b11 - a11 *a33*b11 - a11 *a22 *b11 - 2*a11 *a22*a33*b11
3 2 3 2 2 3 2 3 3
+ 3*a11 *a33 *b11 + 3*a11 *a22 *a33*b11 - 3*a11 *a33 *b11
2 2 3 3 3 4 3
- 3*a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 + a11*a33 *b11
2 3 3 4 3
+ a22 *a33 *b11 - a22*a33 *b11 )
= a product of the elements of: {4,
a11 - a33,
a11 - a22,
2 2 3 2 2 2 2 3
u1 *v2 *( - a11 *a33 *b11 + 2*a11 *a22*a33 *b11 + a11 *a33 *b11
2 2 3 2 3 2 2
- a11*a22 *a33 *b11 - 2*a11*a22*a33 *b11 + a22 *a33 *b11) + u1 *v3
3 2 2 2 2 3 2 2
*( - a11 *a33 *b11 + 2*a11 *a22*a33 *b11 + a11 *a33 *b11 - a11*a22 *a33 *b11
3 2 3
- 2*a11*a22*a33 *b11 + a22 *a33 *b11)
2 3 3 2 3 2 3 3 3
+ u1*u2 *v1*(a11 *a33 - 3*a11 *a22*a33 + 3*a11*a22 *a33 - a22 *a33 ) + u1
3 2 2 2 2 3
*u2*v1*v2*(a11 *a33 *b11 - 2*a11 *a22*a33 *b11 - a11 *a33 *b11
2 2 3 2 3
+ a11*a22 *a33 *b11 + 2*a11*a22*a33 *b11 - a22 *a33 *b11) + u1
2 3 3 2 3 2 4 2 3
*u3 *v1*(a11 *a33 - 2*a11 *a22*a33 - a11 *a33 + a11*a22 *a33
4 2 4 3 2
+ 2*a11*a22*a33 - a22 *a33 ) + u1*u3*v1*v3*(a11 *a33 *b11
2 2 2 3 2 2
- 2*a11 *a22*a33 *b11 - a11 *a33 *b11 + a11*a22 *a33 *b11
3 2 3 3 2 2
+ 2*a11*a22*a33 *b11 - a22 *a33 *b11) + u1*v1 *(a11 *a22*a33*b11
2 2 2 2 2 3 2 2 2 2
- a11 *a33 *b11 - a11*a22 *a33*b11 + a11*a33 *b11 + a22 *a33 *b11
3 2 2 3 2 2 2
- a22*a33 *b11 ) + u1*v1*v2 *(a11 *a33*b11 + 3*a11 *a22*a33*b11
2 2 2 2 2 2 2
- 6*a11 *a33 *b11 - 4*a11*a22 *a33*b11 + 2*a11*a22*a33 *b11
3 2 2 2 2 3 2 2
+ 5*a11*a33 *b11 + 4*a22 *a33 *b11 - 5*a22*a33 *b11 )/4 + u1*v1*v3 *(
3 2 2 2 2 2 2
a11 *a33*b11 + 2*a11 *a22*a33*b11 - 5*a11 *a33 *b11
2 2 2 2 3 2
- 3*a11*a22 *a33*b11 + 2*a11*a22*a33 *b11 + 4*a11*a33 *b11
2 2 2 3 2
+ 3*a22 *a33 *b11 - 4*a22*a33 *b11 )/4
3 3 3 2 3 2 3 3 3
+ u2 *v2*(a11 *a33 - 3*a11 *a22*a33 + 3*a11*a22 *a33 - a22 *a33 )
2 3 3 2 3 2 3 3 3
+ u2 *u3*v3*(a11 *a33 - 3*a11 *a22*a33 + 3*a11*a22 *a33 - a22 *a33 ) + u2
2 3 3 2 3 2 4 2 3
*u3 *v2*(a11 *a33 - 2*a11 *a22*a33 - a11 *a33 + a11*a22 *a33
4 2 4 3 3 2
+ 2*a11*a22*a33 - a22 *a33 ) + u2*v2 *(a11 *a33*b11
2 2 2 2 2 2 2
- a11 *a22*a33*b11 - 2*a11 *a33 *b11 + 2*a11*a22*a33 *b11
3 2 3 2 2 3 2
+ a11*a33 *b11 - a22*a33 *b11 )/4 + u2*v2*v3 *(a11 *a33*b11
2 2 2 2 2 2 2
- 2*a11 *a22*a33*b11 - a11 *a33 *b11 + a11*a22 *a33*b11
2 2 2 2 2 3 3 3
+ 2*a11*a22*a33 *b11 - a22 *a33 *b11 )/4 + u3 *v3*(a11 *a33
2 3 2 4 2 3 4
- 2*a11 *a22*a33 - a11 *a33 + a11*a22 *a33 + 2*a11*a22*a33
2 4 2 2 2 2 2 3
- a22 *a33 ) + u3 *v1 *(a11 *a22*a33 *b11 - a11 *a33 *b11
2 2 4 2 3 4 2
- a11*a22 *a33 *b11 + a11*a33 *b11 + a22 *a33 *b11 - a22*a33 *b11) + u3
2 2 2 2 3 2 2 4
*v2 *(a11 *a22*a33 *b11 - a11 *a33 *b11 - a11*a22 *a33 *b11 + a11*a33 *b11
2 3 4 2 2 2 2
+ a22 *a33 *b11 - a22*a33 *b11) + u3 *v3 *(a11 *a22*a33 *b11
2 3 2 2 4 2 3
- a11 *a33 *b11 - a11*a22 *a33 *b11 + a11*a33 *b11 + a22 *a33 *b11
4 2 2 2 2 2 2
- a22*a33 *b11) + u3*v1 *v3*(a11 *a22*a33*b11 - a11 *a33 *b11
2 2 3 2 2 2 2 3 2
- a11*a22 *a33*b11 + a11*a33 *b11 + a22 *a33 *b11 - a22*a33 *b11 ) +
2 3 2 2 2 2 2 2
u3*v2 *v3*(a11 *a33*b11 + 3*a11 *a22*a33*b11 - 6*a11 *a33 *b11
2 2 2 2 3 2
- 4*a11*a22 *a33*b11 + 2*a11*a22*a33 *b11 + 5*a11*a33 *b11
2 2 2 3 2 3 3 2
+ 4*a22 *a33 *b11 - 5*a22*a33 *b11 )/4 + u3*v3 *(a11 *a33*b11
2 2 2 2 2 2 2
+ 2*a11 *a22*a33*b11 - 5*a11 *a33 *b11 - 3*a11*a22 *a33*b11
2 2 3 2 2 2 2
+ 2*a11*a22*a33 *b11 + 4*a11*a33 *b11 + 3*a22 *a33 *b11
3 2 2 2 2 3 3
- 4*a22*a33 *b11 )/4 + v1 *v2 *(a11*a22 *b11 - 2*a11*a22*a33*b11
2 3 2 3 2 3 3 3 2
+ a11*a33 *b11 - a22 *a33*b11 + 2*a22*a33 *b11 - a33 *b11 )/4 + v1
2 2 3 2 3 3 2 3
*v3 *(a11 *a22*b11 - a11 *a33*b11 - 2*a11*a22*a33*b11 + 2*a11*a33 *b11
2 3 3 3 4 2 3
+ a22*a33 *b11 - a33 *b11 )/4 + v2 *(a11*a22 *b11
3 2 3 2 3 2 3
- 2*a11*a22*a33*b11 + a11*a33 *b11 - a22 *a33*b11 + 2*a22*a33 *b11
3 3 2 2 2 3 2 3 2 3
- a33 *b11 )/4 + v2 *v3 *(a11 *a22*b11 - a11 *a33*b11 + a11*a22 *b11
3 2 3 2 3 2 3
- 4*a11*a22*a33*b11 + 3*a11*a33 *b11 - a22 *a33*b11 + 3*a22*a33 *b11
3 3 4 2 3 2 3
- 2*a33 *b11 )/4 + v3 *(a11 *a22*b11 - a11 *a33*b11
3 2 3 2 3 3 3
- 2*a11*a22*a33*b11 + 2*a11*a33 *b11 + a22*a33 *b11 - a33 *b11 )/4}
{HAM,FI} = 0
3 5 3 4 3 4 4
FI=u1 *v1*(32*a11 *a33 *b11 - 96*a11 *a22*a33 *b11 - 64*a11 *a33 *b11
3 2 3 3 4 3 5
+ 96*a11 *a22 *a33 *b11 + 192*a11 *a22*a33 *b11 + 32*a11 *a33 *b11
2 3 3 2 2 4
- 32*a11 *a22 *a33 *b11 - 192*a11 *a22 *a33 *b11
2 5 3 4 2 5
- 96*a11 *a22*a33 *b11 + 64*a11*a22 *a33 *b11 + 96*a11*a22 *a33 *b11
3 5 2 2 5 4 4 4
- 32*a22 *a33 *b11) + u1 *u2 *(32*a11 *a33 - 128*a11 *a22*a33
4 5 3 2 4 3 5
- 32*a11 *a33 + 192*a11 *a22 *a33 + 128*a11 *a22*a33
2 3 4 2 2 5 4 4
- 128*a11 *a22 *a33 - 192*a11 *a22 *a33 + 32*a11*a22 *a33
3 5 4 5 2 2 5 4
+ 128*a11*a22 *a33 - 32*a22 *a33 ) + u1 *u3 *(32*a11 *a33
4 4 4 5 3 2 4 3 5
- 96*a11 *a22*a33 - 64*a11 *a33 + 96*a11 *a22 *a33 + 192*a11 *a22*a33
3 6 2 3 4 2 2 5
+ 32*a11 *a33 - 32*a11 *a22 *a33 - 192*a11 *a22 *a33
2 6 3 5 2 6 3 6
- 96*a11 *a22*a33 + 64*a11*a22 *a33 + 96*a11*a22 *a33 - 32*a22 *a33 )
2 2 5 2 2 4 3 2 3 2 2 2
+ u1 *v2 *(8*a11 *a33 *b11 - 40*a11 *a33 *b11 - 24*a11 *a22 *a33 *b11
3 3 2 3 4 2
+ 48*a11 *a22*a33 *b11 + 56*a11 *a33 *b11
2 3 2 2 2 2 3 2
+ 16*a11 *a22 *a33 *b11 + 24*a11 *a22 *a33 *b11
2 4 2 2 5 2
- 96*a11 *a22*a33 *b11 - 24*a11 *a33 *b11
3 3 2 2 4 2
- 32*a11*a22 *a33 *b11 + 24*a11*a22 *a33 *b11
5 2 3 4 2 2 5 2
+ 48*a11*a22*a33 *b11 + 16*a22 *a33 *b11 - 24*a22 *a33 *b11 ) +
2 2 5 2 2 4 2 2 4 3 2
u1 *v3 *(8*a11 *a33 *b11 - 8*a11 *a22*a33 *b11 - 32*a11 *a33 *b11
3 2 2 2 3 3 2
- 8*a11 *a22 *a33 *b11 + 48*a11 *a22*a33 *b11
3 4 2 2 3 2 2
+ 40*a11 *a33 *b11 + 8*a11 *a22 *a33 *b11
2 4 2 2 5 2
- 72*a11 *a22*a33 *b11 - 16*a11 *a33 *b11
3 3 2 2 4 2
- 16*a11*a22 *a33 *b11 + 24*a11*a22 *a33 *b11
5 2 3 4 2 2 5 2
+ 32*a11*a22*a33 *b11 + 8*a22 *a33 *b11 - 16*a22 *a33 *b11 ) +
2 4 3 4 4 3 2 3
u1*u3 *v1*( - 32*a11 *a22*a33 *b11 + 32*a11 *a33 *b11 + 64*a11 *a22 *a33 *b11
3 5 2 3 3
- 64*a11 *a33 *b11 - 32*a11 *a22 *a33 *b11
2 2 4 2 5 2 6
- 96*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11 + 32*a11 *a33 *b11
3 4 6 3 5
+ 64*a11*a22 *a33 *b11 - 64*a11*a22*a33 *b11 - 32*a22 *a33 *b11
2 6 4 2 2
+ 32*a22 *a33 *b11) + u1*u3*v1*v3*( - 32*a11 *a22*a33 *b11
4 3 2 3 2 2 2 3 4 2
+ 32*a11 *a33 *b11 + 64*a11 *a22 *a33 *b11 - 64*a11 *a33 *b11
2 3 2 2 2 2 3 2 2 4 2
- 32*a11 *a22 *a33 *b11 - 96*a11 *a22 *a33 *b11 + 96*a11 *a22*a33 *b11
2 5 2 3 3 2 5 2
+ 32*a11 *a33 *b11 + 64*a11*a22 *a33 *b11 - 64*a11*a22*a33 *b11
3 4 2 2 5 2 3 4 3
- 32*a22 *a33 *b11 + 32*a22 *a33 *b11 ) + u1*v1 *( - 8*a11 *a22*a33*b11
4 2 3 3 2 3 3 3 3
+ 8*a11 *a33 *b11 + 32*a11 *a22*a33 *b11 - 32*a11 *a33 *b11
2 3 3 2 2 2 3 2 3 3
+ 8*a11 *a22 *a33*b11 - 24*a11 *a22 *a33 *b11 - 24*a11 *a22*a33 *b11
2 4 3 3 2 3 2 3 3
+ 40*a11 *a33 *b11 - 16*a11*a22 *a33 *b11 + 48*a11*a22 *a33 *b11
4 3 5 3 3 3 3
- 16*a11*a22*a33 *b11 - 16*a11*a33 *b11 + 8*a22 *a33 *b11
2 4 3 5 3 2
- 24*a22 *a33 *b11 + 16*a22*a33 *b11 ) + u1*v1*v2 *(
4 3 4 2 3 3 2 3
- 8*a11 *a22*a33*b11 + 8*a11 *a33 *b11 - 8*a11 *a22 *a33*b11
3 2 3 3 3 3 2 3 3
+ 48*a11 *a22*a33 *b11 - 40*a11 *a33 *b11 + 16*a11 *a22 *a33*b11
2 2 2 3 2 3 3 2 4 3
- 24*a11 *a22 *a33 *b11 - 48*a11 *a22*a33 *b11 + 56*a11 *a33 *b11
3 2 3 2 3 3 4 3
- 32*a11*a22 *a33 *b11 + 72*a11*a22 *a33 *b11 - 16*a11*a22*a33 *b11
5 3 3 3 3 2 4 3
- 24*a11*a33 *b11 + 16*a22 *a33 *b11 - 40*a22 *a33 *b11
5 3 2 4 3
+ 24*a22*a33 *b11 ) + u1*v1*v3 *( - 16*a11 *a22*a33*b11
4 2 3 3 2 3 3 2 3
+ 16*a11 *a33 *b11 + 8*a11 *a22 *a33*b11 + 48*a11 *a22*a33 *b11
3 3 3 2 3 3 2 2 2 3
- 56*a11 *a33 *b11 + 8*a11 *a22 *a33*b11 - 48*a11 *a22 *a33 *b11
2 3 3 2 4 3 3 2 3
- 24*a11 *a22*a33 *b11 + 64*a11 *a33 *b11 - 16*a11*a22 *a33 *b11
2 3 3 4 3 5 3
+ 72*a11*a22 *a33 *b11 - 32*a11*a22*a33 *b11 - 24*a11*a33 *b11
3 3 3 2 4 3 5 3 4
+ 8*a22 *a33 *b11 - 32*a22 *a33 *b11 + 24*a22*a33 *b11 ) + u2 *(
4 4 4 5 3 2 4 3 5
16*a11 *a22*a33 - 16*a11 *a33 - 64*a11 *a22 *a33 + 64*a11 *a22*a33
2 3 4 2 2 5 4 4
+ 96*a11 *a22 *a33 - 96*a11 *a22 *a33 - 64*a11*a22 *a33
3 5 5 4 4 5 2
+ 64*a11*a22 *a33 + 16*a22 *a33 - 16*a22 *a33 ) + u2 *u3*v3*(
4 3 4 4 3 2 3
- 32*a11 *a22*a33 *b11 + 32*a11 *a33 *b11 + 96*a11 *a22 *a33 *b11
3 4 3 5 2 3 3
- 64*a11 *a22*a33 *b11 - 32*a11 *a33 *b11 - 96*a11 *a22 *a33 *b11
2 5 4 3 3 4
+ 96*a11 *a22*a33 *b11 + 32*a11*a22 *a33 *b11 + 64*a11*a22 *a33 *b11
2 5 4 4 3 5 2 2
- 96*a11*a22 *a33 *b11 - 32*a22 *a33 *b11 + 32*a22 *a33 *b11) + u2 *v1 *(
4 2 2 4 3 2 3 2 2 2
- 8*a11 *a22*a33 *b11 + 8*a11 *a33 *b11 + 24*a11 *a22 *a33 *b11
3 3 2 3 4 2 2 3 2 2
- 16*a11 *a22*a33 *b11 - 8*a11 *a33 *b11 - 24*a11 *a22 *a33 *b11
2 4 2 4 2 2 3 3 2
+ 24*a11 *a22*a33 *b11 + 8*a11*a22 *a33 *b11 + 16*a11*a22 *a33 *b11
2 4 2 4 3 2 3 4 2 2 2
- 24*a11*a22 *a33 *b11 - 8*a22 *a33 *b11 + 8*a22 *a33 *b11 ) + u2 *v3 *
4 2 2 4 3 2 3 2 2 2
( - 8*a11 *a22*a33 *b11 + 8*a11 *a33 *b11 + 16*a11 *a22 *a33 *b11
3 4 2 2 3 2 2 2 2 3 2
- 16*a11 *a33 *b11 - 8*a11 *a22 *a33 *b11 - 24*a11 *a22 *a33 *b11
2 4 2 2 5 2 3 3 2
+ 24*a11 *a22*a33 *b11 + 8*a11 *a33 *b11 + 16*a11*a22 *a33 *b11
5 2 3 4 2 2 5 2 4
- 16*a11*a22*a33 *b11 - 8*a22 *a33 *b11 + 8*a22 *a33 *b11 ) + u3 *(
4 4 4 5 3 2 4 3 6
- 16*a11 *a22*a33 + 16*a11 *a33 + 32*a11 *a22 *a33 - 32*a11 *a33
2 3 4 2 2 5 2 6 2 7
- 16*a11 *a22 *a33 - 48*a11 *a22 *a33 + 48*a11 *a22*a33 + 16*a11 *a33
3 5 7 3 6 2 7 3
+ 32*a11*a22 *a33 - 32*a11*a22*a33 - 16*a22 *a33 + 16*a22 *a33 ) + u3
4 3 4 4 3 2 3
*v3*( - 32*a11 *a22*a33 *b11 + 32*a11 *a33 *b11 + 64*a11 *a22 *a33 *b11
3 5 2 3 3 2 2 4
- 64*a11 *a33 *b11 - 32*a11 *a22 *a33 *b11 - 96*a11 *a22 *a33 *b11
2 5 2 6 3 4
+ 96*a11 *a22*a33 *b11 + 32*a11 *a33 *b11 + 64*a11*a22 *a33 *b11
6 3 5 2 6 2 2
- 64*a11*a22*a33 *b11 - 32*a22 *a33 *b11 + 32*a22 *a33 *b11) + u3 *v1 *
4 2 2 4 3 2 3 3 2
( - 8*a11 *a22*a33 *b11 + 8*a11 *a33 *b11 + 32*a11 *a22*a33 *b11
3 4 2 2 3 2 2 2 2 3 2
- 32*a11 *a33 *b11 + 8*a11 *a22 *a33 *b11 - 24*a11 *a22 *a33 *b11
2 4 2 2 5 2 3 3 2
- 24*a11 *a22*a33 *b11 + 40*a11 *a33 *b11 - 16*a11*a22 *a33 *b11
2 4 2 5 2 6 2
+ 48*a11*a22 *a33 *b11 - 16*a11*a22*a33 *b11 - 16*a11*a33 *b11
3 4 2 2 5 2 6 2 2 2
+ 8*a22 *a33 *b11 - 24*a22 *a33 *b11 + 16*a22*a33 *b11 ) + u3 *v2 *(
4 2 2 4 3 2 3 2 2 2
- 8*a11 *a22*a33 *b11 + 8*a11 *a33 *b11 - 8*a11 *a22 *a33 *b11
3 3 2 3 4 2 2 3 2 2
+ 48*a11 *a22*a33 *b11 - 40*a11 *a33 *b11 + 16*a11 *a22 *a33 *b11
2 2 3 2 2 4 2 2 5 2
- 24*a11 *a22 *a33 *b11 - 48*a11 *a22*a33 *b11 + 56*a11 *a33 *b11
3 3 2 2 4 2 5 2
- 32*a11*a22 *a33 *b11 + 72*a11*a22 *a33 *b11 - 16*a11*a22*a33 *b11
6 2 3 4 2 2 5 2
- 24*a11*a33 *b11 + 16*a22 *a33 *b11 - 40*a22 *a33 *b11
6 2 2 2 4 2 2
+ 24*a22*a33 *b11 ) + u3 *v3 *( - 16*a11 *a22*a33 *b11
4 3 2 3 2 2 2 3 3 2
+ 16*a11 *a33 *b11 + 8*a11 *a22 *a33 *b11 + 48*a11 *a22*a33 *b11
3 4 2 2 3 2 2 2 2 3 2
- 56*a11 *a33 *b11 + 8*a11 *a22 *a33 *b11 - 48*a11 *a22 *a33 *b11
2 4 2 2 5 2 3 3 2
- 24*a11 *a22*a33 *b11 + 64*a11 *a33 *b11 - 16*a11*a22 *a33 *b11
2 4 2 5 2 6 2
+ 72*a11*a22 *a33 *b11 - 32*a11*a22*a33 *b11 - 24*a11*a33 *b11
3 4 2 2 5 2 6 2 2
+ 8*a22 *a33 *b11 - 32*a22 *a33 *b11 + 24*a22*a33 *b11 ) + u3*v1 *v3*(
3 2 3 3 2 3 3 3 3
- 16*a11 *a22 *a33*b11 + 32*a11 *a22*a33 *b11 - 16*a11 *a33 *b11
2 3 3 2 3 3 2 4 3
+ 16*a11 *a22 *a33*b11 - 48*a11 *a22*a33 *b11 + 32*a11 *a33 *b11
3 2 3 2 3 3 5 3
- 32*a11*a22 *a33 *b11 + 48*a11*a22 *a33 *b11 - 16*a11*a33 *b11
3 3 3 2 4 3 5 3 2
+ 16*a22 *a33 *b11 - 32*a22 *a33 *b11 + 16*a22*a33 *b11 ) + u3*v2 *v3*(
4 3 4 2 3 3 2 3
- 8*a11 *a22*a33*b11 + 8*a11 *a33 *b11 - 8*a11 *a22 *a33*b11
3 2 3 3 3 3 2 3 3
+ 48*a11 *a22*a33 *b11 - 40*a11 *a33 *b11 + 16*a11 *a22 *a33*b11
2 2 2 3 2 3 3 2 4 3
- 24*a11 *a22 *a33 *b11 - 48*a11 *a22*a33 *b11 + 56*a11 *a33 *b11
3 2 3 2 3 3 4 3
- 32*a11*a22 *a33 *b11 + 72*a11*a22 *a33 *b11 - 16*a11*a22*a33 *b11
5 3 3 3 3 2 4 3
- 24*a11*a33 *b11 + 16*a22 *a33 *b11 - 40*a22 *a33 *b11
5 3 3 4 3 4 2 3
+ 24*a22*a33 *b11 ) + u3*v3 *( - 8*a11 *a22*a33*b11 + 8*a11 *a33 *b11
3 2 3 3 3 3 2 3 3
+ 32*a11 *a22*a33 *b11 - 32*a11 *a33 *b11 + 8*a11 *a22 *a33*b11
2 2 2 3 2 3 3 2 4 3
- 24*a11 *a22 *a33 *b11 - 24*a11 *a22*a33 *b11 + 40*a11 *a33 *b11
3 2 3 2 3 3 4 3
- 16*a11*a22 *a33 *b11 + 48*a11*a22 *a33 *b11 - 16*a11*a22*a33 *b11
5 3 3 3 3 2 4 3
- 16*a11*a33 *b11 + 8*a22 *a33 *b11 - 24*a22 *a33 *b11
5 3 2 2 4 4 4 4
+ 16*a22*a33 *b11 ) + v1 *v2 *( - 2*a11 *a22*b11 + 2*a11 *a33*b11
3 2 4 3 4 3 2 4
+ 2*a11 *a22 *b11 + 4*a11 *a22*a33*b11 - 6*a11 *a33 *b11
2 3 4 2 2 4 2 2 4
- 4*a11 *a22 *b11 + 6*a11 *a22 *a33*b11 - 12*a11 *a22*a33 *b11
2 3 4 3 4 2 2 4
+ 10*a11 *a33 *b11 + 8*a11*a22 *a33*b11 - 18*a11*a22 *a33 *b11
3 4 4 4 3 2 4
+ 20*a11*a22*a33 *b11 - 10*a11*a33 *b11 - 4*a22 *a33 *b11
2 3 4 4 4 5 4 2 2
+ 10*a22 *a33 *b11 - 10*a22*a33 *b11 + 4*a33 *b11 ) + v1 *v3 *(
4 4 4 4 3 4
- 2*a11 *a22*b11 + 2*a11 *a33*b11 + 8*a11 *a22*a33*b11
3 2 4 2 3 4 2 2 4
- 8*a11 *a33 *b11 - 2*a11 *a22 *b11 + 6*a11 *a22 *a33*b11
2 2 4 2 3 4 3 4
- 18*a11 *a22*a33 *b11 + 14*a11 *a33 *b11 + 4*a11*a22 *a33*b11
2 2 4 3 4 4 4
- 12*a11*a22 *a33 *b11 + 20*a11*a22*a33 *b11 - 12*a11*a33 *b11
3 2 4 2 3 4 4 4 5 4
- 2*a22 *a33 *b11 + 6*a22 *a33 *b11 - 8*a22*a33 *b11 + 4*a33 *b11 ) +
4 4 4 4 4 3 4
v2 *( - a11 *a22*b11 + a11 *a33*b11 + 4*a11 *a22*a33*b11
3 2 4 2 3 4 2 2 4
- 4*a11 *a33 *b11 - 4*a11 *a22 *b11 + 12*a11 *a22 *a33*b11
2 2 4 2 3 4 3 4
- 18*a11 *a22*a33 *b11 + 10*a11 *a33 *b11 + 8*a11*a22 *a33*b11
2 2 4 3 4 4 4
- 24*a11*a22 *a33 *b11 + 28*a11*a22*a33 *b11 - 12*a11*a33 *b11
3 2 4 2 3 4 4 4 5 4
- 4*a22 *a33 *b11 + 12*a22 *a33 *b11 - 13*a22*a33 *b11 + 5*a33 *b11
2 2 4 4 4 4 3 2 4
) + v2 *v3 *( - 4*a11 *a22*b11 + 4*a11 *a33*b11 - 2*a11 *a22 *b11
3 4 3 2 4 2 3 4
+ 20*a11 *a22*a33*b11 - 18*a11 *a33 *b11 - 4*a11 *a22 *b11
2 2 4 2 2 4 2 3 4
+ 18*a11 *a22 *a33*b11 - 48*a11 *a22*a33 *b11 + 34*a11 *a33 *b11
3 4 2 2 4 3 4
+ 8*a11*a22 *a33*b11 - 30*a11*a22 *a33 *b11 + 52*a11*a22*a33 *b11
4 4 3 2 4 2 3 4
- 30*a11*a33 *b11 - 4*a22 *a33 *b11 + 14*a22 *a33 *b11
4 4 5 4 4 4 4
- 20*a22*a33 *b11 + 10*a33 *b11 ) + v3 *( - 3*a11 *a22*b11
4 4 3 4 3 2 4
+ 3*a11 *a33*b11 + 12*a11 *a22*a33*b11 - 12*a11 *a33 *b11
2 3 4 2 2 4 2 2 4
- a11 *a22 *b11 + 3*a11 *a22 *a33*b11 - 21*a11 *a22*a33 *b11
2 3 4 3 4 2 2 4
+ 19*a11 *a33 *b11 + 2*a11*a22 *a33*b11 - 6*a11*a22 *a33 *b11
3 4 4 4 3 2 4
+ 18*a11*a22*a33 *b11 - 14*a11*a33 *b11 - a22 *a33 *b11
2 3 4 4 4 5 4
+ 3*a22 *a33 *b11 - 6*a22*a33 *b11 + 4*a33 *b11 )
Stop of a subroutine.
Number of garbage collections exceeds max_gc_fac.
The factorization crashed!
which the program can not factorize further.
{HAM,FI} = 0
5 4 4 2
FI=u1*v1*(4*a11 *a33*b11 - 16*a11 *a22*a33*b11 - 4*a11 *a33 *b11
3 2 3 2 2 3
+ 24*a11 *a22 *a33*b11 + 16*a11 *a22*a33 *b11 - 16*a11 *a22 *a33*b11
2 2 2 4 3 2
- 24*a11 *a22 *a33 *b11 + 4*a11*a22 *a33*b11 + 16*a11*a22 *a33 *b11
4 2 2 5 2 4 2
- 4*a22 *a33 *b11) + u2 *(4*a11 *a33 - 20*a11 *a22*a33
3 2 2 2 3 2 4 2 5 2
+ 40*a11 *a22 *a33 - 40*a11 *a22 *a33 + 20*a11*a22 *a33 - 4*a22 *a33 )
2 5 2 4 2 4 3 3 2 2
+ u3 *(4*a11 *a33 - 16*a11 *a22*a33 - 4*a11 *a33 + 24*a11 *a22 *a33
3 3 2 3 2 2 2 3
+ 16*a11 *a22*a33 - 16*a11 *a22 *a33 - 24*a11 *a22 *a33
4 2 3 3 4 3 2 5 2
+ 4*a11*a22 *a33 + 16*a11*a22 *a33 - 4*a22 *a33 ) + v2 *(a11 *b11
4 2 4 2 3 2 2
- 3*a11 *a22*b11 - 2*a11 *a33*b11 + 3*a11 *a22 *b11
3 2 3 2 2 2 3 2
+ 6*a11 *a22*a33*b11 + a11 *a33 *b11 - a11 *a22 *b11
2 2 2 2 2 2 3 2
- 6*a11 *a22 *a33*b11 - 3*a11 *a22*a33 *b11 + 2*a11*a22 *a33*b11
2 2 2 3 2 2 2 5 2
+ 3*a11*a22 *a33 *b11 - a22 *a33 *b11 ) + v3 *(a11 *b11
4 2 4 2 3 2 2
- 4*a11 *a22*b11 - a11 *a33*b11 + 6*a11 *a22 *b11
3 2 2 3 2 2 2 2
+ 4*a11 *a22*a33*b11 - 4*a11 *a22 *b11 - 6*a11 *a22 *a33*b11
4 2 3 2 4 2
+ a11*a22 *b11 + 4*a11*a22 *a33*b11 - a22 *a33*b11 )
= a product of the elements of: {4,
a11 - a22,
a11 - a22,
a11 - a22,
2 2 2
u1*v1*(a11 *a33*b11 - a11*a22*a33*b11 - a11*a33 *b11 + a22*a33 *b11)
2 2 2 2 2 2
+ u2 *(a11 *a33 - 2*a11*a22*a33 + a22 *a33 )
2 2 2 2 3 3
+ u3 *(a11 *a33 - a11*a22*a33 - a11*a33 + a22*a33 )
2 2 2 2 2
2 a11 *b11 - 2*a11*a33*b11 + a33 *b11
+ v2 *----------------------------------------
4
2 2 2 2 2
2 a11 *b11 - a11*a22*b11 - a11*a33*b11 + a22*a33*b11
+ v3 *--------------------------------------------------------}
4
{HAM,FI} = 0
And again in machine readable form:
HAM=u1**2*a11 + u1*v1*b11 + u2**2*a22 + u3**2*a33 + u3*v3*( - a11*a22*b11 + a11*
a33*b11)/(a11*a33 - a22*a33) + v2**2*(a11**2*a22*b11**2 - a11*a22**2*b11**2 -
a11*a33**2*b11**2 + a22*a33**2*b11**2)/(4*a11**2*a33**2 - 8*a11*a22*a33**2 + 4*
a22**2*a33**2) + v3**2*(a11**2*a33*b11**2 - a11*a22**2*b11**2 - a11*a33**2*b11**
2 + a22**2*a33*b11**2)/(4*a11**2*a33**2 - 8*a11*a22*a33**2 + 4*a22**2*a33**2)$
FI=u1**2*v2**2*(16*a11**5*a33**2*b11**2 - 48*a11**4*a22*a33**2*b11**2 - 32*a11**
4*a33**3*b11**2 + 48*a11**3*a22**2*a33**2*b11**2 + 96*a11**3*a22*a33**3*b11**2 +
16*a11**3*a33**4*b11**2 - 16*a11**2*a22**3*a33**2*b11**2 - 96*a11**2*a22**2*a33
**3*b11**2 - 48*a11**2*a22*a33**4*b11**2 + 32*a11*a22**3*a33**3*b11**2 + 48*a11*
a22**2*a33**4*b11**2 - 16*a22**3*a33**4*b11**2) + u1**2*v3**2*(16*a11**5*a33**2*
b11**2 - 48*a11**4*a22*a33**2*b11**2 - 32*a11**4*a33**3*b11**2 + 48*a11**3*a22**
2*a33**2*b11**2 + 96*a11**3*a22*a33**3*b11**2 + 16*a11**3*a33**4*b11**2 - 16*a11
**2*a22**3*a33**2*b11**2 - 96*a11**2*a22**2*a33**3*b11**2 - 48*a11**2*a22*a33**4
*b11**2 + 32*a11*a22**3*a33**3*b11**2 + 48*a11*a22**2*a33**4*b11**2 - 16*a22**3*
a33**4*b11**2) + u1*u2*v1*v2*( - 32*a11**5*a33**2*b11**2 + 96*a11**4*a22*a33**2*
b11**2 + 64*a11**4*a33**3*b11**2 - 96*a11**3*a22**2*a33**2*b11**2 - 192*a11**3*
a22*a33**3*b11**2 - 32*a11**3*a33**4*b11**2 + 32*a11**2*a22**3*a33**2*b11**2 +
192*a11**2*a22**2*a33**3*b11**2 + 96*a11**2*a22*a33**4*b11**2 - 64*a11*a22**3*
a33**3*b11**2 - 96*a11*a22**2*a33**4*b11**2 + 32*a22**3*a33**4*b11**2) + u1*u3*
v1*v3*( - 32*a11**5*a33**2*b11**2 + 96*a11**4*a22*a33**2*b11**2 + 64*a11**4*a33
**3*b11**2 - 96*a11**3*a22**2*a33**2*b11**2 - 192*a11**3*a22*a33**3*b11**2 - 32*
a11**3*a33**4*b11**2 + 32*a11**2*a22**3*a33**2*b11**2 + 192*a11**2*a22**2*a33**3
*b11**2 + 96*a11**2*a22*a33**4*b11**2 - 64*a11*a22**3*a33**3*b11**2 - 96*a11*a22
**2*a33**4*b11**2 + 32*a22**3*a33**4*b11**2) + u1*v1**3*( - 8*a11**5*a33*b11**3
+ 40*a11**4*a33**2*b11**3 + 24*a11**3*a22**2*a33*b11**3 - 48*a11**3*a22*a33**2*
b11**3 - 56*a11**3*a33**3*b11**3 - 16*a11**2*a22**3*a33*b11**3 - 24*a11**2*a22**
2*a33**2*b11**3 + 96*a11**2*a22*a33**3*b11**3 + 24*a11**2*a33**4*b11**3 + 32*a11
*a22**3*a33**2*b11**3 - 24*a11*a22**2*a33**3*b11**3 - 48*a11*a22*a33**4*b11**3 -
16*a22**3*a33**3*b11**3 + 24*a22**2*a33**4*b11**3) + u1*v1*v2**2*( - 8*a11**5*
a33*b11**3 + 40*a11**4*a33**2*b11**3 + 24*a11**3*a22**2*a33*b11**3 - 48*a11**3*
a22*a33**2*b11**3 - 56*a11**3*a33**3*b11**3 - 16*a11**2*a22**3*a33*b11**3 - 24*
a11**2*a22**2*a33**2*b11**3 + 96*a11**2*a22*a33**3*b11**3 + 24*a11**2*a33**4*b11
**3 + 32*a11*a22**3*a33**2*b11**3 - 24*a11*a22**2*a33**3*b11**3 - 48*a11*a22*a33
**4*b11**3 - 16*a22**3*a33**3*b11**3 + 24*a22**2*a33**4*b11**3) + u1*v1*v3**2*(
- 8*a11**5*a33*b11**3 + 40*a11**4*a33**2*b11**3 + 24*a11**3*a22**2*a33*b11**3 -
48*a11**3*a22*a33**2*b11**3 - 56*a11**3*a33**3*b11**3 - 16*a11**2*a22**3*a33*b11
**3 - 24*a11**2*a22**2*a33**2*b11**3 + 96*a11**2*a22*a33**3*b11**3 + 24*a11**2*
a33**4*b11**3 + 32*a11*a22**3*a33**2*b11**3 - 24*a11*a22**2*a33**3*b11**3 - 48*
a11*a22*a33**4*b11**3 - 16*a22**3*a33**3*b11**3 + 24*a22**2*a33**4*b11**3) + u2
**4*(16*a11**5*a33**4 - 80*a11**4*a22*a33**4 + 160*a11**3*a22**2*a33**4 - 160*
a11**2*a22**3*a33**4 + 80*a11*a22**4*a33**4 - 16*a22**5*a33**4) + u2**3*v2*( -
32*a11**5*a33**3*b11 + 128*a11**4*a22*a33**3*b11 + 32*a11**4*a33**4*b11 - 192*
a11**3*a22**2*a33**3*b11 - 128*a11**3*a22*a33**4*b11 + 128*a11**2*a22**3*a33**3*
b11 + 192*a11**2*a22**2*a33**4*b11 - 32*a11*a22**4*a33**3*b11 - 128*a11*a22**3*
a33**4*b11 + 32*a22**4*a33**4*b11) + u2**2*u3**2*(32*a11**5*a33**4 - 128*a11**4*
a22*a33**4 - 32*a11**4*a33**5 + 192*a11**3*a22**2*a33**4 + 128*a11**3*a22*a33**5
- 128*a11**2*a22**3*a33**4 - 192*a11**2*a22**2*a33**5 + 32*a11*a22**4*a33**4 +
128*a11*a22**3*a33**5 - 32*a22**4*a33**5) + u2**2*u3*v3*( - 32*a11**5*a33**3*b11
+ 128*a11**4*a22*a33**3*b11 + 32*a11**4*a33**4*b11 - 192*a11**3*a22**2*a33**3*
b11 - 128*a11**3*a22*a33**4*b11 + 128*a11**2*a22**3*a33**3*b11 + 192*a11**2*a22
**2*a33**4*b11 - 32*a11*a22**4*a33**3*b11 - 128*a11*a22**3*a33**4*b11 + 32*a22**
4*a33**4*b11) + u2**2*v1**2*( - 8*a11**5*a33**2*b11**2 + 24*a11**4*a22*a33**2*
b11**2 + 16*a11**4*a33**3*b11**2 - 24*a11**3*a22**2*a33**2*b11**2 - 48*a11**3*
a22*a33**3*b11**2 - 8*a11**3*a33**4*b11**2 + 8*a11**2*a22**3*a33**2*b11**2 + 48*
a11**2*a22**2*a33**3*b11**2 + 24*a11**2*a22*a33**4*b11**2 - 16*a11*a22**3*a33**3
*b11**2 - 24*a11*a22**2*a33**4*b11**2 + 8*a22**3*a33**4*b11**2) + u2**2*v3**2*(
- 8*a11**4*a22*a33**2*b11**2 + 8*a11**4*a33**3*b11**2 + 24*a11**3*a22**2*a33**2*
b11**2 - 16*a11**3*a22*a33**3*b11**2 - 8*a11**3*a33**4*b11**2 - 24*a11**2*a22**3
*a33**2*b11**2 + 24*a11**2*a22*a33**4*b11**2 + 8*a11*a22**4*a33**2*b11**2 + 16*
a11*a22**3*a33**3*b11**2 - 24*a11*a22**2*a33**4*b11**2 - 8*a22**4*a33**3*b11**2
+ 8*a22**3*a33**4*b11**2) + u2*u3**2*v2*( - 32*a11**5*a33**3*b11 + 96*a11**4*a22
*a33**3*b11 + 64*a11**4*a33**4*b11 - 96*a11**3*a22**2*a33**3*b11 - 192*a11**3*
a22*a33**4*b11 - 32*a11**3*a33**5*b11 + 32*a11**2*a22**3*a33**3*b11 + 192*a11**2
*a22**2*a33**4*b11 + 96*a11**2*a22*a33**5*b11 - 64*a11*a22**3*a33**4*b11 - 96*
a11*a22**2*a33**5*b11 + 32*a22**3*a33**5*b11) + u2*v2**3*( - 8*a11**5*a33*b11**3
+ 16*a11**4*a22*a33*b11**3 + 24*a11**4*a33**2*b11**3 - 8*a11**3*a22**2*a33*b11
**3 - 48*a11**3*a22*a33**2*b11**3 - 24*a11**3*a33**3*b11**3 + 24*a11**2*a22**2*
a33**2*b11**3 + 48*a11**2*a22*a33**3*b11**3 + 8*a11**2*a33**4*b11**3 - 24*a11*
a22**2*a33**3*b11**3 - 16*a11*a22*a33**4*b11**3 + 8*a22**2*a33**4*b11**3) + u2*
v2*v3**2*( - 8*a11**5*a33*b11**3 + 24*a11**4*a22*a33*b11**3 + 16*a11**4*a33**2*
b11**3 - 24*a11**3*a22**2*a33*b11**3 - 48*a11**3*a22*a33**2*b11**3 - 8*a11**3*
a33**3*b11**3 + 8*a11**2*a22**3*a33*b11**3 + 48*a11**2*a22**2*a33**2*b11**3 + 24
*a11**2*a22*a33**3*b11**3 - 16*a11*a22**3*a33**2*b11**3 - 24*a11*a22**2*a33**3*
b11**3 + 8*a22**3*a33**3*b11**3) + u3**4*(16*a11**5*a33**4 - 48*a11**4*a22*a33**
4 - 32*a11**4*a33**5 + 48*a11**3*a22**2*a33**4 + 96*a11**3*a22*a33**5 + 16*a11**
3*a33**6 - 16*a11**2*a22**3*a33**4 - 96*a11**2*a22**2*a33**5 - 48*a11**2*a22*a33
**6 + 32*a11*a22**3*a33**5 + 48*a11*a22**2*a33**6 - 16*a22**3*a33**6) + u3**3*v3
*( - 32*a11**5*a33**3*b11 + 96*a11**4*a22*a33**3*b11 + 64*a11**4*a33**4*b11 - 96
*a11**3*a22**2*a33**3*b11 - 192*a11**3*a22*a33**4*b11 - 32*a11**3*a33**5*b11 +
32*a11**2*a22**3*a33**3*b11 + 192*a11**2*a22**2*a33**4*b11 + 96*a11**2*a22*a33**
5*b11 - 64*a11*a22**3*a33**4*b11 - 96*a11*a22**2*a33**5*b11 + 32*a22**3*a33**5*
b11) + u3**2*v1**2*( - 8*a11**5*a33**2*b11**2 + 40*a11**4*a33**3*b11**2 + 24*a11
**3*a22**2*a33**2*b11**2 - 48*a11**3*a22*a33**3*b11**2 - 56*a11**3*a33**4*b11**2
- 16*a11**2*a22**3*a33**2*b11**2 - 24*a11**2*a22**2*a33**3*b11**2 + 96*a11**2*
a22*a33**4*b11**2 + 24*a11**2*a33**5*b11**2 + 32*a11*a22**3*a33**3*b11**2 - 24*
a11*a22**2*a33**4*b11**2 - 48*a11*a22*a33**5*b11**2 - 16*a22**3*a33**4*b11**2 +
24*a22**2*a33**5*b11**2) + u3**2*v2**2*( - 16*a11**4*a22*a33**2*b11**2 + 16*a11
**4*a33**3*b11**2 + 32*a11**3*a22**2*a33**2*b11**2 - 32*a11**3*a33**4*b11**2 -
16*a11**2*a22**3*a33**2*b11**2 - 48*a11**2*a22**2*a33**3*b11**2 + 48*a11**2*a22*
a33**4*b11**2 + 16*a11**2*a33**5*b11**2 + 32*a11*a22**3*a33**3*b11**2 - 32*a11*
a22*a33**5*b11**2 - 16*a22**3*a33**4*b11**2 + 16*a22**2*a33**5*b11**2) + u3**2*
v3**2*( - 24*a11**4*a22*a33**2*b11**2 + 24*a11**4*a33**3*b11**2 + 48*a11**3*a22
**2*a33**2*b11**2 - 48*a11**3*a33**4*b11**2 - 24*a11**2*a22**3*a33**2*b11**2 -
72*a11**2*a22**2*a33**3*b11**2 + 72*a11**2*a22*a33**4*b11**2 + 24*a11**2*a33**5*
b11**2 + 48*a11*a22**3*a33**3*b11**2 - 48*a11*a22*a33**5*b11**2 - 24*a22**3*a33
**4*b11**2 + 24*a22**2*a33**5*b11**2) + u3*v1**2*v3*( - 16*a11**4*a22*a33*b11**3
+ 16*a11**4*a33**2*b11**3 + 32*a11**3*a22**2*a33*b11**3 - 32*a11**3*a33**3*b11
**3 - 16*a11**2*a22**3*a33*b11**3 - 48*a11**2*a22**2*a33**2*b11**3 + 48*a11**2*
a22*a33**3*b11**3 + 16*a11**2*a33**4*b11**3 + 32*a11*a22**3*a33**2*b11**3 - 32*
a11*a22*a33**4*b11**3 - 16*a22**3*a33**3*b11**3 + 16*a22**2*a33**4*b11**3) + u3*
v2**2*v3*( - 8*a11**5*a33*b11**3 + 40*a11**4*a33**2*b11**3 + 24*a11**3*a22**2*
a33*b11**3 - 48*a11**3*a22*a33**2*b11**3 - 56*a11**3*a33**3*b11**3 - 16*a11**2*
a22**3*a33*b11**3 - 24*a11**2*a22**2*a33**2*b11**3 + 96*a11**2*a22*a33**3*b11**3
+ 24*a11**2*a33**4*b11**3 + 32*a11*a22**3*a33**2*b11**3 - 24*a11*a22**2*a33**3*
b11**3 - 48*a11*a22*a33**4*b11**3 - 16*a22**3*a33**3*b11**3 + 24*a22**2*a33**4*
b11**3) + u3*v3**3*( - 8*a11**5*a33*b11**3 + 8*a11**4*a22*a33*b11**3 + 32*a11**4
*a33**2*b11**3 + 8*a11**3*a22**2*a33*b11**3 - 48*a11**3*a22*a33**2*b11**3 - 40*
a11**3*a33**3*b11**3 - 8*a11**2*a22**3*a33*b11**3 + 72*a11**2*a22*a33**3*b11**3
+ 16*a11**2*a33**4*b11**3 + 16*a11*a22**3*a33**2*b11**3 - 24*a11*a22**2*a33**3*
b11**3 - 32*a11*a22*a33**4*b11**3 - 8*a22**3*a33**3*b11**3 + 16*a22**2*a33**4*
b11**3) + v1**2*v2**2*( - 2*a11**5*b11**4 + 2*a11**4*a22*b11**4 + 8*a11**4*a33*
b11**4 - 4*a11**3*a22**2*b11**4 - 16*a11**3*a33**2*b11**4 + 4*a11**2*a22**3*b11
**4 + 16*a11**2*a33**3*b11**4 - 8*a11*a22**3*a33*b11**4 + 12*a11*a22**2*a33**2*
b11**4 - 8*a11*a22*a33**3*b11**4 - 6*a11*a33**4*b11**4 + 4*a22**3*a33**2*b11**4
- 8*a22**2*a33**3*b11**4 + 6*a22*a33**4*b11**4) + v1**2*v3**2*( - 2*a11**5*b11**
4 + 10*a11**4*a33*b11**4 + 2*a11**3*a22**2*b11**4 - 4*a11**3*a22*a33*b11**4 - 18
*a11**3*a33**2*b11**4 - 6*a11**2*a22**2*a33*b11**4 + 12*a11**2*a22*a33**2*b11**4
+ 14*a11**2*a33**3*b11**4 + 6*a11*a22**2*a33**2*b11**4 - 12*a11*a22*a33**3*b11
**4 - 4*a11*a33**4*b11**4 - 2*a22**2*a33**3*b11**4 + 4*a22*a33**4*b11**4) + v2**
4*( - a11**5*b11**4 + a11**4*a22*b11**4 + 4*a11**4*a33*b11**4 - 4*a11**3*a22**2*
b11**4 + 4*a11**3*a22*a33*b11**4 - 10*a11**3*a33**2*b11**4 + 4*a11**2*a22**3*b11
**4 - 6*a11**2*a22*a33**2*b11**4 + 12*a11**2*a33**3*b11**4 - 8*a11*a22**3*a33*
b11**4 + 12*a11*a22**2*a33**2*b11**4 - 4*a11*a22*a33**3*b11**4 - 5*a11*a33**4*
b11**4 + 4*a22**3*a33**2*b11**4 - 8*a22**2*a33**3*b11**4 + 5*a22*a33**4*b11**4)
+ v2**2*v3**2*( - 2*a11**5*b11**4 - 2*a11**4*a22*b11**4 + 12*a11**4*a33*b11**4 +
8*a11**3*a22*a33*b11**4 - 28*a11**3*a33**2*b11**4 + 4*a11**2*a22**3*b11**4 - 12
*a11**2*a22**2*a33*b11**4 + 28*a11**2*a33**3*b11**4 - 8*a11*a22**3*a33*b11**4 +
24*a11*a22**2*a33**2*b11**4 - 16*a11*a22*a33**3*b11**4 - 10*a11*a33**4*b11**4 +
4*a22**3*a33**2*b11**4 - 12*a22**2*a33**3*b11**4 + 10*a22*a33**4*b11**4) + v3**4
*( - a11**5*b11**4 - 3*a11**4*a22*b11**4 + 8*a11**4*a33*b11**4 + 5*a11**3*a22**2
*b11**4 + 2*a11**3*a22*a33*b11**4 - 17*a11**3*a33**2*b11**4 - a11**2*a22**3*b11
**4 - 12*a11**2*a22**2*a33*b11**4 + 9*a11**2*a22*a33**2*b11**4 + 14*a11**2*a33**
3*b11**4 + 2*a11*a22**3*a33*b11**4 + 9*a11*a22**2*a33**2*b11**4 - 12*a11*a22*a33
**3*b11**4 - 4*a11*a33**4*b11**4 - a22**3*a33**2*b11**4 - 2*a22**2*a33**3*b11**4
+ 4*a22*a33**4*b11**4)$
FI=u1**2*v2**2*( - 4*a11**5*a33**2 + 16*a11**4*a22*a33**2 + 4*a11**4*a33**3 - 24
*a11**3*a22**2*a33**2 - 16*a11**3*a22*a33**3 + 16*a11**2*a22**3*a33**2 + 24*a11
**2*a22**2*a33**3 - 4*a11*a22**4*a33**2 - 16*a11*a22**3*a33**3 + 4*a22**4*a33**3
) + u1**2*v3**2*( - 4*a11**5*a33**2 + 16*a11**4*a22*a33**2 + 4*a11**4*a33**3 -
24*a11**3*a22**2*a33**2 - 16*a11**3*a22*a33**3 + 16*a11**2*a22**3*a33**2 + 24*
a11**2*a22**2*a33**3 - 4*a11*a22**4*a33**2 - 16*a11*a22**3*a33**3 + 4*a22**4*a33
**3) + u1*u2*v1*v2*(8*a11**5*a33**2 - 32*a11**4*a22*a33**2 - 8*a11**4*a33**3 +
48*a11**3*a22**2*a33**2 + 32*a11**3*a22*a33**3 - 32*a11**2*a22**3*a33**2 - 48*
a11**2*a22**2*a33**3 + 8*a11*a22**4*a33**2 + 32*a11*a22**3*a33**3 - 8*a22**4*a33
**3) + u1*u3*v1*v3*(8*a11**5*a33**2 - 32*a11**4*a22*a33**2 - 8*a11**4*a33**3 +
48*a11**3*a22**2*a33**2 + 32*a11**3*a22*a33**3 - 32*a11**2*a22**3*a33**2 - 48*
a11**2*a22**2*a33**3 + 8*a11*a22**4*a33**2 + 32*a11*a22**3*a33**3 - 8*a22**4*a33
**3) + u1*v1**3*( - 4*a11**5*a33*b11 + 16*a11**4*a22*a33*b11 + 4*a11**4*a33**2*
b11 - 24*a11**3*a22**2*a33*b11 - 16*a11**3*a22*a33**2*b11 + 16*a11**2*a22**3*a33
*b11 + 24*a11**2*a22**2*a33**2*b11 - 4*a11*a22**4*a33*b11 - 16*a11*a22**3*a33**2
*b11 + 4*a22**4*a33**2*b11) + u1*v1*v2**2*( - 4*a11**5*a33*b11 + 16*a11**4*a22*
a33*b11 + 4*a11**4*a33**2*b11 - 24*a11**3*a22**2*a33*b11 - 16*a11**3*a22*a33**2*
b11 + 16*a11**2*a22**3*a33*b11 + 24*a11**2*a22**2*a33**2*b11 - 4*a11*a22**4*a33*
b11 - 16*a11*a22**3*a33**2*b11 + 4*a22**4*a33**2*b11) + u1*v1*v3**2*( - 4*a11**5
*a33*b11 + 16*a11**4*a22*a33*b11 + 4*a11**4*a33**2*b11 - 24*a11**3*a22**2*a33*
b11 - 16*a11**3*a22*a33**2*b11 + 16*a11**2*a22**3*a33*b11 + 24*a11**2*a22**2*a33
**2*b11 - 4*a11*a22**4*a33*b11 - 16*a11*a22**3*a33**2*b11 + 4*a22**4*a33**2*b11)
+ u2**2*v1**2*( - 4*a11**5*a33**2 + 16*a11**4*a22*a33**2 + 4*a11**4*a33**3 - 24
*a11**3*a22**2*a33**2 - 16*a11**3*a22*a33**3 + 16*a11**2*a22**3*a33**2 + 24*a11
**2*a22**2*a33**3 - 4*a11*a22**4*a33**2 - 16*a11*a22**3*a33**3 + 4*a22**4*a33**3
) + u2**2*v3**2*( - 4*a11**5*a33**2 + 16*a11**4*a22*a33**2 + 4*a11**4*a33**3 -
24*a11**3*a22**2*a33**2 - 16*a11**3*a22*a33**3 + 16*a11**2*a22**3*a33**2 + 24*
a11**2*a22**2*a33**3 - 4*a11*a22**4*a33**2 - 16*a11*a22**3*a33**3 + 4*a22**4*a33
**3) + u2*u3*v2*v3*(8*a11**5*a33**2 - 32*a11**4*a22*a33**2 - 8*a11**4*a33**3 +
48*a11**3*a22**2*a33**2 + 32*a11**3*a22*a33**3 - 32*a11**2*a22**3*a33**2 - 48*
a11**2*a22**2*a33**3 + 8*a11*a22**4*a33**2 + 32*a11*a22**3*a33**3 - 8*a22**4*a33
**3) + u3**2*v1**2*( - 4*a11**5*a33**2 + 16*a11**4*a22*a33**2 + 4*a11**4*a33**3
- 24*a11**3*a22**2*a33**2 - 16*a11**3*a22*a33**3 + 16*a11**2*a22**3*a33**2 + 24*
a11**2*a22**2*a33**3 - 4*a11*a22**4*a33**2 - 16*a11*a22**3*a33**3 + 4*a22**4*a33
**3) + u3**2*v2**2*( - 4*a11**5*a33**2 + 16*a11**4*a22*a33**2 + 4*a11**4*a33**3
- 24*a11**3*a22**2*a33**2 - 16*a11**3*a22*a33**3 + 16*a11**2*a22**3*a33**2 + 24*
a11**2*a22**2*a33**3 - 4*a11*a22**4*a33**2 - 16*a11*a22**3*a33**3 + 4*a22**4*a33
**3) + u3*v1**2*v3*(4*a11**4*a22*a33*b11 - 4*a11**4*a33**2*b11 - 12*a11**3*a22**
2*a33*b11 + 8*a11**3*a22*a33**2*b11 + 4*a11**3*a33**3*b11 + 12*a11**2*a22**3*a33
*b11 - 12*a11**2*a22*a33**3*b11 - 4*a11*a22**4*a33*b11 - 8*a11*a22**3*a33**2*b11
+ 12*a11*a22**2*a33**3*b11 + 4*a22**4*a33**2*b11 - 4*a22**3*a33**3*b11) + u3*v2
**2*v3*(4*a11**4*a22*a33*b11 - 4*a11**4*a33**2*b11 - 12*a11**3*a22**2*a33*b11 +
8*a11**3*a22*a33**2*b11 + 4*a11**3*a33**3*b11 + 12*a11**2*a22**3*a33*b11 - 12*
a11**2*a22*a33**3*b11 - 4*a11*a22**4*a33*b11 - 8*a11*a22**3*a33**2*b11 + 12*a11*
a22**2*a33**3*b11 + 4*a22**4*a33**2*b11 - 4*a22**3*a33**3*b11) + u3*v3**3*(4*a11
**4*a22*a33*b11 - 4*a11**4*a33**2*b11 - 12*a11**3*a22**2*a33*b11 + 8*a11**3*a22*
a33**2*b11 + 4*a11**3*a33**3*b11 + 12*a11**2*a22**3*a33*b11 - 12*a11**2*a22*a33
**3*b11 - 4*a11*a22**4*a33*b11 - 8*a11*a22**3*a33**2*b11 + 12*a11*a22**2*a33**3*
b11 + 4*a22**4*a33**2*b11 - 4*a22**3*a33**3*b11) + v1**2*v2**2*( - a11**5*b11**2
+ 2*a11**4*a22*b11**2 + 3*a11**4*a33*b11**2 - 8*a11**3*a22*a33*b11**2 - 2*a11**
3*a33**2*b11**2 - 2*a11**2*a22**3*b11**2 + 6*a11**2*a22**2*a33*b11**2 + 6*a11**2
*a22*a33**2*b11**2 + a11*a22**4*b11**2 - 6*a11*a22**2*a33**2*b11**2 - a22**4*a33
*b11**2 + 2*a22**3*a33**2*b11**2) + v1**2*v3**2*( - a11**5*b11**2 + 4*a11**4*a22
*b11**2 + a11**4*a33*b11**2 - 5*a11**3*a22**2*b11**2 - 6*a11**3*a22*a33*b11**2 +
a11**3*a33**2*b11**2 + 2*a11**2*a22**3*b11**2 + 9*a11**2*a22**2*a33*b11**2 -
a11**2*a33**3*b11**2 - 4*a11*a22**3*a33*b11**2 - 3*a11*a22**2*a33**2*b11**2 + 2*
a11*a22*a33**3*b11**2 + 2*a22**3*a33**2*b11**2 - a22**2*a33**3*b11**2) + v2**4*(
- a11**5*b11**2 + 2*a11**4*a22*b11**2 + 3*a11**4*a33*b11**2 - 8*a11**3*a22*a33*
b11**2 - 2*a11**3*a33**2*b11**2 - 2*a11**2*a22**3*b11**2 + 6*a11**2*a22**2*a33*
b11**2 + 6*a11**2*a22*a33**2*b11**2 + a11*a22**4*b11**2 - 6*a11*a22**2*a33**2*
b11**2 - a22**4*a33*b11**2 + 2*a22**3*a33**2*b11**2) + v2**2*v3**2*( - 2*a11**5*
b11**2 + 6*a11**4*a22*b11**2 + 4*a11**4*a33*b11**2 - 5*a11**3*a22**2*b11**2 - 14
*a11**3*a22*a33*b11**2 - a11**3*a33**2*b11**2 + 15*a11**2*a22**2*a33*b11**2 + 6*
a11**2*a22*a33**2*b11**2 - a11**2*a33**3*b11**2 + a11*a22**4*b11**2 - 4*a11*a22
**3*a33*b11**2 - 9*a11*a22**2*a33**2*b11**2 + 2*a11*a22*a33**3*b11**2 - a22**4*
a33*b11**2 + 4*a22**3*a33**2*b11**2 - a22**2*a33**3*b11**2) + v3**4*( - a11**5*
b11**2 + 4*a11**4*a22*b11**2 + a11**4*a33*b11**2 - 5*a11**3*a22**2*b11**2 - 6*
a11**3*a22*a33*b11**2 + a11**3*a33**2*b11**2 + 2*a11**2*a22**3*b11**2 + 9*a11**2
*a22**2*a33*b11**2 - a11**2*a33**3*b11**2 - 4*a11*a22**3*a33*b11**2 - 3*a11*a22
**2*a33**2*b11**2 + 2*a11*a22*a33**3*b11**2 + 2*a22**3*a33**2*b11**2 - a22**2*
a33**3*b11**2)$
FI=u1**2*v2**2*( - 4*a11**5*a33**2*b11 + 12*a11**4*a22*a33**2*b11 + 8*a11**4*a33
**3*b11 - 12*a11**3*a22**2*a33**2*b11 - 24*a11**3*a22*a33**3*b11 - 4*a11**3*a33
**4*b11 + 4*a11**2*a22**3*a33**2*b11 + 24*a11**2*a22**2*a33**3*b11 + 12*a11**2*
a22*a33**4*b11 - 8*a11*a22**3*a33**3*b11 - 12*a11*a22**2*a33**4*b11 + 4*a22**3*
a33**4*b11) + u1**2*v3**2*( - 4*a11**5*a33**2*b11 + 12*a11**4*a22*a33**2*b11 + 8
*a11**4*a33**3*b11 - 12*a11**3*a22**2*a33**2*b11 - 24*a11**3*a22*a33**3*b11 - 4*
a11**3*a33**4*b11 + 4*a11**2*a22**3*a33**2*b11 + 24*a11**2*a22**2*a33**3*b11 +
12*a11**2*a22*a33**4*b11 - 8*a11*a22**3*a33**3*b11 - 12*a11*a22**2*a33**4*b11 +
4*a22**3*a33**4*b11) + u1*u2**2*v1*(4*a11**5*a33**3 - 16*a11**4*a22*a33**3 - 4*
a11**4*a33**4 + 24*a11**3*a22**2*a33**3 + 16*a11**3*a22*a33**4 - 16*a11**2*a22**
3*a33**3 - 24*a11**2*a22**2*a33**4 + 4*a11*a22**4*a33**3 + 16*a11*a22**3*a33**4
- 4*a22**4*a33**4) + u1*u2*v1*v2*(4*a11**5*a33**2*b11 - 12*a11**4*a22*a33**2*b11
- 8*a11**4*a33**3*b11 + 12*a11**3*a22**2*a33**2*b11 + 24*a11**3*a22*a33**3*b11
+ 4*a11**3*a33**4*b11 - 4*a11**2*a22**3*a33**2*b11 - 24*a11**2*a22**2*a33**3*b11
- 12*a11**2*a22*a33**4*b11 + 8*a11*a22**3*a33**3*b11 + 12*a11*a22**2*a33**4*b11
- 4*a22**3*a33**4*b11) + u1*u3**2*v1*(4*a11**5*a33**3 - 12*a11**4*a22*a33**3 -
8*a11**4*a33**4 + 12*a11**3*a22**2*a33**3 + 24*a11**3*a22*a33**4 + 4*a11**3*a33
**5 - 4*a11**2*a22**3*a33**3 - 24*a11**2*a22**2*a33**4 - 12*a11**2*a22*a33**5 +
8*a11*a22**3*a33**4 + 12*a11*a22**2*a33**5 - 4*a22**3*a33**5) + u1*u3*v1*v3*(4*
a11**5*a33**2*b11 - 12*a11**4*a22*a33**2*b11 - 8*a11**4*a33**3*b11 + 12*a11**3*
a22**2*a33**2*b11 + 24*a11**3*a22*a33**3*b11 + 4*a11**3*a33**4*b11 - 4*a11**2*
a22**3*a33**2*b11 - 24*a11**2*a22**2*a33**3*b11 - 12*a11**2*a22*a33**4*b11 + 8*
a11*a22**3*a33**3*b11 + 12*a11*a22**2*a33**4*b11 - 4*a22**3*a33**4*b11) + u1*v1
**3*(4*a11**4*a22*a33*b11**2 - 4*a11**4*a33**2*b11**2 - 8*a11**3*a22**2*a33*b11
**2 + 8*a11**3*a33**3*b11**2 + 4*a11**2*a22**3*a33*b11**2 + 12*a11**2*a22**2*a33
**2*b11**2 - 12*a11**2*a22*a33**3*b11**2 - 4*a11**2*a33**4*b11**2 - 8*a11*a22**3
*a33**2*b11**2 + 8*a11*a22*a33**4*b11**2 + 4*a22**3*a33**3*b11**2 - 4*a22**2*a33
**4*b11**2) + u1*v1*v2**2*(a11**5*a33*b11**2 + 2*a11**4*a22*a33*b11**2 - 7*a11**
4*a33**2*b11**2 - 7*a11**3*a22**2*a33*b11**2 + 6*a11**3*a22*a33**2*b11**2 + 11*
a11**3*a33**3*b11**2 + 4*a11**2*a22**3*a33*b11**2 + 9*a11**2*a22**2*a33**2*b11**
2 - 18*a11**2*a22*a33**3*b11**2 - 5*a11**2*a33**4*b11**2 - 8*a11*a22**3*a33**2*
b11**2 + 3*a11*a22**2*a33**3*b11**2 + 10*a11*a22*a33**4*b11**2 + 4*a22**3*a33**3
*b11**2 - 5*a22**2*a33**4*b11**2) + u1*v1*v3**2*(a11**5*a33*b11**2 + a11**4*a22*
a33*b11**2 - 6*a11**4*a33**2*b11**2 - 5*a11**3*a22**2*a33*b11**2 + 6*a11**3*a22*
a33**2*b11**2 + 9*a11**3*a33**3*b11**2 + 3*a11**2*a22**3*a33*b11**2 + 6*a11**2*
a22**2*a33**2*b11**2 - 15*a11**2*a22*a33**3*b11**2 - 4*a11**2*a33**4*b11**2 - 6*
a11*a22**3*a33**2*b11**2 + 3*a11*a22**2*a33**3*b11**2 + 8*a11*a22*a33**4*b11**2
+ 3*a22**3*a33**3*b11**2 - 4*a22**2*a33**4*b11**2) + u2**3*v2*(4*a11**5*a33**3 -
16*a11**4*a22*a33**3 - 4*a11**4*a33**4 + 24*a11**3*a22**2*a33**3 + 16*a11**3*
a22*a33**4 - 16*a11**2*a22**3*a33**3 - 24*a11**2*a22**2*a33**4 + 4*a11*a22**4*
a33**3 + 16*a11*a22**3*a33**4 - 4*a22**4*a33**4) + u2**2*u3*v3*(4*a11**5*a33**3
- 16*a11**4*a22*a33**3 - 4*a11**4*a33**4 + 24*a11**3*a22**2*a33**3 + 16*a11**3*
a22*a33**4 - 16*a11**2*a22**3*a33**3 - 24*a11**2*a22**2*a33**4 + 4*a11*a22**4*
a33**3 + 16*a11*a22**3*a33**4 - 4*a22**4*a33**4) + u2*u3**2*v2*(4*a11**5*a33**3
- 12*a11**4*a22*a33**3 - 8*a11**4*a33**4 + 12*a11**3*a22**2*a33**3 + 24*a11**3*
a22*a33**4 + 4*a11**3*a33**5 - 4*a11**2*a22**3*a33**3 - 24*a11**2*a22**2*a33**4
- 12*a11**2*a22*a33**5 + 8*a11*a22**3*a33**4 + 12*a11*a22**2*a33**5 - 4*a22**3*
a33**5) + u2*v2**3*(a11**5*a33*b11**2 - 2*a11**4*a22*a33*b11**2 - 3*a11**4*a33**
2*b11**2 + a11**3*a22**2*a33*b11**2 + 6*a11**3*a22*a33**2*b11**2 + 3*a11**3*a33
**3*b11**2 - 3*a11**2*a22**2*a33**2*b11**2 - 6*a11**2*a22*a33**3*b11**2 - a11**2
*a33**4*b11**2 + 3*a11*a22**2*a33**3*b11**2 + 2*a11*a22*a33**4*b11**2 - a22**2*
a33**4*b11**2) + u2*v2*v3**2*(a11**5*a33*b11**2 - 3*a11**4*a22*a33*b11**2 - 2*
a11**4*a33**2*b11**2 + 3*a11**3*a22**2*a33*b11**2 + 6*a11**3*a22*a33**2*b11**2 +
a11**3*a33**3*b11**2 - a11**2*a22**3*a33*b11**2 - 6*a11**2*a22**2*a33**2*b11**2
- 3*a11**2*a22*a33**3*b11**2 + 2*a11*a22**3*a33**2*b11**2 + 3*a11*a22**2*a33**3
*b11**2 - a22**3*a33**3*b11**2) + u3**3*v3*(4*a11**5*a33**3 - 12*a11**4*a22*a33
**3 - 8*a11**4*a33**4 + 12*a11**3*a22**2*a33**3 + 24*a11**3*a22*a33**4 + 4*a11**
3*a33**5 - 4*a11**2*a22**3*a33**3 - 24*a11**2*a22**2*a33**4 - 12*a11**2*a22*a33
**5 + 8*a11*a22**3*a33**4 + 12*a11*a22**2*a33**5 - 4*a22**3*a33**5) + u3**2*v1**
2*(4*a11**4*a22*a33**2*b11 - 4*a11**4*a33**3*b11 - 8*a11**3*a22**2*a33**2*b11 +
8*a11**3*a33**4*b11 + 4*a11**2*a22**3*a33**2*b11 + 12*a11**2*a22**2*a33**3*b11 -
12*a11**2*a22*a33**4*b11 - 4*a11**2*a33**5*b11 - 8*a11*a22**3*a33**3*b11 + 8*
a11*a22*a33**5*b11 + 4*a22**3*a33**4*b11 - 4*a22**2*a33**5*b11) + u3**2*v2**2*(4
*a11**4*a22*a33**2*b11 - 4*a11**4*a33**3*b11 - 8*a11**3*a22**2*a33**2*b11 + 8*
a11**3*a33**4*b11 + 4*a11**2*a22**3*a33**2*b11 + 12*a11**2*a22**2*a33**3*b11 -
12*a11**2*a22*a33**4*b11 - 4*a11**2*a33**5*b11 - 8*a11*a22**3*a33**3*b11 + 8*a11
*a22*a33**5*b11 + 4*a22**3*a33**4*b11 - 4*a22**2*a33**5*b11) + u3**2*v3**2*(4*
a11**4*a22*a33**2*b11 - 4*a11**4*a33**3*b11 - 8*a11**3*a22**2*a33**2*b11 + 8*a11
**3*a33**4*b11 + 4*a11**2*a22**3*a33**2*b11 + 12*a11**2*a22**2*a33**3*b11 - 12*
a11**2*a22*a33**4*b11 - 4*a11**2*a33**5*b11 - 8*a11*a22**3*a33**3*b11 + 8*a11*
a22*a33**5*b11 + 4*a22**3*a33**4*b11 - 4*a22**2*a33**5*b11) + u3*v1**2*v3*(4*a11
**4*a22*a33*b11**2 - 4*a11**4*a33**2*b11**2 - 8*a11**3*a22**2*a33*b11**2 + 8*a11
**3*a33**3*b11**2 + 4*a11**2*a22**3*a33*b11**2 + 12*a11**2*a22**2*a33**2*b11**2
- 12*a11**2*a22*a33**3*b11**2 - 4*a11**2*a33**4*b11**2 - 8*a11*a22**3*a33**2*b11
**2 + 8*a11*a22*a33**4*b11**2 + 4*a22**3*a33**3*b11**2 - 4*a22**2*a33**4*b11**2)
+ u3*v2**2*v3*(a11**5*a33*b11**2 + 2*a11**4*a22*a33*b11**2 - 7*a11**4*a33**2*
b11**2 - 7*a11**3*a22**2*a33*b11**2 + 6*a11**3*a22*a33**2*b11**2 + 11*a11**3*a33
**3*b11**2 + 4*a11**2*a22**3*a33*b11**2 + 9*a11**2*a22**2*a33**2*b11**2 - 18*a11
**2*a22*a33**3*b11**2 - 5*a11**2*a33**4*b11**2 - 8*a11*a22**3*a33**2*b11**2 + 3*
a11*a22**2*a33**3*b11**2 + 10*a11*a22*a33**4*b11**2 + 4*a22**3*a33**3*b11**2 - 5
*a22**2*a33**4*b11**2) + u3*v3**3*(a11**5*a33*b11**2 + a11**4*a22*a33*b11**2 - 6
*a11**4*a33**2*b11**2 - 5*a11**3*a22**2*a33*b11**2 + 6*a11**3*a22*a33**2*b11**2
+ 9*a11**3*a33**3*b11**2 + 3*a11**2*a22**3*a33*b11**2 + 6*a11**2*a22**2*a33**2*
b11**2 - 15*a11**2*a22*a33**3*b11**2 - 4*a11**2*a33**4*b11**2 - 6*a11*a22**3*a33
**2*b11**2 + 3*a11*a22**2*a33**3*b11**2 + 8*a11*a22*a33**4*b11**2 + 3*a22**3*a33
**3*b11**2 - 4*a22**2*a33**4*b11**2) + v1**2*v2**2*(a11**3*a22**2*b11**3 - 2*a11
**3*a22*a33*b11**3 + a11**3*a33**2*b11**3 - a11**2*a22**3*b11**3 + 3*a11**2*a22*
a33**2*b11**3 - 2*a11**2*a33**3*b11**3 + 2*a11*a22**3*a33*b11**3 - 3*a11*a22**2*
a33**2*b11**3 + a11*a33**4*b11**3 - a22**3*a33**2*b11**3 + 2*a22**2*a33**3*b11**
3 - a22*a33**4*b11**3) + v1**2*v3**2*(a11**4*a22*b11**3 - a11**4*a33*b11**3 -
a11**3*a22**2*b11**3 - 2*a11**3*a22*a33*b11**3 + 3*a11**3*a33**2*b11**3 + 3*a11
**2*a22**2*a33*b11**3 - 3*a11**2*a33**3*b11**3 - 3*a11*a22**2*a33**2*b11**3 + 2*
a11*a22*a33**3*b11**3 + a11*a33**4*b11**3 + a22**2*a33**3*b11**3 - a22*a33**4*
b11**3) + v2**4*(a11**3*a22**2*b11**3 - 2*a11**3*a22*a33*b11**3 + a11**3*a33**2*
b11**3 - a11**2*a22**3*b11**3 + 3*a11**2*a22*a33**2*b11**3 - 2*a11**2*a33**3*b11
**3 + 2*a11*a22**3*a33*b11**3 - 3*a11*a22**2*a33**2*b11**3 + a11*a33**4*b11**3 -
a22**3*a33**2*b11**3 + 2*a22**2*a33**3*b11**3 - a22*a33**4*b11**3) + v2**2*v3**
2*(a11**4*a22*b11**3 - a11**4*a33*b11**3 - 4*a11**3*a22*a33*b11**3 + 4*a11**3*
a33**2*b11**3 - a11**2*a22**3*b11**3 + 3*a11**2*a22**2*a33*b11**3 + 3*a11**2*a22
*a33**2*b11**3 - 5*a11**2*a33**3*b11**3 + 2*a11*a22**3*a33*b11**3 - 6*a11*a22**2
*a33**2*b11**3 + 2*a11*a22*a33**3*b11**3 + 2*a11*a33**4*b11**3 - a22**3*a33**2*
b11**3 + 3*a22**2*a33**3*b11**3 - 2*a22*a33**4*b11**3) + v3**4*(a11**4*a22*b11**
3 - a11**4*a33*b11**3 - a11**3*a22**2*b11**3 - 2*a11**3*a22*a33*b11**3 + 3*a11**
3*a33**2*b11**3 + 3*a11**2*a22**2*a33*b11**3 - 3*a11**2*a33**3*b11**3 - 3*a11*
a22**2*a33**2*b11**3 + 2*a11*a22*a33**3*b11**3 + a11*a33**4*b11**3 + a22**2*a33
**3*b11**3 - a22*a33**4*b11**3)$
FI=u1**3*v1*(32*a11**5*a33**3*b11 - 96*a11**4*a22*a33**3*b11 - 64*a11**4*a33**4*
b11 + 96*a11**3*a22**2*a33**3*b11 + 192*a11**3*a22*a33**4*b11 + 32*a11**3*a33**5
*b11 - 32*a11**2*a22**3*a33**3*b11 - 192*a11**2*a22**2*a33**4*b11 - 96*a11**2*
a22*a33**5*b11 + 64*a11*a22**3*a33**4*b11 + 96*a11*a22**2*a33**5*b11 - 32*a22**3
*a33**5*b11) + u1**2*u2**2*(32*a11**5*a33**4 - 128*a11**4*a22*a33**4 - 32*a11**4
*a33**5 + 192*a11**3*a22**2*a33**4 + 128*a11**3*a22*a33**5 - 128*a11**2*a22**3*
a33**4 - 192*a11**2*a22**2*a33**5 + 32*a11*a22**4*a33**4 + 128*a11*a22**3*a33**5
- 32*a22**4*a33**5) + u1**2*u3**2*(32*a11**5*a33**4 - 96*a11**4*a22*a33**4 - 64
*a11**4*a33**5 + 96*a11**3*a22**2*a33**4 + 192*a11**3*a22*a33**5 + 32*a11**3*a33
**6 - 32*a11**2*a22**3*a33**4 - 192*a11**2*a22**2*a33**5 - 96*a11**2*a22*a33**6
+ 64*a11*a22**3*a33**5 + 96*a11*a22**2*a33**6 - 32*a22**3*a33**6) + u1**2*v2**2*
(8*a11**5*a33**2*b11**2 - 40*a11**4*a33**3*b11**2 - 24*a11**3*a22**2*a33**2*b11
**2 + 48*a11**3*a22*a33**3*b11**2 + 56*a11**3*a33**4*b11**2 + 16*a11**2*a22**3*
a33**2*b11**2 + 24*a11**2*a22**2*a33**3*b11**2 - 96*a11**2*a22*a33**4*b11**2 -
24*a11**2*a33**5*b11**2 - 32*a11*a22**3*a33**3*b11**2 + 24*a11*a22**2*a33**4*b11
**2 + 48*a11*a22*a33**5*b11**2 + 16*a22**3*a33**4*b11**2 - 24*a22**2*a33**5*b11
**2) + u1**2*v3**2*(8*a11**5*a33**2*b11**2 - 8*a11**4*a22*a33**2*b11**2 - 32*a11
**4*a33**3*b11**2 - 8*a11**3*a22**2*a33**2*b11**2 + 48*a11**3*a22*a33**3*b11**2
+ 40*a11**3*a33**4*b11**2 + 8*a11**2*a22**3*a33**2*b11**2 - 72*a11**2*a22*a33**4
*b11**2 - 16*a11**2*a33**5*b11**2 - 16*a11*a22**3*a33**3*b11**2 + 24*a11*a22**2*
a33**4*b11**2 + 32*a11*a22*a33**5*b11**2 + 8*a22**3*a33**4*b11**2 - 16*a22**2*
a33**5*b11**2) + u1*u3**2*v1*( - 32*a11**4*a22*a33**3*b11 + 32*a11**4*a33**4*b11
+ 64*a11**3*a22**2*a33**3*b11 - 64*a11**3*a33**5*b11 - 32*a11**2*a22**3*a33**3*
b11 - 96*a11**2*a22**2*a33**4*b11 + 96*a11**2*a22*a33**5*b11 + 32*a11**2*a33**6*
b11 + 64*a11*a22**3*a33**4*b11 - 64*a11*a22*a33**6*b11 - 32*a22**3*a33**5*b11 +
32*a22**2*a33**6*b11) + u1*u3*v1*v3*( - 32*a11**4*a22*a33**2*b11**2 + 32*a11**4*
a33**3*b11**2 + 64*a11**3*a22**2*a33**2*b11**2 - 64*a11**3*a33**4*b11**2 - 32*
a11**2*a22**3*a33**2*b11**2 - 96*a11**2*a22**2*a33**3*b11**2 + 96*a11**2*a22*a33
**4*b11**2 + 32*a11**2*a33**5*b11**2 + 64*a11*a22**3*a33**3*b11**2 - 64*a11*a22*
a33**5*b11**2 - 32*a22**3*a33**4*b11**2 + 32*a22**2*a33**5*b11**2) + u1*v1**3*(
- 8*a11**4*a22*a33*b11**3 + 8*a11**4*a33**2*b11**3 + 32*a11**3*a22*a33**2*b11**3
- 32*a11**3*a33**3*b11**3 + 8*a11**2*a22**3*a33*b11**3 - 24*a11**2*a22**2*a33**
2*b11**3 - 24*a11**2*a22*a33**3*b11**3 + 40*a11**2*a33**4*b11**3 - 16*a11*a22**3
*a33**2*b11**3 + 48*a11*a22**2*a33**3*b11**3 - 16*a11*a22*a33**4*b11**3 - 16*a11
*a33**5*b11**3 + 8*a22**3*a33**3*b11**3 - 24*a22**2*a33**4*b11**3 + 16*a22*a33**
5*b11**3) + u1*v1*v2**2*( - 8*a11**4*a22*a33*b11**3 + 8*a11**4*a33**2*b11**3 - 8
*a11**3*a22**2*a33*b11**3 + 48*a11**3*a22*a33**2*b11**3 - 40*a11**3*a33**3*b11**
3 + 16*a11**2*a22**3*a33*b11**3 - 24*a11**2*a22**2*a33**2*b11**3 - 48*a11**2*a22
*a33**3*b11**3 + 56*a11**2*a33**4*b11**3 - 32*a11*a22**3*a33**2*b11**3 + 72*a11*
a22**2*a33**3*b11**3 - 16*a11*a22*a33**4*b11**3 - 24*a11*a33**5*b11**3 + 16*a22
**3*a33**3*b11**3 - 40*a22**2*a33**4*b11**3 + 24*a22*a33**5*b11**3) + u1*v1*v3**
2*( - 16*a11**4*a22*a33*b11**3 + 16*a11**4*a33**2*b11**3 + 8*a11**3*a22**2*a33*
b11**3 + 48*a11**3*a22*a33**2*b11**3 - 56*a11**3*a33**3*b11**3 + 8*a11**2*a22**3
*a33*b11**3 - 48*a11**2*a22**2*a33**2*b11**3 - 24*a11**2*a22*a33**3*b11**3 + 64*
a11**2*a33**4*b11**3 - 16*a11*a22**3*a33**2*b11**3 + 72*a11*a22**2*a33**3*b11**3
- 32*a11*a22*a33**4*b11**3 - 24*a11*a33**5*b11**3 + 8*a22**3*a33**3*b11**3 - 32
*a22**2*a33**4*b11**3 + 24*a22*a33**5*b11**3) + u2**4*(16*a11**4*a22*a33**4 - 16
*a11**4*a33**5 - 64*a11**3*a22**2*a33**4 + 64*a11**3*a22*a33**5 + 96*a11**2*a22
**3*a33**4 - 96*a11**2*a22**2*a33**5 - 64*a11*a22**4*a33**4 + 64*a11*a22**3*a33
**5 + 16*a22**5*a33**4 - 16*a22**4*a33**5) + u2**2*u3*v3*( - 32*a11**4*a22*a33**
3*b11 + 32*a11**4*a33**4*b11 + 96*a11**3*a22**2*a33**3*b11 - 64*a11**3*a22*a33**
4*b11 - 32*a11**3*a33**5*b11 - 96*a11**2*a22**3*a33**3*b11 + 96*a11**2*a22*a33**
5*b11 + 32*a11*a22**4*a33**3*b11 + 64*a11*a22**3*a33**4*b11 - 96*a11*a22**2*a33
**5*b11 - 32*a22**4*a33**4*b11 + 32*a22**3*a33**5*b11) + u2**2*v1**2*( - 8*a11**
4*a22*a33**2*b11**2 + 8*a11**4*a33**3*b11**2 + 24*a11**3*a22**2*a33**2*b11**2 -
16*a11**3*a22*a33**3*b11**2 - 8*a11**3*a33**4*b11**2 - 24*a11**2*a22**3*a33**2*
b11**2 + 24*a11**2*a22*a33**4*b11**2 + 8*a11*a22**4*a33**2*b11**2 + 16*a11*a22**
3*a33**3*b11**2 - 24*a11*a22**2*a33**4*b11**2 - 8*a22**4*a33**3*b11**2 + 8*a22**
3*a33**4*b11**2) + u2**2*v3**2*( - 8*a11**4*a22*a33**2*b11**2 + 8*a11**4*a33**3*
b11**2 + 16*a11**3*a22**2*a33**2*b11**2 - 16*a11**3*a33**4*b11**2 - 8*a11**2*a22
**3*a33**2*b11**2 - 24*a11**2*a22**2*a33**3*b11**2 + 24*a11**2*a22*a33**4*b11**2
+ 8*a11**2*a33**5*b11**2 + 16*a11*a22**3*a33**3*b11**2 - 16*a11*a22*a33**5*b11
**2 - 8*a22**3*a33**4*b11**2 + 8*a22**2*a33**5*b11**2) + u3**4*( - 16*a11**4*a22
*a33**4 + 16*a11**4*a33**5 + 32*a11**3*a22**2*a33**4 - 32*a11**3*a33**6 - 16*a11
**2*a22**3*a33**4 - 48*a11**2*a22**2*a33**5 + 48*a11**2*a22*a33**6 + 16*a11**2*
a33**7 + 32*a11*a22**3*a33**5 - 32*a11*a22*a33**7 - 16*a22**3*a33**6 + 16*a22**2
*a33**7) + u3**3*v3*( - 32*a11**4*a22*a33**3*b11 + 32*a11**4*a33**4*b11 + 64*a11
**3*a22**2*a33**3*b11 - 64*a11**3*a33**5*b11 - 32*a11**2*a22**3*a33**3*b11 - 96*
a11**2*a22**2*a33**4*b11 + 96*a11**2*a22*a33**5*b11 + 32*a11**2*a33**6*b11 + 64*
a11*a22**3*a33**4*b11 - 64*a11*a22*a33**6*b11 - 32*a22**3*a33**5*b11 + 32*a22**2
*a33**6*b11) + u3**2*v1**2*( - 8*a11**4*a22*a33**2*b11**2 + 8*a11**4*a33**3*b11
**2 + 32*a11**3*a22*a33**3*b11**2 - 32*a11**3*a33**4*b11**2 + 8*a11**2*a22**3*
a33**2*b11**2 - 24*a11**2*a22**2*a33**3*b11**2 - 24*a11**2*a22*a33**4*b11**2 +
40*a11**2*a33**5*b11**2 - 16*a11*a22**3*a33**3*b11**2 + 48*a11*a22**2*a33**4*b11
**2 - 16*a11*a22*a33**5*b11**2 - 16*a11*a33**6*b11**2 + 8*a22**3*a33**4*b11**2 -
24*a22**2*a33**5*b11**2 + 16*a22*a33**6*b11**2) + u3**2*v2**2*( - 8*a11**4*a22*
a33**2*b11**2 + 8*a11**4*a33**3*b11**2 - 8*a11**3*a22**2*a33**2*b11**2 + 48*a11
**3*a22*a33**3*b11**2 - 40*a11**3*a33**4*b11**2 + 16*a11**2*a22**3*a33**2*b11**2
- 24*a11**2*a22**2*a33**3*b11**2 - 48*a11**2*a22*a33**4*b11**2 + 56*a11**2*a33
**5*b11**2 - 32*a11*a22**3*a33**3*b11**2 + 72*a11*a22**2*a33**4*b11**2 - 16*a11*
a22*a33**5*b11**2 - 24*a11*a33**6*b11**2 + 16*a22**3*a33**4*b11**2 - 40*a22**2*
a33**5*b11**2 + 24*a22*a33**6*b11**2) + u3**2*v3**2*( - 16*a11**4*a22*a33**2*b11
**2 + 16*a11**4*a33**3*b11**2 + 8*a11**3*a22**2*a33**2*b11**2 + 48*a11**3*a22*
a33**3*b11**2 - 56*a11**3*a33**4*b11**2 + 8*a11**2*a22**3*a33**2*b11**2 - 48*a11
**2*a22**2*a33**3*b11**2 - 24*a11**2*a22*a33**4*b11**2 + 64*a11**2*a33**5*b11**2
- 16*a11*a22**3*a33**3*b11**2 + 72*a11*a22**2*a33**4*b11**2 - 32*a11*a22*a33**5
*b11**2 - 24*a11*a33**6*b11**2 + 8*a22**3*a33**4*b11**2 - 32*a22**2*a33**5*b11**
2 + 24*a22*a33**6*b11**2) + u3*v1**2*v3*( - 16*a11**3*a22**2*a33*b11**3 + 32*a11
**3*a22*a33**2*b11**3 - 16*a11**3*a33**3*b11**3 + 16*a11**2*a22**3*a33*b11**3 -
48*a11**2*a22*a33**3*b11**3 + 32*a11**2*a33**4*b11**3 - 32*a11*a22**3*a33**2*b11
**3 + 48*a11*a22**2*a33**3*b11**3 - 16*a11*a33**5*b11**3 + 16*a22**3*a33**3*b11
**3 - 32*a22**2*a33**4*b11**3 + 16*a22*a33**5*b11**3) + u3*v2**2*v3*( - 8*a11**4
*a22*a33*b11**3 + 8*a11**4*a33**2*b11**3 - 8*a11**3*a22**2*a33*b11**3 + 48*a11**
3*a22*a33**2*b11**3 - 40*a11**3*a33**3*b11**3 + 16*a11**2*a22**3*a33*b11**3 - 24
*a11**2*a22**2*a33**2*b11**3 - 48*a11**2*a22*a33**3*b11**3 + 56*a11**2*a33**4*
b11**3 - 32*a11*a22**3*a33**2*b11**3 + 72*a11*a22**2*a33**3*b11**3 - 16*a11*a22*
a33**4*b11**3 - 24*a11*a33**5*b11**3 + 16*a22**3*a33**3*b11**3 - 40*a22**2*a33**
4*b11**3 + 24*a22*a33**5*b11**3) + u3*v3**3*( - 8*a11**4*a22*a33*b11**3 + 8*a11
**4*a33**2*b11**3 + 32*a11**3*a22*a33**2*b11**3 - 32*a11**3*a33**3*b11**3 + 8*
a11**2*a22**3*a33*b11**3 - 24*a11**2*a22**2*a33**2*b11**3 - 24*a11**2*a22*a33**3
*b11**3 + 40*a11**2*a33**4*b11**3 - 16*a11*a22**3*a33**2*b11**3 + 48*a11*a22**2*
a33**3*b11**3 - 16*a11*a22*a33**4*b11**3 - 16*a11*a33**5*b11**3 + 8*a22**3*a33**
3*b11**3 - 24*a22**2*a33**4*b11**3 + 16*a22*a33**5*b11**3) + v1**2*v2**2*( - 2*
a11**4*a22*b11**4 + 2*a11**4*a33*b11**4 + 2*a11**3*a22**2*b11**4 + 4*a11**3*a22*
a33*b11**4 - 6*a11**3*a33**2*b11**4 - 4*a11**2*a22**3*b11**4 + 6*a11**2*a22**2*
a33*b11**4 - 12*a11**2*a22*a33**2*b11**4 + 10*a11**2*a33**3*b11**4 + 8*a11*a22**
3*a33*b11**4 - 18*a11*a22**2*a33**2*b11**4 + 20*a11*a22*a33**3*b11**4 - 10*a11*
a33**4*b11**4 - 4*a22**3*a33**2*b11**4 + 10*a22**2*a33**3*b11**4 - 10*a22*a33**4
*b11**4 + 4*a33**5*b11**4) + v1**2*v3**2*( - 2*a11**4*a22*b11**4 + 2*a11**4*a33*
b11**4 + 8*a11**3*a22*a33*b11**4 - 8*a11**3*a33**2*b11**4 - 2*a11**2*a22**3*b11
**4 + 6*a11**2*a22**2*a33*b11**4 - 18*a11**2*a22*a33**2*b11**4 + 14*a11**2*a33**
3*b11**4 + 4*a11*a22**3*a33*b11**4 - 12*a11*a22**2*a33**2*b11**4 + 20*a11*a22*
a33**3*b11**4 - 12*a11*a33**4*b11**4 - 2*a22**3*a33**2*b11**4 + 6*a22**2*a33**3*
b11**4 - 8*a22*a33**4*b11**4 + 4*a33**5*b11**4) + v2**4*( - a11**4*a22*b11**4 +
a11**4*a33*b11**4 + 4*a11**3*a22*a33*b11**4 - 4*a11**3*a33**2*b11**4 - 4*a11**2*
a22**3*b11**4 + 12*a11**2*a22**2*a33*b11**4 - 18*a11**2*a22*a33**2*b11**4 + 10*
a11**2*a33**3*b11**4 + 8*a11*a22**3*a33*b11**4 - 24*a11*a22**2*a33**2*b11**4 +
28*a11*a22*a33**3*b11**4 - 12*a11*a33**4*b11**4 - 4*a22**3*a33**2*b11**4 + 12*
a22**2*a33**3*b11**4 - 13*a22*a33**4*b11**4 + 5*a33**5*b11**4) + v2**2*v3**2*( -
4*a11**4*a22*b11**4 + 4*a11**4*a33*b11**4 - 2*a11**3*a22**2*b11**4 + 20*a11**3*
a22*a33*b11**4 - 18*a11**3*a33**2*b11**4 - 4*a11**2*a22**3*b11**4 + 18*a11**2*
a22**2*a33*b11**4 - 48*a11**2*a22*a33**2*b11**4 + 34*a11**2*a33**3*b11**4 + 8*
a11*a22**3*a33*b11**4 - 30*a11*a22**2*a33**2*b11**4 + 52*a11*a22*a33**3*b11**4 -
30*a11*a33**4*b11**4 - 4*a22**3*a33**2*b11**4 + 14*a22**2*a33**3*b11**4 - 20*
a22*a33**4*b11**4 + 10*a33**5*b11**4) + v3**4*( - 3*a11**4*a22*b11**4 + 3*a11**4
*a33*b11**4 + 12*a11**3*a22*a33*b11**4 - 12*a11**3*a33**2*b11**4 - a11**2*a22**3
*b11**4 + 3*a11**2*a22**2*a33*b11**4 - 21*a11**2*a22*a33**2*b11**4 + 19*a11**2*
a33**3*b11**4 + 2*a11*a22**3*a33*b11**4 - 6*a11*a22**2*a33**2*b11**4 + 18*a11*
a22*a33**3*b11**4 - 14*a11*a33**4*b11**4 - a22**3*a33**2*b11**4 + 3*a22**2*a33**
3*b11**4 - 6*a22*a33**4*b11**4 + 4*a33**5*b11**4)$
FI=u1*v1*(4*a11**5*a33*b11 - 16*a11**4*a22*a33*b11 - 4*a11**4*a33**2*b11 + 24*
a11**3*a22**2*a33*b11 + 16*a11**3*a22*a33**2*b11 - 16*a11**2*a22**3*a33*b11 - 24
*a11**2*a22**2*a33**2*b11 + 4*a11*a22**4*a33*b11 + 16*a11*a22**3*a33**2*b11 - 4*
a22**4*a33**2*b11) + u2**2*(4*a11**5*a33**2 - 20*a11**4*a22*a33**2 + 40*a11**3*
a22**2*a33**2 - 40*a11**2*a22**3*a33**2 + 20*a11*a22**4*a33**2 - 4*a22**5*a33**2
) + u3**2*(4*a11**5*a33**2 - 16*a11**4*a22*a33**2 - 4*a11**4*a33**3 + 24*a11**3*
a22**2*a33**2 + 16*a11**3*a22*a33**3 - 16*a11**2*a22**3*a33**2 - 24*a11**2*a22**
2*a33**3 + 4*a11*a22**4*a33**2 + 16*a11*a22**3*a33**3 - 4*a22**4*a33**3) + v2**2
*(a11**5*b11**2 - 3*a11**4*a22*b11**2 - 2*a11**4*a33*b11**2 + 3*a11**3*a22**2*
b11**2 + 6*a11**3*a22*a33*b11**2 + a11**3*a33**2*b11**2 - a11**2*a22**3*b11**2 -
6*a11**2*a22**2*a33*b11**2 - 3*a11**2*a22*a33**2*b11**2 + 2*a11*a22**3*a33*b11
**2 + 3*a11*a22**2*a33**2*b11**2 - a22**3*a33**2*b11**2) + v3**2*(a11**5*b11**2
- 4*a11**4*a22*b11**2 - a11**4*a33*b11**2 + 6*a11**3*a22**2*b11**2 + 4*a11**3*
a22*a33*b11**2 - 4*a11**2*a22**3*b11**2 - 6*a11**2*a22**2*a33*b11**2 + a11*a22**
4*b11**2 + 4*a11*a22**3*a33*b11**2 - a22**4*a33*b11**2)$