Solution 10 to problem e3c2new
Expressions |
Parameters |
Relevance |
Back to problem e3c2new
Expressions
The solution is given through the following expressions:
b22=0
b31=0
b32=0
c12=0
c13=0
c22=0
c23=0
1 2
- ---*b33
4
c33=-------------
a22
1 1
- ---*a22*b33*n1 - ---*a33*b33*n1
2 2
m1=------------------------------------
2
a22 - a22*a33
1 1
- ---*a22*b33*n2 - ---*a33*b33*n2
2 2
m2=------------------------------------
2
a22 - a22*a33
- b33*n3
m3=-----------
a22 - a33
- b33*n3*q1
r6=----------------
2
a22 - a22*a33
1
- ---*b33*n2*q1
2
r5=------------------
2
a22 - a22*a33
1
- ---*b33*n1*q1
2
r4=------------------
2
a22 - a22*a33
n3*q1
r3=-----------
a22 - a33
n2*q1
r2=-----------
a22 - a33
n1*q1
r1=-----------
a22 - a33
q20=0
q19=0
q17=0
q16=0
1 2
---*b33 *q1
4
q15=-------------
2
a22
q14=0
q13=0
- b33*q1
q12=-----------
a22
q11=0
1 2
---*b33 *q1
4
q10=-------------
2
a22
q9=0
q8=0
- b33*q1
q7=-----------
a22
q5=0
q4=0
q3=q1
q2=0
- b33*k26*n3
p56=----------------
2
a22 - a22*a33
1
- ---*b33*k26*n2
2
p55=-------------------
2
a22 - a22*a33
1
- ---*b33*k26*n1
2
p54=-------------------
2
a22 - a22*a33
k26*n3
p53=-----------
a22 - a33
k26*n2
p52=-----------
a22 - a33
k26*n1
p51=-----------
a22 - a33
- b33*k26*n3
p50=----------------
2
a22 - a22*a33
p49=0
1
- ---*b33*k23*n2
2
p48=-------------------
2
a22 - a22*a33
p47=0
p46=0
- b33*k26*n3
p45=----------------
2
a22 - a22*a33
1
- ---*b33*k23*n1
2
p44=-------------------
2
a22 - a22*a33
p43=0
p42=0
- k23*n3
p41=-----------
a22 - a33
- k23*n2
p40=-----------
a22 - a33
- k23*n1
p39=-----------
a22 - a33
p38=0
p37=0
p36=0
1
- ---*b33*k26*n2
2
p35=-------------------
2
a22 - a22*a33
1
- ---*b33*k26*n1
2
p34=-------------------
2
a22 - a22*a33
k26*n3
p33=-----------
a22 - a33
4 3 1 3 2 2
p32=(a22 *k26*n2 - 3*a22 *a33*k26*n2 - ---*a22 *b33*k23*n2 + 3*a22 *a33 *k26*n2
2
3 2 3 3 2
+ ---*a22 *a33*b33*k23*n2 - a22*a33 *k26*n2 - ---*a22*a33 *b33*k23*n2
2 2
1 3 5 4 3 2 2 3
+ ---*a33 *b33*k23*n2)/(a22 - 4*a22 *a33 + 6*a22 *a33 - 4*a22 *a33
2
4
+ a22*a33 )
k26*n1
p31=-----------
a22 - a33
1
- ---*b33*k26*n2
2
p30=-------------------
2
a22 - a22*a33
p29=0
1
- ---*b33*k23*n1
2
p28=-------------------
2
a22 - a22*a33
1
- ---*b33*k23*n2
2
p27=-------------------
2
a22 - a22*a33
p26=0
- k23*n3
p25=-----------
a22 - a33
p24=0
- k23*n2
p23=-----------
a22 - a33
- k23*n1
p22=-----------
a22 - a33
p21=0
1
- ---*b33*k26*n1
2
p20=-------------------
2
a22 - a22*a33
k26*n3
p19=-----------
a22 - a33
k26*n2
p18=-----------
a22 - a33
1
a22*k26*n1 - ---*b33*k23*n1
2
p17=-----------------------------
2
a22 - a22*a33
p16=0
p15=0
- k23*n3
p14=-----------
a22 - a33
p13=0
- k23*n2
p12=-----------
a22 - a33
- k23*n1
p11=-----------
a22 - a33
p10=0
p9=0
p8=0
p7=0
p6=0
p5=0
p4=0
p3=0
p2=0
p1=0
k125=0
k124=0
k122=0
k121=0
1 2
---*b33 *k26
4
k120=--------------
2
a22
k119=0
k118=0
- b33*k26
k117=------------
a22
k116=0
1 2
---*b33 *k26
4
k115=--------------
2
a22
k114=0
k113=0
- b33*k26
k112=------------
a22
k110=0
k109=0
k108=k26
k107=0
k106=k26
k105=0
k104=0
k103=0
k102=0
k101=0
k100=0
k99=0
k98=0
k97=0
k96=0
k95=0
k94=0
k93=0
k92=0
k91=0
k90=0
k89=0
k88=0
k87=0
k86=0
k85=0
k84=0
k83=0
k82=0
k81=0
k80=k23
k79=0
k78=0
k77=0
k76=0
k75=0
k74=0
k73=0
k72=0
k71=0
1 2
---*b33 *k26
4
k70=--------------
2
a22
k69=0
k68=0
- b33*k26
k67=------------
a22
k66=0
1 2
---*b33 *k26
2
k65=--------------
2
a22
k64=0
k63=0
- b33*k26
k62=------------
a22
k61=0
k60=0
k59=0
k58=k26
k57=0
k56=k26
k55=0
k54=0
- b33*k26
k53=------------
a22
k52=0
k51=0
k50=0
k49=0
k48=0
k47=0
k46=0
k45=0
k44=k23
k43=0
k42=0
k41=0
k40=0
k38=0
k37=0
k36=0
1 2
---*b33 *k26
4
k35=--------------
2
a22
k34=0
k33=0
- b33*k26
k32=------------
a22
k30=0
k29=0
k28=k26
k27=0
k25=0
k24=0
k22=0
k21=0
k20=0
k19=0
k18=0
k17=0
k16=0
k14=0
k13=0
k12=0
k11=0
k10=0
k9=0
k8=0
k7=0
k6=0
k5=0
k4=0
k3=0
k2=0
k1=0
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:
n3,q1,k26,k23,n1,b33,n2,a33,a22
Relevance for the application:
The following expression INT is a first
integral for the Hamiltonian HAM:
3 2 3 2 2 2 2 2 2 2
HAM=(4*a22 *u1 + 4*a22 *u2 - 4*a22 *a33*u1 - 4*a22 *a33*u2 + 4*a22 *a33*u3
2 2 2 2
+ 4*a22 *b33*u3*v3 + 4*a22 *n1*u1 + 4*a22 *n2*u2 + 4*a22 *n3*u3
2 2
- 4*a22*a33 *u3 - 4*a22*a33*b33*u3*v3 - 4*a22*a33*n1*u1 - 4*a22*a33*n2*u2
2 2
- 4*a22*a33*n3*u3 - a22*b33 *v3 - 2*a22*b33*n1*v1 - 2*a22*b33*n2*v2
2 2
- 4*a22*b33*n3*v3 + a33*b33 *v3 - 2*a33*b33*n1*v1 - 2*a33*b33*n2*v2)/(4
*a22*(a22 - a33))
3 2 3 2 3 3
INT=(4*a22 *k23*u1*u3 *v1 + 4*a22 *k23*u2*u3 *v2 + 4*a22 *k23*u3 *v3
3 2 2 3 2 2 3 2 2
+ 4*a22 *k26*u1 *v1 + 4*a22 *k26*u1 *v2 + 4*a22 *k26*u1 *v3
3 2 2 3 2 2 3 2 2
+ 4*a22 *k26*u2 *v1 + 4*a22 *k26*u2 *v2 + 4*a22 *k26*u2 *v3
3 2 3 2 2 2
+ 4*a22 *q1*u1 + 4*a22 *q1*u2 - 4*a22 *a33*k23*u1*u3 *v1
2 2 2 3
- 4*a22 *a33*k23*u2*u3 *v2 - 4*a22 *a33*k23*u3 *v3
2 2 2 2 2 2 2 2 2
- 4*a22 *a33*k26*u1 *v1 - 4*a22 *a33*k26*u1 *v2 - 4*a22 *a33*k26*u1 *v3
2 2 2 2 2 2 2 2 2
- 4*a22 *a33*k26*u2 *v1 - 4*a22 *a33*k26*u2 *v2 - 4*a22 *a33*k26*u2 *v3
2 2 2 2 2 3
- 4*a22 *a33*q1*u1 - 4*a22 *a33*q1*u2 - 4*a22 *b33*k26*u1*v1
2 2 2 2
- 4*a22 *b33*k26*u1*v1*v2 - 4*a22 *b33*k26*u1*v1*v3
2 2 2 3
- 4*a22 *b33*k26*u2*v1 *v2 - 4*a22 *b33*k26*u2*v2
2 2 2 2
- 4*a22 *b33*k26*u2*v2*v3 - 4*a22 *b33*q1*u1*v1 - 4*a22 *b33*q1*u2*v2
2 2 2 2
- 4*a22 *k23*n1*u1 *v1 - 4*a22 *k23*n1*u1*u2*v2 - 4*a22 *k23*n1*u1*u3*v3
2 2 2 2
- 4*a22 *k23*n2*u1*u2*v1 - 4*a22 *k23*n2*u2 *v2 - 4*a22 *k23*n2*u2*u3*v3
2 2 2 2
- 4*a22 *k23*n3*u1*u3*v1 - 4*a22 *k23*n3*u2*u3*v2 - 4*a22 *k23*n3*u3 *v3
2 2 2 2 2 2
+ 4*a22 *k26*n1*u1*v1 + 4*a22 *k26*n1*u1*v2 + 4*a22 *k26*n1*u1*v3
2 2 2 2 2 2
+ 4*a22 *k26*n2*u2*v1 + 4*a22 *k26*n2*u2*v2 + 4*a22 *k26*n2*u2*v3
2 2 2 2 2 2
+ 4*a22 *k26*n3*u3*v1 + 4*a22 *k26*n3*u3*v2 + 4*a22 *k26*n3*u3*v3
2 2 2
+ 4*a22 *n1*q1*u1 + 4*a22 *n2*q1*u2 + 4*a22 *n3*q1*u3
3 2
+ 4*a22*a33*b33*k26*u1*v1 + 4*a22*a33*b33*k26*u1*v1*v2
2 2
+ 4*a22*a33*b33*k26*u1*v1*v3 + 4*a22*a33*b33*k26*u2*v1 *v2
3 2
+ 4*a22*a33*b33*k26*u2*v2 + 4*a22*a33*b33*k26*u2*v2*v3
2 4
+ 4*a22*a33*b33*q1*u1*v1 + 4*a22*a33*b33*q1*u2*v2 + a22*b33 *k26*v1
2 2 2 2 2 2 2 4
+ 2*a22*b33 *k26*v1 *v2 + a22*b33 *k26*v1 *v3 + a22*b33 *k26*v2
2 2 2 2 2 2 2
+ a22*b33 *k26*v2 *v3 + a22*b33 *q1*v1 + a22*b33 *q1*v2
2
- 2*a22*b33*k23*n1*u1*v1 - 2*a22*b33*k23*n1*u2*v1*v2
- 2*a22*b33*k23*n1*u3*v1*v3 - 2*a22*b33*k23*n2*u1*v1*v2
2
- 2*a22*b33*k23*n2*u2*v2 - 2*a22*b33*k23*n2*u3*v2*v3
3 2 2
- 2*a22*b33*k26*n1*v1 - 2*a22*b33*k26*n1*v1*v2 - 2*a22*b33*k26*n1*v1*v3
2 3 2
- 2*a22*b33*k26*n2*v1 *v2 - 2*a22*b33*k26*n2*v2 - 2*a22*b33*k26*n2*v2*v3
2 2 3
- 4*a22*b33*k26*n3*v1 *v3 - 4*a22*b33*k26*n3*v2 *v3 - 4*a22*b33*k26*n3*v3
- 2*a22*b33*n1*q1*v1 - 2*a22*b33*n2*q1*v2 - 4*a22*b33*n3*q1*v3
2 4 2 2 2 2 2 2
- a33*b33 *k26*v1 - 2*a33*b33 *k26*v1 *v2 - a33*b33 *k26*v1 *v3
2 4 2 2 2 2 2
- a33*b33 *k26*v2 - a33*b33 *k26*v2 *v3 - a33*b33 *q1*v1
2 2 2
- a33*b33 *q1*v2 )/(4*a22 *(a22 - a33))
And again in machine readable form:
HAM=(4*a22**3*u1**2 + 4*a22**3*u2**2 - 4*a22**2*a33*u1**2 - 4*a22**2*a33*u2**2 +
4*a22**2*a33*u3**2 + 4*a22**2*b33*u3*v3 + 4*a22**2*n1*u1 + 4*a22**2*n2*u2 + 4*
a22**2*n3*u3 - 4*a22*a33**2*u3**2 - 4*a22*a33*b33*u3*v3 - 4*a22*a33*n1*u1 - 4*
a22*a33*n2*u2 - 4*a22*a33*n3*u3 - a22*b33**2*v3**2 - 2*a22*b33*n1*v1 - 2*a22*b33
*n2*v2 - 4*a22*b33*n3*v3 + a33*b33**2*v3**2 - 2*a33*b33*n1*v1 - 2*a33*b33*n2*v2)
/(4*a22*(a22 - a33))$
INT=(4*a22**3*k23*u1*u3**2*v1 + 4*a22**3*k23*u2*u3**2*v2 + 4*a22**3*k23*u3**3*v3
+ 4*a22**3*k26*u1**2*v1**2 + 4*a22**3*k26*u1**2*v2**2 + 4*a22**3*k26*u1**2*v3**
2 + 4*a22**3*k26*u2**2*v1**2 + 4*a22**3*k26*u2**2*v2**2 + 4*a22**3*k26*u2**2*v3
**2 + 4*a22**3*q1*u1**2 + 4*a22**3*q1*u2**2 - 4*a22**2*a33*k23*u1*u3**2*v1 - 4*
a22**2*a33*k23*u2*u3**2*v2 - 4*a22**2*a33*k23*u3**3*v3 - 4*a22**2*a33*k26*u1**2*
v1**2 - 4*a22**2*a33*k26*u1**2*v2**2 - 4*a22**2*a33*k26*u1**2*v3**2 - 4*a22**2*
a33*k26*u2**2*v1**2 - 4*a22**2*a33*k26*u2**2*v2**2 - 4*a22**2*a33*k26*u2**2*v3**
2 - 4*a22**2*a33*q1*u1**2 - 4*a22**2*a33*q1*u2**2 - 4*a22**2*b33*k26*u1*v1**3 -
4*a22**2*b33*k26*u1*v1*v2**2 - 4*a22**2*b33*k26*u1*v1*v3**2 - 4*a22**2*b33*k26*
u2*v1**2*v2 - 4*a22**2*b33*k26*u2*v2**3 - 4*a22**2*b33*k26*u2*v2*v3**2 - 4*a22**
2*b33*q1*u1*v1 - 4*a22**2*b33*q1*u2*v2 - 4*a22**2*k23*n1*u1**2*v1 - 4*a22**2*k23
*n1*u1*u2*v2 - 4*a22**2*k23*n1*u1*u3*v3 - 4*a22**2*k23*n2*u1*u2*v1 - 4*a22**2*
k23*n2*u2**2*v2 - 4*a22**2*k23*n2*u2*u3*v3 - 4*a22**2*k23*n3*u1*u3*v1 - 4*a22**2
*k23*n3*u2*u3*v2 - 4*a22**2*k23*n3*u3**2*v3 + 4*a22**2*k26*n1*u1*v1**2 + 4*a22**
2*k26*n1*u1*v2**2 + 4*a22**2*k26*n1*u1*v3**2 + 4*a22**2*k26*n2*u2*v1**2 + 4*a22
**2*k26*n2*u2*v2**2 + 4*a22**2*k26*n2*u2*v3**2 + 4*a22**2*k26*n3*u3*v1**2 + 4*
a22**2*k26*n3*u3*v2**2 + 4*a22**2*k26*n3*u3*v3**2 + 4*a22**2*n1*q1*u1 + 4*a22**2
*n2*q1*u2 + 4*a22**2*n3*q1*u3 + 4*a22*a33*b33*k26*u1*v1**3 + 4*a22*a33*b33*k26*
u1*v1*v2**2 + 4*a22*a33*b33*k26*u1*v1*v3**2 + 4*a22*a33*b33*k26*u2*v1**2*v2 + 4*
a22*a33*b33*k26*u2*v2**3 + 4*a22*a33*b33*k26*u2*v2*v3**2 + 4*a22*a33*b33*q1*u1*
v1 + 4*a22*a33*b33*q1*u2*v2 + a22*b33**2*k26*v1**4 + 2*a22*b33**2*k26*v1**2*v2**
2 + a22*b33**2*k26*v1**2*v3**2 + a22*b33**2*k26*v2**4 + a22*b33**2*k26*v2**2*v3
**2 + a22*b33**2*q1*v1**2 + a22*b33**2*q1*v2**2 - 2*a22*b33*k23*n1*u1*v1**2 - 2*
a22*b33*k23*n1*u2*v1*v2 - 2*a22*b33*k23*n1*u3*v1*v3 - 2*a22*b33*k23*n2*u1*v1*v2
- 2*a22*b33*k23*n2*u2*v2**2 - 2*a22*b33*k23*n2*u3*v2*v3 - 2*a22*b33*k26*n1*v1**3
- 2*a22*b33*k26*n1*v1*v2**2 - 2*a22*b33*k26*n1*v1*v3**2 - 2*a22*b33*k26*n2*v1**
2*v2 - 2*a22*b33*k26*n2*v2**3 - 2*a22*b33*k26*n2*v2*v3**2 - 4*a22*b33*k26*n3*v1
**2*v3 - 4*a22*b33*k26*n3*v2**2*v3 - 4*a22*b33*k26*n3*v3**3 - 2*a22*b33*n1*q1*v1
- 2*a22*b33*n2*q1*v2 - 4*a22*b33*n3*q1*v3 - a33*b33**2*k26*v1**4 - 2*a33*b33**2
*k26*v1**2*v2**2 - a33*b33**2*k26*v1**2*v3**2 - a33*b33**2*k26*v2**4 - a33*b33**
2*k26*v2**2*v3**2 - a33*b33**2*q1*v1**2 - a33*b33**2*q1*v2**2)/(4*a22**2*(a22 -
a33))$