Solution 49 to problem e3quant


Expressions | Parameters | Relevance | Back to problem e3quant

Expressions

The solution is given through the following expressions:

a22= - a33


b22=0


b31=0


b32=0


c12=0


c13=0


c22=0


c23=0


      1     2
     ---*b33
      4
c33=----------
       a33


n1=0


n3=0


m1=0


m2=0


m3=0


r6=0


r4=0


r3=0


         2                     1
    2*a33 *r5 - 2*b33*k1*n2 - ---*b33*k12*n2
                               2
r2=------------------------------------------
                    a33*b33


r1=0


q20=0


           4           3     2    2          1     2    2
q19=( - a33 *b33*r5 - ---*a33 *b33 *k1*n2 + ---*a33 *b33 *k12*n2
                       4                     4

         1      2      3    1      2       3      4
      + ----*b33 *k1*n2  + ----*b33 *k12*n2 )/(a33 *n2)
         16                 16


      1     2
     ---*b33 *k12*n2
      4
q18=-----------------
             3
          a33


q17=0


           2           3     2          1     2
      - a33 *b33*r5 - ---*b33 *k1*n2 + ---*b33 *k12*n2
                       4                4
q16=---------------------------------------------------
                             2
                          a33 *n2


q14=0


         1                1
      - ---*a33*k38*n2 - ---*b33*k1*n2
         2                2
q13=-----------------------------------
                      2
                   a33


q11=0


q10=

         4         3             2                2               1           3
  - 4*a33 *r5 - a33 *k38*n2 - a33 *b33*k1*n2 + a33 *b33*k12*n2 + ---*b33*k1*n2
                                                                  4
--------------------------------------------------------------------------------
                                       3
                                    a33 *n2


q9=0


q8=0


            4         2                2               1           3
q7=( - 4*a33 *r5 + a33 *b33*k1*n2 + a33 *b33*k12*n2 + ---*b33*k1*n2
                                                       4

        1            3      2
     - ---*b33*k12*n2 )/(a33 *b33*n2)
        4


     1
    ---*b33*k12*n2
     2
q6=----------------
            2
         a33


q5=0


            2
     - 4*a33 *r5 - a33*k38*n2 - b33*k1*n2 + b33*k12*n2
q4=----------------------------------------------------
                          a33*n2


    k1*n2
q3=-------
     a33


q2=0


            2
     - 4*a33 *r5 + 5*b33*k1*n2 + b33*k12*n2
q1=-----------------------------------------
                    b33*n2


p56=0


         1     2               1     3          1     3
      - ---*a33 *b33*k16*n2 - ---*b33 *k1*n2 + ---*b33 *k12*n2
         4                     8                8
p55=-----------------------------------------------------------
                                  4
                               a33


p54=0


         1     2               1     3          1      3
      - ---*a33 *b33*k16*n2 - ---*b33 *k1*n2 + ----*b33 *k12*n2
         4                     4                16
p53=------------------------------------------------------------
                                   4
                                a33


p52=0


p51=0


p50=0


p49=0


         1     2               1     3          1      3
      - ---*a33 *b33*k16*n2 - ---*b33 *k1*n2 + ----*b33 *k12*n2
         4                     4                16
p48=------------------------------------------------------------
                                   4
                                a33


p47=0


p46=0


         1                    1     2          1     2
      - ---*a33*b33*k38*n2 + ---*b33 *k1*n2 + ---*b33 *k12*n2
         4                    2                4
p45=----------------------------------------------------------
                                  3
                               a33


p44=0


p43=0


        2       1     2
     b33 *k1 - ---*b33 *k12
                4
p42=------------------------
                 2
              a33


p41=0


p40=0


         1
      - ---*b33*k12*n2
         4
p39=-------------------
              2
           a33


p38=0


p37=0


         1     2           1     2
      - ---*a33 *k16*n2 - ---*b33 *k1*n2
         2                 4
p36=-------------------------------------
                       3
                    a33


p35=0


         1     2           1                    1     2          1     2
      - ---*a33 *k16*n2 - ---*a33*b33*k38*n2 - ---*b33 *k1*n2 + ---*b33 *k12*n2
         2                 4                    2                8
p34=----------------------------------------------------------------------------
                                           3
                                        a33


      1     2
     ---*b33 *k12
      4
p33=--------------
            2
         a33


p32=0


         1     2           1     2          1     2
      - ---*a33 *k16*n2 - ---*b33 *k1*n2 + ---*b33 *k12*n2
         2                 2                8
p31=-------------------------------------------------------
                                3
                             a33


      1                            1
     ---*a33*k38*n2 + b33*k1*n2 - ---*b33*k12*n2
      2                            2
p30=---------------------------------------------
                           2
                        a33


p29=0


     2*a33*k38 + 2*b33*k1
p28=----------------------
             a33


         1
      - ---*k12*n2
         2
p27=---------------
          a33


p26=0


      1                1               1
     ---*a33*k38*n2 - ---*b33*k1*n2 - ---*b33*k12*n2
      2                2               4
p25=-------------------------------------------------
                             2
                          a33


p24=0


p23=0


                 1
      - k1*n2 + ---*k12*n2
                 2
p22=-----------------------
              a33


p21=0


         1     2
      - ---*b33 *k12
         4
p20=-----------------
             2
          a33


p19=0


p18=0


         1
      - ---*b33*k38*n2
         4
p17=-------------------
              2
           a33


p16=0


p15=0


      - 2*a33*k38 - 2*b33*k1
p14=-------------------------
               a33


p13=0


p12=0


p11=0


p10=0


     1
    ---*k38*n2
     2
p9=------------
       a33


p8= - 4*k1


p7=0


p6=0


        1               1
     - ---*b33*k1*n2 - ---*b33*k12*n2
        2               4
p5=-----------------------------------
                     2
                  a33


p4=0


p3=0


                1
     - k1*n2 + ---*k12*n2
                2
p2=-----------------------
             a33


p1=0


k125=0


       1     2    2        1     4       1      4
      ---*a33 *b33 *k16 + ---*b33 *k1 - ----*b33 *k12
       4                   8             16
k124=-------------------------------------------------
                              4
                           a33


k123=0


       1     2    2        3      4       1      4
      ---*a33 *b33 *k16 + ----*b33 *k1 - ----*b33 *k12
       4                   16             16
k122=--------------------------------------------------
                               4
                            a33


k121=0


k120=0


k119=0


k118=0


       1     2    2        1     4       1      4
      ---*a33 *b33 *k16 + ---*b33 *k1 - ----*b33 *k12
       4                   8             16
k117=-------------------------------------------------
                              4
                           a33


k116=0


       1     2    2        3     4       1     4
      ---*a33 *b33 *k16 + ---*b33 *k1 - ---*b33 *k12
       2                   8             8
k115=------------------------------------------------
                              4
                           a33


k114=0


k113=0


       1     2    2        3      4       1      4
      ---*a33 *b33 *k16 + ----*b33 *k1 - ----*b33 *k12
       4                   16             16
k112=--------------------------------------------------
                               4
                            a33


k110=0


          1     3
       - ---*b33 *k1
          2
k109=----------------
              3
           a33


k108=0


k107=0


k106=0


k105=0


          1     3
       - ---*b33 *k1
          2
k104=----------------
              3
           a33


k103=0


k102=0


k100=0


k99=0


k98=0


k97=0


     a33*k38 - b33*k12
k95=-------------------
            a33


k94=0


k93=0


k91=0


        2            1     3       1     3
     a33 *b33*k16 + ---*b33 *k1 - ---*b33 *k12
                     2             4
k90=-------------------------------------------
                          3
                       a33


k89=0


        2            1     3       1     3
     a33 *b33*k16 + ---*b33 *k1 - ---*b33 *k12
                     2             4
k88=-------------------------------------------
                          3
                       a33


k87=0


k86=0


k85=0


k84=0


        2            1     3       1     3
     a33 *b33*k16 + ---*b33 *k1 - ---*b33 *k12
                     2             4
k83=-------------------------------------------
                          3
                       a33


k82=0


k81=0


             2       1     2
      - 2*b33 *k1 + ---*b33 *k12
                     2
k80=-----------------------------
                   2
                a33


k79=0


k78=0


k77=0


k76=0


k75=0


k74=k38


k73=0


k72=0


        2        1     2       1     2
     a33 *k16 + ---*b33 *k1 - ---*b33 *k12
                 2             4
k71=---------------------------------------
                        2
                     a33


k70=0


k69=k16


k68=0


k67=0


        2          2       1     2
     a33 *k16 + b33 *k1 - ---*b33 *k12
                           4
k66=-----------------------------------
                      2
                   a33


      - 2*b33*k1
k65=-------------
         a33


k64=0


k63=0


k62=k12


k61=0


k59=0


k58=0


k57=k1


k56=0


k55=0


k54=0


k53=0


        2            1     3       1     3
     a33 *b33*k16 + ---*b33 *k1 - ---*b33 *k12
                     2             4
k52=-------------------------------------------
                          3
                       a33


k51=0


        2            1     3       1     3
     a33 *b33*k16 + ---*b33 *k1 - ---*b33 *k12
                     2             4
k50=-------------------------------------------
                          3
                       a33


k49=0


k48=0


        2            1     3       1     3
     a33 *b33*k16 + ---*b33 *k1 - ---*b33 *k12
                     2             4
k47=-------------------------------------------
                          3
                       a33


k46=0


k45=0


k44=0


             2       1     2
      - 2*b33 *k1 + ---*b33 *k12
                     2
k43=-----------------------------
                   2
                a33


k42=0


k41=0


k40=0


k39=0


k37=0


k36=0


k35=0


k34=0


k33=0


             2       1     2
      - 2*b33 *k1 + ---*b33 *k12
                     2
k32=-----------------------------
                   2
                a33


k31=0


k30=0


k29=0


k28=0


k27=0


k26=0


k25=0


k24=0


k23=0


k22=0


        2        1     2       1     2
     a33 *k16 + ---*b33 *k1 - ---*b33 *k12
                 2             4
k21=---------------------------------------
                        2
                     a33


k20=0


        2          2       1     2
     a33 *k16 + b33 *k1 - ---*b33 *k12
                           4
k19=-----------------------------------
                      2
                   a33


k18=0


k17=0


      - 2*b33*k1
k15=-------------
         a33


k14=0


k13=0


k11=0


k10=0


k9=0


k8=0


k7=2*k1


k6=0


k5=0


k4=0


k3=0


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:
 k38,k16,k12,k1,r5,b33,n2,a33

Relevance for the application:

The following expression INT is a first

integral for the Hamiltonian HAM:

             2   2        2   2        2   2
HAM=( - 4*a33 *u1  - 4*a33 *u2  + 4*a33 *u3  + 4*a33*b33*u3*v3 + 4*a33*n2*u2

           2   2
      + b33 *v3 )/(4*a33)

              6      2         6                  5
INT=( - 64*a33 *r5*u1  - 64*a33 *r5*u1*v1 - 64*a33 *b33*r5*u1*u3

              5          2         5                  4    2
      - 64*a33 *b33*r5*v1  + 32*a33 *n2*r5*u2 - 16*a33 *b33 *r5*u1*v3

              4    2                  4             4
      - 16*a33 *b33 *r5*v1*v3 + 16*a33 *b33*k1*n2*u1

              4             2               4             2
      + 32*a33 *b33*k1*n2*u1 *u2*u3 + 80*a33 *b33*k1*n2*u1

              4                   2         4                2
      + 16*a33 *b33*k1*n2*u1*u2*v2  - 64*a33 *b33*k1*n2*u1*u3

              4                         4              2
      + 16*a33 *b33*k1*n2*u1*v1 + 16*a33 *b33*k12*n2*u1

              4                    2         4
      + 16*a33 *b33*k12*n2*u1*v1*v2  + 16*a33 *b33*k12*n2*u1*v1

              4              2   2         4              3
      + 16*a33 *b33*k12*n2*u2 *u3  + 16*a33 *b33*k16*n2*u1 *v1

              4              3            4
      + 16*a33 *b33*k16*n2*u1 *v3 + 16*a33 *b33*k16*n2*u1*u2*u3*v1

              4                 3         4              3
      + 16*a33 *b33*k16*n2*u1*v2  + 16*a33 *b33*k16*n2*u2 *v1

              4                 3         4                 2
      + 16*a33 *b33*k16*n2*v1*v2  + 16*a33 *b33*k38*n2*u1*u2 *v2

              4                             4
      - 32*a33 *b33*k38*n2*u1*u3*v1 - 16*a33 *b33*k38*n2*u1*u3

              4              3            4
      + 16*a33 *b33*k38*n2*u2 *v3 + 16*a33 *b33*k38*n2*u2*u3*v2*v3

              4                             4              2
      + 32*a33 *b33*k38*n2*u2*v1*v2 - 16*a33 *b33*k38*n2*v1

              4                      3    2
      + 16*a33 *b33*n2*r5*v2 - 32*a33 *b33 *k1*n2*u1*u3*v1

              3    2                     3    2
      - 16*a33 *b33 *k1*n2*u1*u3 + 32*a33 *b33 *k1*n2*u2*v1*v2

              3    2         4         3    2         2   2
      - 32*a33 *b33 *k1*n2*u3  - 32*a33 *b33 *k1*n2*v1 *v2

              3    2         2         3    2
      - 16*a33 *b33 *k1*n2*v1  + 16*a33 *b33 *k12*n2*u1*u3

              3    2                            3    2          2
      - 16*a33 *b33 *k12*n2*u2*u3*v2*v3 + 16*a33 *b33 *k12*n2*v1

              3    2                            3    2             2
      + 16*a33 *b33 *k16*n2*u1*u2*v1*v2 + 16*a33 *b33 *k16*n2*u1*v1 *v2

              3    2          2               3    2
      + 16*a33 *b33 *k16*n2*u2 *v1*v3 + 16*a33 *b33 *k16*n2*u2*u3*v1*v2

              3    2             2            3    2          3
      + 16*a33 *b33 *k16*n2*u2*v1 *v3 + 16*a33 *b33 *k16*n2*v1 *v3

              3          2   2            3          2
      - 16*a33 *b33*k1*n2 *u1 *u2 - 16*a33 *b33*k1*n2 *u1*u2*v2

              3          2   2         3          2
      + 16*a33 *b33*k1*n2 *u2  - 32*a33 *b33*k1*n2 *u2

             3           2   2           3           2
      + 8*a33 *b33*k12*n2 *u1 *u2 + 8*a33 *b33*k12*n2 *u1*u2*v2

             3           2                 3           2
      - 8*a33 *b33*k12*n2 *u1*v1*v2 - 8*a33 *b33*k12*n2 *u2

             3           2   2           3           2      2
      - 8*a33 *b33*k16*n2 *u1 *v3 - 8*a33 *b33*k16*n2 *u1*v2

             3           2      2        3           2      2
      - 8*a33 *b33*k16*n2 *v1*v2  + 8*a33 *b33*k38*n2 *u2*u3

             3           2                 3           2
      + 8*a33 *b33*k38*n2 *u2*u3*v2 - 8*a33 *b33*k38*n2 *u3*v2

             3           2   2           2    3         3
      + 8*a33 *b33*k38*n2 *v1 *v2 + 8*a33 *b33 *k1*n2*u1 *v3

             2    3                           2    3            2
      + 8*a33 *b33 *k1*n2*u1*u2*u3*v1 - 32*a33 *b33 *k1*n2*u1*u3 *v2

              2    3            3         2    3
      - 32*a33 *b33 *k1*n2*u1*v1  + 16*a33 *b33 *k1*n2*u1*v1*v3

              2    3            3         2    3
      + 16*a33 *b33 *k1*n2*u1*v2  - 12*a33 *b33 *k1*n2*u1*v3

              2    3         3            2    3         3
      + 16*a33 *b33 *k1*n2*u2 *v1 - 32*a33 *b33 *k1*n2*u3 *v3

              2    3                    2    3          3
      - 12*a33 *b33 *k1*n2*v1*v3 - 4*a33 *b33 *k12*n2*u1 *v3

             2    3                           2    3             2
      - 4*a33 *b33 *k12*n2*u1*u2*u3*v1 + 8*a33 *b33 *k12*n2*u1*u3 *v2

             2    3             3        2    3
      + 8*a33 *b33 *k12*n2*u1*v1  - 4*a33 *b33 *k12*n2*u1*v1*v3

             2    3             3        2    3
      - 4*a33 *b33 *k12*n2*u1*v2  + 4*a33 *b33 *k12*n2*u1*v3

             2    3          3           2    3          3
      - 4*a33 *b33 *k12*n2*u2 *v1 + 8*a33 *b33 *k12*n2*u3 *v3

             2    3             2        2    3          3
      + 4*a33 *b33 *k12*n2*u3*v2  - 4*a33 *b33 *k12*n2*v1

             2    3                     2    3                2
      + 4*a33 *b33 *k12*n2*v1*v3 + 4*a33 *b33 *k16*n2*u1*v1*v3

             2    3                2        2    3             3
      + 4*a33 *b33 *k16*n2*u2*v2*v3  + 4*a33 *b33 *k16*n2*u2*v3

             2    3          2   2        2    3             3
      + 8*a33 *b33 *k16*n2*v1 *v3  + 4*a33 *b33 *k16*n2*v1*v3

             2    2      2   2           2    2      2
      - 8*a33 *b33 *k1*n2 *u1 *u3 - 8*a33 *b33 *k1*n2 *u2*u3*v2

             2    2      2               2    2      2   2
      - 8*a33 *b33 *k1*n2 *u3*v2 + 16*a33 *b33 *k1*n2 *v1 *v2

             2    2       2   2           2    2       2
      - 4*a33 *b33 *k12*n2 *u1 *u3 - 4*a33 *b33 *k12*n2 *u1*u3*v3

             2    2       2                 2    2       2   2
      - 4*a33 *b33 *k12*n2 *u2*u3*v2 + 8*a33 *b33 *k12*n2 *u3

             2    2       2   2           2    2       2
      - 8*a33 *b33 *k12*n2 *v1 *v2 - 4*a33 *b33 *k16*n2 *u3*v2*v3

             2    2       2      2        2    2       2      2
      - 4*a33 *b33 *k16*n2 *u3*v3  - 4*a33 *b33 *k16*n2 *v2*v3

             2    2       2      2        2    2       2   2
      - 4*a33 *b33 *k38*n2 *u1*v1  - 4*a33 *b33 *k38*n2 *v1 *v3

             2    2       2      2        2          3
      - 4*a33 *b33 *k38*n2 *v1*v2  + 4*a33 *b33*k1*n2 *u1*v1

             2           3                  4
      - 4*a33 *b33*k12*n2 *u1*v1 + 8*a33*b33 *k1*n2*u1*u2*v1*v2

                 4               2            4            2
      - 8*a33*b33 *k1*n2*u1*u3*v3  + 8*a33*b33 *k1*n2*u1*v1 *v2

                 4         2                  4
      + 8*a33*b33 *k1*n2*u2 *v1*v3 + 8*a33*b33 *k1*n2*u2*u3*v1*v2

                 4            2               4         3
      + 8*a33*b33 *k1*n2*u2*v1 *v3 + 8*a33*b33 *k1*n2*v1 *v3

                 4            2               4
      - 8*a33*b33 *k1*n2*v1*v2 *v3 - 4*a33*b33 *k12*n2*u1*u2*v1*v2

                 4             2               4          2
      - 4*a33*b33 *k12*n2*u1*v1 *v2 - 4*a33*b33 *k12*n2*u2 *v1*v3

                 4                               4             2
      - 4*a33*b33 *k12*n2*u2*u3*v1*v2 - 4*a33*b33 *k12*n2*u2*v1 *v3

                 4          3               3      2   2
      - 4*a33*b33 *k12*n2*v1 *v3 - 4*a33*b33 *k1*n2 *u1 *v3

                 3      2      2            3      2   2
      - 8*a33*b33 *k1*n2 *u1*v2  + 8*a33*b33 *k1*n2 *v1 *v3

                 3      2      2            3       2      2
      - 8*a33*b33 *k1*n2 *v1*v2  + 2*a33*b33 *k12*n2 *u1*v2

                 3       2                  3       2   2
      + 4*a33*b33 *k12*n2 *u3*v3 + 4*a33*b33 *k12*n2 *v1 *v3

                 3       2      2            2      3   2
      + 2*a33*b33 *k12*n2 *v1*v2  + 4*a33*b33 *k1*n2 *v1

             5               2        5               2        5            3
      + 3*b33 *k1*n2*u1*v1*v3  + 2*b33 *k1*n2*u2*v2*v3  + 3*b33 *k1*n2*u2*v3

             5         2   2        5            3      5                2
      + 6*b33 *k1*n2*v1 *v3  + 2*b33 *k1*n2*v1*v3  - b33 *k12*n2*u1*v1*v3

           5                2      5             3        5          2   2
      - b33 *k12*n2*u2*v2*v3  - b33 *k12*n2*u2*v3  - 2*b33 *k12*n2*v1 *v3

           5             3        4      2                 4      2      2
      - b33 *k12*n2*v1*v3  - 4*b33 *k1*n2 *u3*v2*v3 - 4*b33 *k1*n2 *u3*v3

             4      2      2      4       2               4       2      2
      - 2*b33 *k1*n2 *v2*v3  + b33 *k12*n2 *u3*v2*v3 + b33 *k12*n2 *u3*v3

             4       2      2      3      3            3       3               4
      + 2*b33 *k12*n2 *v2*v3  + b33 *k1*n2 *v1*v3 + b33 *k12*n2 *v1*v3)/(16*a33

       *b33*n2)



And again in machine readable form:



HAM=( - 4*a33**2*u1**2 - 4*a33**2*u2**2 + 4*a33**2*u3**2 + 4*a33*b33*u3*v3 + 4*
a33*n2*u2 + b33**2*v3**2)/(4*a33)$

INT=( - 64*a33**6*r5*u1**2 - 64*a33**6*r5*u1*v1 - 64*a33**5*b33*r5*u1*u3 - 64*
a33**5*b33*r5*v1**2 + 32*a33**5*n2*r5*u2 - 16*a33**4*b33**2*r5*u1*v3 - 16*a33**4
*b33**2*r5*v1*v3 + 16*a33**4*b33*k1*n2*u1**4 + 32*a33**4*b33*k1*n2*u1**2*u2*u3 +
 80*a33**4*b33*k1*n2*u1**2 + 16*a33**4*b33*k1*n2*u1*u2*v2**2 - 64*a33**4*b33*k1*
n2*u1*u3**2 + 16*a33**4*b33*k1*n2*u1*v1 + 16*a33**4*b33*k12*n2*u1**2 + 16*a33**4
*b33*k12*n2*u1*v1*v2**2 + 16*a33**4*b33*k12*n2*u1*v1 + 16*a33**4*b33*k12*n2*u2**
2*u3**2 + 16*a33**4*b33*k16*n2*u1**3*v1 + 16*a33**4*b33*k16*n2*u1**3*v3 + 16*a33
**4*b33*k16*n2*u1*u2*u3*v1 + 16*a33**4*b33*k16*n2*u1*v2**3 + 16*a33**4*b33*k16*
n2*u2**3*v1 + 16*a33**4*b33*k16*n2*v1*v2**3 + 16*a33**4*b33*k38*n2*u1*u2**2*v2 -
 32*a33**4*b33*k38*n2*u1*u3*v1 - 16*a33**4*b33*k38*n2*u1*u3 + 16*a33**4*b33*k38*
n2*u2**3*v3 + 16*a33**4*b33*k38*n2*u2*u3*v2*v3 + 32*a33**4*b33*k38*n2*u2*v1*v2 -
 16*a33**4*b33*k38*n2*v1**2 + 16*a33**4*b33*n2*r5*v2 - 32*a33**3*b33**2*k1*n2*u1
*u3*v1 - 16*a33**3*b33**2*k1*n2*u1*u3 + 32*a33**3*b33**2*k1*n2*u2*v1*v2 - 32*a33
**3*b33**2*k1*n2*u3**4 - 32*a33**3*b33**2*k1*n2*v1**2*v2**2 - 16*a33**3*b33**2*
k1*n2*v1**2 + 16*a33**3*b33**2*k12*n2*u1*u3 - 16*a33**3*b33**2*k12*n2*u2*u3*v2*
v3 + 16*a33**3*b33**2*k12*n2*v1**2 + 16*a33**3*b33**2*k16*n2*u1*u2*v1*v2 + 16*
a33**3*b33**2*k16*n2*u1*v1**2*v2 + 16*a33**3*b33**2*k16*n2*u2**2*v1*v3 + 16*a33
**3*b33**2*k16*n2*u2*u3*v1*v2 + 16*a33**3*b33**2*k16*n2*u2*v1**2*v3 + 16*a33**3*
b33**2*k16*n2*v1**3*v3 - 16*a33**3*b33*k1*n2**2*u1**2*u2 - 16*a33**3*b33*k1*n2**
2*u1*u2*v2 + 16*a33**3*b33*k1*n2**2*u2**2 - 32*a33**3*b33*k1*n2**2*u2 + 8*a33**3
*b33*k12*n2**2*u1**2*u2 + 8*a33**3*b33*k12*n2**2*u1*u2*v2 - 8*a33**3*b33*k12*n2
**2*u1*v1*v2 - 8*a33**3*b33*k12*n2**2*u2 - 8*a33**3*b33*k16*n2**2*u1**2*v3 - 8*
a33**3*b33*k16*n2**2*u1*v2**2 - 8*a33**3*b33*k16*n2**2*v1*v2**2 + 8*a33**3*b33*
k38*n2**2*u2*u3**2 + 8*a33**3*b33*k38*n2**2*u2*u3*v2 - 8*a33**3*b33*k38*n2**2*u3
*v2 + 8*a33**3*b33*k38*n2**2*v1**2*v2 + 8*a33**2*b33**3*k1*n2*u1**3*v3 + 8*a33**
2*b33**3*k1*n2*u1*u2*u3*v1 - 32*a33**2*b33**3*k1*n2*u1*u3**2*v2 - 32*a33**2*b33
**3*k1*n2*u1*v1**3 + 16*a33**2*b33**3*k1*n2*u1*v1*v3 + 16*a33**2*b33**3*k1*n2*u1
*v2**3 - 12*a33**2*b33**3*k1*n2*u1*v3 + 16*a33**2*b33**3*k1*n2*u2**3*v1 - 32*a33
**2*b33**3*k1*n2*u3**3*v3 - 12*a33**2*b33**3*k1*n2*v1*v3 - 4*a33**2*b33**3*k12*
n2*u1**3*v3 - 4*a33**2*b33**3*k12*n2*u1*u2*u3*v1 + 8*a33**2*b33**3*k12*n2*u1*u3
**2*v2 + 8*a33**2*b33**3*k12*n2*u1*v1**3 - 4*a33**2*b33**3*k12*n2*u1*v1*v3 - 4*
a33**2*b33**3*k12*n2*u1*v2**3 + 4*a33**2*b33**3*k12*n2*u1*v3 - 4*a33**2*b33**3*
k12*n2*u2**3*v1 + 8*a33**2*b33**3*k12*n2*u3**3*v3 + 4*a33**2*b33**3*k12*n2*u3*v2
**2 - 4*a33**2*b33**3*k12*n2*v1**3 + 4*a33**2*b33**3*k12*n2*v1*v3 + 4*a33**2*b33
**3*k16*n2*u1*v1*v3**2 + 4*a33**2*b33**3*k16*n2*u2*v2*v3**2 + 4*a33**2*b33**3*
k16*n2*u2*v3**3 + 8*a33**2*b33**3*k16*n2*v1**2*v3**2 + 4*a33**2*b33**3*k16*n2*v1
*v3**3 - 8*a33**2*b33**2*k1*n2**2*u1**2*u3 - 8*a33**2*b33**2*k1*n2**2*u2*u3*v2 -
 8*a33**2*b33**2*k1*n2**2*u3*v2 + 16*a33**2*b33**2*k1*n2**2*v1**2*v2 - 4*a33**2*
b33**2*k12*n2**2*u1**2*u3 - 4*a33**2*b33**2*k12*n2**2*u1*u3*v3 - 4*a33**2*b33**2
*k12*n2**2*u2*u3*v2 + 8*a33**2*b33**2*k12*n2**2*u3**2 - 8*a33**2*b33**2*k12*n2**
2*v1**2*v2 - 4*a33**2*b33**2*k16*n2**2*u3*v2*v3 - 4*a33**2*b33**2*k16*n2**2*u3*
v3**2 - 4*a33**2*b33**2*k16*n2**2*v2*v3**2 - 4*a33**2*b33**2*k38*n2**2*u1*v1**2 
- 4*a33**2*b33**2*k38*n2**2*v1**2*v3 - 4*a33**2*b33**2*k38*n2**2*v1*v2**2 + 4*
a33**2*b33*k1*n2**3*u1*v1 - 4*a33**2*b33*k12*n2**3*u1*v1 + 8*a33*b33**4*k1*n2*u1
*u2*v1*v2 - 8*a33*b33**4*k1*n2*u1*u3*v3**2 + 8*a33*b33**4*k1*n2*u1*v1**2*v2 + 8*
a33*b33**4*k1*n2*u2**2*v1*v3 + 8*a33*b33**4*k1*n2*u2*u3*v1*v2 + 8*a33*b33**4*k1*
n2*u2*v1**2*v3 + 8*a33*b33**4*k1*n2*v1**3*v3 - 8*a33*b33**4*k1*n2*v1*v2**2*v3 - 
4*a33*b33**4*k12*n2*u1*u2*v1*v2 - 4*a33*b33**4*k12*n2*u1*v1**2*v2 - 4*a33*b33**4
*k12*n2*u2**2*v1*v3 - 4*a33*b33**4*k12*n2*u2*u3*v1*v2 - 4*a33*b33**4*k12*n2*u2*
v1**2*v3 - 4*a33*b33**4*k12*n2*v1**3*v3 - 4*a33*b33**3*k1*n2**2*u1**2*v3 - 8*a33
*b33**3*k1*n2**2*u1*v2**2 + 8*a33*b33**3*k1*n2**2*v1**2*v3 - 8*a33*b33**3*k1*n2
**2*v1*v2**2 + 2*a33*b33**3*k12*n2**2*u1*v2**2 + 4*a33*b33**3*k12*n2**2*u3*v3 + 
4*a33*b33**3*k12*n2**2*v1**2*v3 + 2*a33*b33**3*k12*n2**2*v1*v2**2 + 4*a33*b33**2
*k1*n2**3*v1**2 + 3*b33**5*k1*n2*u1*v1*v3**2 + 2*b33**5*k1*n2*u2*v2*v3**2 + 3*
b33**5*k1*n2*u2*v3**3 + 6*b33**5*k1*n2*v1**2*v3**2 + 2*b33**5*k1*n2*v1*v3**3 - 
b33**5*k12*n2*u1*v1*v3**2 - b33**5*k12*n2*u2*v2*v3**2 - b33**5*k12*n2*u2*v3**3 -
 2*b33**5*k12*n2*v1**2*v3**2 - b33**5*k12*n2*v1*v3**3 - 4*b33**4*k1*n2**2*u3*v2*
v3 - 4*b33**4*k1*n2**2*u3*v3**2 - 2*b33**4*k1*n2**2*v2*v3**2 + b33**4*k12*n2**2*
u3*v2*v3 + b33**4*k12*n2**2*u3*v3**2 + 2*b33**4*k12*n2**2*v2*v3**2 + b33**3*k1*
n2**3*v1*v3 + b33**3*k12*n2**3*v1*v3)/(16*a33**4*b33*n2)$