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