Solution 22 to problem e3quant


Expressions | Parameters | Relevance | Back to problem e3quant

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


n2=0


n3=0


        1                1
     - ---*a22*b33*n1 - ---*a33*b33*n1
        2                2
m1=------------------------------------
                 2
              a22  - a22*a33


m2=0


m3=0


r6=0


r5=0


r3=0


r2=0


            2
     - 2*a22 *r4 + 2*a22*a33*r4 + 6*b33*k1*n1
r1=-------------------------------------------
                a22*b33 - a33*b33


      1     2
     ---*b33 *k12*n1
      2
q20=-----------------
        3      2
     a22  - a22 *a33


         1     2           1                    3     2          1     2
      - ---*a22 *b33*r4 + ---*a22*a33*b33*r4 - ---*b33 *k1*n1 + ---*b33 *k12*n1
         2                 2                    4                4
q19=----------------------------------------------------------------------------
                                         2
                                      a22 *n1


q18=0


q17=0


         1     4           3     3               3     2    2
q16=( - ---*a22 *b33*r4 + ---*a22 *a33*b33*r4 - ---*a22 *a33 *b33*r4
         2                 2                     2

         3     2    2          1     2    2           1         3
      - ---*a22 *b33 *k1*n1 + ---*a22 *b33 *k12*n1 + ---*a22*a33 *b33*r4
         4                     4                      2

         3             2          1             2           3     2    2
      + ---*a22*a33*b33 *k1*n1 - ---*a22*a33*b33 *k12*n1 - ---*a33 *b33 *k1*n1
         2                        2                         4

         1     2    2           1     2      3    1     2       3      4
      + ---*a33 *b33 *k12*n1 + ---*b33 *k1*n1  + ---*b33 *k12*n1 )/(a22 *n1
         4                      4                 4

           3             2    2
    - 2*a22 *a33*n1 + a22 *a33 *n1)


       b33*k1*n1
q14=----------------
        2
     a22  - a22*a33


q13=0


      - b33*k12*n1
q11=----------------
        2
     a22  - a22*a33


          2
     2*a22 *r4 - 2*a22*a33*r4 - a22*k38*n1 + b33*k1*n1 - b33*k12*n1
q10=----------------------------------------------------------------
                                 a22*n1


q9=0


     - 2*k1*n1 + 2*k12*n1
q8=-----------------------
          a22 - a33


            2
     - 2*a22 *r4 + 2*a22*a33*r4 + b33*k1*n1 + b33*k12*n1
q7=------------------------------------------------------
                           b33*n1


q6=0


q5=0


         4           3             3               2    2           2
q4=(2*a22 *r4 - 6*a22 *a33*r4 - a22 *k38*n1 + 6*a22 *a33 *r4 + 2*a22 *a33*k38*n1

          2                2                       3             2
     + a22 *b33*k1*n1 - a22 *b33*k12*n1 - 2*a22*a33 *r4 - a22*a33 *k38*n1

                                                       2
     - 2*a22*a33*b33*k1*n1 + 2*a22*a33*b33*k12*n1 + a33 *b33*k1*n1

          2                       3      3           2                 2
     - a33 *b33*k12*n1 - b33*k1*n1 )/(a22 *n1 - 2*a22 *a33*n1 + a22*a33 *n1)


q3=0


q2=0


            5           4               3                3
q1=( - 2*a22 *r4 + 4*a22 *a33*r4 + 5*a22 *b33*k1*n1 + a22 *b33*k12*n1

            2    3           2                    2
     - 4*a22 *a33 *r4 - 5*a22 *a33*b33*k1*n1 - a22 *a33*b33*k12*n1

                4               2                    2
     + 2*a22*a33 *r4 - 5*a22*a33 *b33*k1*n1 - a22*a33 *b33*k12*n1

                    3                 3        3                3
     + a22*b33*k1*n1  - a22*b33*k12*n1  + 5*a33 *b33*k1*n1 + a33 *b33*k12*n1

                    3                 3      3             2
     + a33*b33*k1*n1  - a33*b33*k12*n1 )/(a22 *b33*n1 - a22 *a33*b33*n1

             2             3
    - a22*a33 *b33*n1 + a33 *b33*n1)


p56=0


p55=0


p54=0


p53=0


         1     3               1     2                   1         3
p52=( - ---*a22 *b33*k16*n1 + ---*a22 *a33*b33*k16*n1 - ---*a22*b33 *k1*n1
         2                     2                         4

         1         3           1         3          1         3             5
      + ---*a22*b33 *k12*n1 + ---*a33*b33 *k1*n1 - ---*a33*b33 *k12*n1)/(a22
         4                     4                    4

           4          3    2
    - 2*a22 *a33 + a22 *a33 )


p51=0


         1     3               1     2                   1         3
p50=( - ---*a22 *b33*k16*n1 + ---*a22 *a33*b33*k16*n1 - ---*a22*b33 *k1*n1
         2                     2                         2

         1         3           1         3          1         3             5
      + ---*a22*b33 *k12*n1 + ---*a33*b33 *k1*n1 - ---*a33*b33 *k12*n1)/(a22
         8                     2                    8

           4          3    2
    - 2*a22 *a33 + a22 *a33 )


p49=0


p48=0


         1     3               1     2                   1         3
p47=( - ---*a22 *b33*k16*n1 + ---*a22 *a33*b33*k16*n1 - ---*a22*b33 *k1*n1
         2                     2                         2

         1         3           1         3          1         3             5
      + ---*a22*b33 *k12*n1 + ---*a33*b33 *k1*n1 - ---*a33*b33 *k12*n1)/(a22
         8                     2                    8

           4          3    2
    - 2*a22 *a33 + a22 *a33 )


p46=0


p45=0


p44=0


         1     2               1                              2
p43=( - ---*a22 *b33*k38*n1 + ---*a22*a33*b33*k38*n1 - a22*b33 *k1*n1
         2                     2

         1         2                 2          1         2             4
      - ---*a22*b33 *k12*n1 + a33*b33 *k1*n1 + ---*a33*b33 *k12*n1)/(a22
         2                                      2

           3          2    2
    - 2*a22 *a33 + a22 *a33 )


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


p41=0


p40=0


p39=0


         1
      - ---*b33*k12*n1
         2
p38=-------------------
         2
      a22  - a22*a33


p37=0


p36=0


p35=0


p34=0


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


         1
      - ---*b33*k38*n1
         2
p32=-------------------
         2
      a22  - a22*a33


p31=0


p30=0


p29=0


     2*a22*k38 - 2*b33*k1
p28=----------------------
             a22


p27=0


p26=0


p25=0


                     1
      - b33*k1*n1 - ---*b33*k12*n1
                     2
p24=-------------------------------
               2
            a22  - a22*a33


p23=0


p22=0


        3             2               1         2          1         2
     a22 *k16*n1 - a22 *a33*k16*n1 + ---*a22*b33 *k1*n1 - ---*a33*b33 *k1*n1
                                      2                    2
p21=-------------------------------------------------------------------------
                             4        3          2    2
                          a22  - 2*a22 *a33 + a22 *a33


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


        3             2                     2          1         2
p19=(a22 *k16*n1 - a22 *a33*k16*n1 + a22*b33 *k1*n1 - ---*a22*b33 *k12*n1
                                                       4

               2          1         2             4        3          2    2
      - a33*b33 *k1*n1 + ---*a33*b33 *k12*n1)/(a22  - 2*a22 *a33 + a22 *a33 )
                          4


p18=0


p17=0


        3             2               1     2
p16=(a22 *k16*n1 - a22 *a33*k16*n1 - ---*a22 *b33*k38*n1
                                      2

         1                              2          1         2
      + ---*a22*a33*b33*k38*n1 + a22*b33 *k1*n1 - ---*a22*b33 *k12*n1
         2                                         4

               2          1         2             4        3          2    2
      - a33*b33 *k1*n1 + ---*a33*b33 *k12*n1)/(a22  - 2*a22 *a33 + a22 *a33 )
                          4


           2
p15=( - a22 *k38*n1 + a22*a33*k38*n1 + 2*a22*b33*k1*n1 - a22*b33*k12*n1

                                              3        2              2
      - 2*a33*b33*k1*n1 + a33*b33*k12*n1)/(a22  - 2*a22 *a33 + a22*a33 )


      - 2*a22*k38 + 2*b33*k1
p14=-------------------------
               a22


p13=0


      k12*n1
p12=-----------
     a22 - a33


p11=0


      - k38*n1
p10=-----------
     a22 - a33


p9=0


p8= - 4*k1


    2*k1*n1 - k12*n1
p7=------------------
       a22 - a33


p6=0


p5=0


          2                                            1
p4=( - a22 *k38*n1 + a22*a33*k38*n1 - a22*b33*k1*n1 - ---*a22*b33*k12*n1
                                                       2

                        1                      3        2              2
     + a33*b33*k1*n1 + ---*a33*b33*k12*n1)/(a22  - 2*a22 *a33 + a22*a33 )
                        2


p3=0


p2=0


    2*k1*n1 - k12*n1
p1=------------------
       a22 - a33


k125=0


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


k123=0


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


k121=0


k120=0


k119=0


k118=0


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


k116=0


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


k114=0


k113=0


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


k110=0


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


k108=0


k107=0


k106=0


k105=0


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


k103=0


k102=0


k100=0


k99=0


k98=0


k97=0


     a22*k38 + b33*k12
k95=-------------------
            a22


k94=0


k93=0


k91=0


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


k89=0


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


k87=0


k86=0


k85=0


k84=0


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


k82=0


k81=0


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


k79=0


k78=0


k77=0


k76=0


k75=0


k74=k38


k73=0


k72=0


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


k70=0


k69=k16


k68=0


k67=0


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


     2*b33*k1
k65=----------
       a22


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
      - a22 *b33*k16 - ---*b33 *k1 + ---*b33 *k12
                        2             4
k52=----------------------------------------------
                            3
                         a22


k51=0


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


k49=0


k48=0


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


k46=0


k45=0


k44=0


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


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
                a22


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
     a22 *k16 + ---*b33 *k1 - ---*b33 *k12
                 2             4
k21=---------------------------------------
                        2
                     a22


k20=0


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


k18=0


k17=0


     2*b33*k1
k15=----------
       a22


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:
 b33,r4,k38,k16,k12,k1,n1,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*a33 *u3  - 4*a22*a33*b33*u3*v3

                                 2   2                            2   2
      - 4*a22*a33*n1*u1 - a22*b33 *v3  - 2*a22*b33*n1*v1 + a33*b33 *v3

      - 2*a33*b33*n1*v1)/(4*a22*(a22 - a33))

              9      2         9                  8          2
INT=( - 32*a22 *r4*u1  - 32*a22 *r4*u1*v1 + 64*a22 *a33*r4*u1

              8                      8                      8          2
      + 64*a22 *a33*r4*u1*v1 + 32*a22 *b33*r4*u1*u3 + 32*a22 *b33*r4*v1

              8                  7                          7              2
      - 32*a22 *n1*r4*u1 - 64*a22 *a33*b33*r4*u1*u3 - 64*a22 *a33*b33*r4*v1

              7                     7    2                 7    2
      + 32*a22 *a33*n1*r4*u1 - 8*a22 *b33 *r4*u1*v3 - 8*a22 *b33 *r4*v1*v3

              7             4         7             2
      + 16*a22 *b33*k1*n1*u1  + 32*a22 *b33*k1*n1*u1 *u2*u3

              7             2         7                   2
      + 80*a22 *b33*k1*n1*u1  + 16*a22 *b33*k1*n1*u1*u2*v2

              7                2         7
      - 64*a22 *b33*k1*n1*u1*u3  + 16*a22 *b33*k1*n1*u1*v1

              7              2         7                    2
      + 16*a22 *b33*k12*n1*u1  + 16*a22 *b33*k12*n1*u1*v1*v2

              7                          7              2   2
      + 16*a22 *b33*k12*n1*u1*v1 + 16*a22 *b33*k12*n1*u2 *u3

              7              3            7              3
      + 16*a22 *b33*k16*n1*u1 *v1 + 16*a22 *b33*k16*n1*u1 *v3

              7                                7                 3
      + 16*a22 *b33*k16*n1*u1*u2*u3*v1 + 16*a22 *b33*k16*n1*u1*v2

              7              3            7                 3
      + 16*a22 *b33*k16*n1*u2 *v1 + 16*a22 *b33*k16*n1*v1*v2

              7                 2            7
      + 16*a22 *b33*k38*n1*u1*u2 *v2 - 32*a22 *b33*k38*n1*u1*u3*v1

              7                          7              3
      - 16*a22 *b33*k38*n1*u1*u3 + 16*a22 *b33*k38*n1*u2 *v3

              7                                7
      + 16*a22 *b33*k38*n1*u2*u3*v2*v3 + 32*a22 *b33*k38*n1*u2*v1*v2

              7              2         7                      6    3      2
      - 16*a22 *b33*k38*n1*v1  + 16*a22 *b33*n1*r4*v1 - 64*a22 *a33 *r4*u1

              6    3                  6    2
      - 64*a22 *a33 *r4*u1*v1 + 32*a22 *a33 *n1*r4*u1

              6        2                  6        2
      + 16*a22 *a33*b33 *r4*u1*v3 + 16*a22 *a33*b33 *r4*v1*v3

              6                 4         6                 2
      - 16*a22 *a33*b33*k1*n1*u1  - 32*a22 *a33*b33*k1*n1*u1 *u2*u3

              6                 2         6                       2
      - 80*a22 *a33*b33*k1*n1*u1  - 16*a22 *a33*b33*k1*n1*u1*u2*v2

              6                    2         6
      + 64*a22 *a33*b33*k1*n1*u1*u3  - 16*a22 *a33*b33*k1*n1*u1*v1

              6                  2         6                        2
      - 16*a22 *a33*b33*k12*n1*u1  - 16*a22 *a33*b33*k12*n1*u1*v1*v2

              6                              6                  2   2
      - 16*a22 *a33*b33*k12*n1*u1*v1 - 16*a22 *a33*b33*k12*n1*u2 *u3

              6                  3            6                  3
      - 16*a22 *a33*b33*k16*n1*u1 *v1 - 16*a22 *a33*b33*k16*n1*u1 *v3

              6                                    6                     3
      - 16*a22 *a33*b33*k16*n1*u1*u2*u3*v1 - 16*a22 *a33*b33*k16*n1*u1*v2

              6                  3            6                     3
      - 16*a22 *a33*b33*k16*n1*u2 *v1 - 16*a22 *a33*b33*k16*n1*v1*v2

              6                     2            6
      - 16*a22 *a33*b33*k38*n1*u1*u2 *v2 + 32*a22 *a33*b33*k38*n1*u1*u3*v1

              6                              6                  3
      + 16*a22 *a33*b33*k38*n1*u1*u3 - 16*a22 *a33*b33*k38*n1*u2 *v3

              6                                    6
      - 16*a22 *a33*b33*k38*n1*u2*u3*v2*v3 - 32*a22 *a33*b33*k38*n1*u2*v1*v2

              6                  2         6
      + 16*a22 *a33*b33*k38*n1*v1  - 16*a22 *a33*b33*n1*r4*v1

              6    2                        6    2
      + 32*a22 *b33 *k1*n1*u1*u3*v1 + 16*a22 *b33 *k1*n1*u1*u3

              6    2                        6    2         4
      - 32*a22 *b33 *k1*n1*u2*v1*v2 + 32*a22 *b33 *k1*n1*u3

              6    2         2   2         6    2         2
      + 32*a22 *b33 *k1*n1*v1 *v2  + 16*a22 *b33 *k1*n1*v1

              6    2                      6    2
      - 16*a22 *b33 *k12*n1*u1*u3 + 16*a22 *b33 *k12*n1*u2*u3*v2*v3

              6    2          2         6    2
      - 16*a22 *b33 *k12*n1*v1  - 16*a22 *b33 *k16*n1*u1*u2*v1*v2

              6    2             2            6    2          2
      - 16*a22 *b33 *k16*n1*u1*v1 *v2 - 16*a22 *b33 *k16*n1*u2 *v1*v3

              6    2                            6    2             2
      - 16*a22 *b33 *k16*n1*u2*u3*v1*v2 - 16*a22 *b33 *k16*n1*u2*v1 *v3

              6    2          3            6          2   3
      - 16*a22 *b33 *k16*n1*v1 *v3 + 32*a22 *b33*k1*n1 *u1

              6          2            6          2   2
      + 96*a22 *b33*k1*n1 *u1 + 32*a22 *b33*k1*n1 *u2 *u3

              6          2               6           2   3
      - 32*a22 *b33*k1*n1 *u2*v1 - 16*a22 *b33*k12*n1 *u1

              6           2                  6           2   2
      + 16*a22 *b33*k12*n1 *u1*u2*v1 - 16*a22 *b33*k12*n1 *u2 *u3

              6           2               6           2   2
      + 32*a22 *b33*k12*n1 *u2*v1 + 16*a22 *b33*k16*n1 *u1 *v2

              6           2   2            6           2      2
      + 16*a22 *b33*k16*n1 *u3 *v1 + 16*a22 *b33*k16*n1 *u3*v1

              6           2   3         6           2
      - 16*a22 *b33*k38*n1 *u2  - 16*a22 *b33*k38*n1 *u2*u3*v1

              6           2   3         5    4      2         5    4
      - 16*a22 *b33*k38*n1 *u3  + 32*a22 *a33 *r4*u1  + 32*a22 *a33 *r4*u1*v1

              5    3                      5    3          2
      + 64*a22 *a33 *b33*r4*u1*u3 + 64*a22 *a33 *b33*r4*v1

              5    3                  5    2             4
      - 32*a22 *a33 *n1*r4*u1 - 16*a22 *a33 *b33*k1*n1*u1

              5    2             2               5    2             2
      - 32*a22 *a33 *b33*k1*n1*u1 *u2*u3 - 80*a22 *a33 *b33*k1*n1*u1

              5    2                   2         5    2                2
      - 16*a22 *a33 *b33*k1*n1*u1*u2*v2  + 64*a22 *a33 *b33*k1*n1*u1*u3

              5    2                         5    2              2
      - 16*a22 *a33 *b33*k1*n1*u1*v1 - 16*a22 *a33 *b33*k12*n1*u1

              5    2                    2         5    2
      - 16*a22 *a33 *b33*k12*n1*u1*v1*v2  - 16*a22 *a33 *b33*k12*n1*u1*v1

              5    2              2   2         5    2              3
      - 16*a22 *a33 *b33*k12*n1*u2 *u3  - 16*a22 *a33 *b33*k16*n1*u1 *v1

              5    2              3            5    2
      - 16*a22 *a33 *b33*k16*n1*u1 *v3 - 16*a22 *a33 *b33*k16*n1*u1*u2*u3*v1

              5    2                 3         5    2              3
      - 16*a22 *a33 *b33*k16*n1*u1*v2  - 16*a22 *a33 *b33*k16*n1*u2 *v1

              5    2                 3         5    2                 2
      - 16*a22 *a33 *b33*k16*n1*v1*v2  - 16*a22 *a33 *b33*k38*n1*u1*u2 *v2

              5    2                             5    2
      + 32*a22 *a33 *b33*k38*n1*u1*u3*v1 + 16*a22 *a33 *b33*k38*n1*u1*u3

              5    2              3            5    2
      - 16*a22 *a33 *b33*k38*n1*u2 *v3 - 16*a22 *a33 *b33*k38*n1*u2*u3*v2*v3

              5    2                             5    2              2
      - 32*a22 *a33 *b33*k38*n1*u2*v1*v2 + 16*a22 *a33 *b33*k38*n1*v1

              5    2                      5        2
      - 16*a22 *a33 *b33*n1*r4*v1 - 32*a22 *a33*b33 *k1*n1*u1*u3*v1

              5        2                     5        2
      - 16*a22 *a33*b33 *k1*n1*u1*u3 + 32*a22 *a33*b33 *k1*n1*u2*v1*v2

              5        2         4         5        2         2   2
      - 32*a22 *a33*b33 *k1*n1*u3  - 32*a22 *a33*b33 *k1*n1*v1 *v2

              5        2         2         5        2
      - 16*a22 *a33*b33 *k1*n1*v1  + 16*a22 *a33*b33 *k12*n1*u1*u3

              5        2                            5        2          2
      - 16*a22 *a33*b33 *k12*n1*u2*u3*v2*v3 + 16*a22 *a33*b33 *k12*n1*v1

              5        2                            5        2             2
      + 16*a22 *a33*b33 *k16*n1*u1*u2*v1*v2 + 16*a22 *a33*b33 *k16*n1*u1*v1 *v2

              5        2          2               5        2
      + 16*a22 *a33*b33 *k16*n1*u2 *v1*v3 + 16*a22 *a33*b33 *k16*n1*u2*u3*v1*v2

              5        2             2            5        2          3
      + 16*a22 *a33*b33 *k16*n1*u2*v1 *v3 + 16*a22 *a33*b33 *k16*n1*v1 *v3

             5    3         3           5    3
      + 8*a22 *b33 *k1*n1*u1 *v3 + 8*a22 *b33 *k1*n1*u1*u2*u3*v1

              5    3            2            5    3            3
      - 32*a22 *b33 *k1*n1*u1*u3 *v2 - 32*a22 *b33 *k1*n1*u1*v1

              5    3                        5    3            3
      + 16*a22 *b33 *k1*n1*u1*v1*v3 + 16*a22 *b33 *k1*n1*u1*v2

              5    3                     5    3         3
      - 12*a22 *b33 *k1*n1*u1*v3 + 16*a22 *b33 *k1*n1*u2 *v1

              5    3         3            5    3
      - 32*a22 *b33 *k1*n1*u3 *v3 - 12*a22 *b33 *k1*n1*v1*v3

             5    3          3           5    3
      - 4*a22 *b33 *k12*n1*u1 *v3 - 4*a22 *b33 *k12*n1*u1*u2*u3*v1

             5    3             2           5    3             3
      + 8*a22 *b33 *k12*n1*u1*u3 *v2 + 8*a22 *b33 *k12*n1*u1*v1

             5    3                        5    3             3
      - 4*a22 *b33 *k12*n1*u1*v1*v3 - 4*a22 *b33 *k12*n1*u1*v2

             5    3                     5    3          3
      + 4*a22 *b33 *k12*n1*u1*v3 - 4*a22 *b33 *k12*n1*u2 *v1

             5    3          3           5    3             2
      + 8*a22 *b33 *k12*n1*u3 *v3 + 4*a22 *b33 *k12*n1*u3*v2

             5    3          3        5    3
      - 4*a22 *b33 *k12*n1*v1  + 4*a22 *b33 *k12*n1*v1*v3

             5    3                2        5    3                2
      + 4*a22 *b33 *k16*n1*u1*v1*v3  + 4*a22 *b33 *k16*n1*u2*v2*v3

             5    3             3        5    3          2   2
      + 4*a22 *b33 *k16*n1*u2*v3  + 8*a22 *b33 *k16*n1*v1 *v3

             5    3             3         5    2      2
      + 4*a22 *b33 *k16*n1*v1*v3  - 16*a22 *b33 *k1*n1 *u1*u3*v2

              5    2      2   3         5    2      2
      - 16*a22 *b33 *k1*n1 *u2  + 32*a22 *b33 *k1*n1 *u2*u3*v1

              5    2      2              5    2       2
      + 16*a22 *b33 *k1*n1 *v1*v2 - 8*a22 *b33 *k12*n1 *u1*u3*v2

              5    2       2              5    2       2   3
      - 16*a22 *b33 *k12*n1 *u1*v2 - 8*a22 *b33 *k12*n1 *u2

             5    2       2   2            5    2       2
      - 8*a22 *b33 *k12*n1 *u2 *v3 - 16*a22 *b33 *k12*n1 *u2*u3*v1

             5    2       2                 5    2       2      2
      - 8*a22 *b33 *k16*n1 *u2*v2*v3 - 8*a22 *b33 *k16*n1 *u2*v3

             5    2       2   2           5    2       2
      - 8*a22 *b33 *k16*n1 *v2 *v3 - 8*a22 *b33 *k38*n1 *u2*v1*v3

             5    2       2      2        5    2       2   2
      - 8*a22 *b33 *k38*n1 *u2*v2  - 8*a22 *b33 *k38*n1 *u3 *v1

              5          3   2         5           3   2
      + 16*a22 *b33*k1*n1 *u1  - 16*a22 *b33*k12*n1 *u1

              4    4                      4    4          2
      - 32*a22 *a33 *b33*r4*u1*u3 - 32*a22 *a33 *b33*r4*v1

              4    3    2                  4    3    2
      - 16*a22 *a33 *b33 *r4*u1*v3 - 16*a22 *a33 *b33 *r4*v1*v3

              4    3             4         4    3             2
      + 16*a22 *a33 *b33*k1*n1*u1  + 32*a22 *a33 *b33*k1*n1*u1 *u2*u3

              4    3             2         4    3                   2
      + 80*a22 *a33 *b33*k1*n1*u1  + 16*a22 *a33 *b33*k1*n1*u1*u2*v2

              4    3                2         4    3
      - 64*a22 *a33 *b33*k1*n1*u1*u3  + 16*a22 *a33 *b33*k1*n1*u1*v1

              4    3              2         4    3                    2
      + 16*a22 *a33 *b33*k12*n1*u1  + 16*a22 *a33 *b33*k12*n1*u1*v1*v2

              4    3                          4    3              2   2
      + 16*a22 *a33 *b33*k12*n1*u1*v1 + 16*a22 *a33 *b33*k12*n1*u2 *u3

              4    3              3            4    3              3
      + 16*a22 *a33 *b33*k16*n1*u1 *v1 + 16*a22 *a33 *b33*k16*n1*u1 *v3

              4    3                                4    3                 3
      + 16*a22 *a33 *b33*k16*n1*u1*u2*u3*v1 + 16*a22 *a33 *b33*k16*n1*u1*v2

              4    3              3            4    3                 3
      + 16*a22 *a33 *b33*k16*n1*u2 *v1 + 16*a22 *a33 *b33*k16*n1*v1*v2

              4    3                 2            4    3
      + 16*a22 *a33 *b33*k38*n1*u1*u2 *v2 - 32*a22 *a33 *b33*k38*n1*u1*u3*v1

              4    3                          4    3              3
      - 16*a22 *a33 *b33*k38*n1*u1*u3 + 16*a22 *a33 *b33*k38*n1*u2 *v3

              4    3                                4    3
      + 16*a22 *a33 *b33*k38*n1*u2*u3*v2*v3 + 32*a22 *a33 *b33*k38*n1*u2*v1*v2

              4    3              2         4    3
      - 16*a22 *a33 *b33*k38*n1*v1  + 16*a22 *a33 *b33*n1*r4*v1

              4    2    2                        4    2    2
      - 32*a22 *a33 *b33 *k1*n1*u1*u3*v1 - 16*a22 *a33 *b33 *k1*n1*u1*u3

              4    2    2                        4    2    2         4
      + 32*a22 *a33 *b33 *k1*n1*u2*v1*v2 - 32*a22 *a33 *b33 *k1*n1*u3

              4    2    2         2   2         4    2    2         2
      - 32*a22 *a33 *b33 *k1*n1*v1 *v2  - 16*a22 *a33 *b33 *k1*n1*v1

              4    2    2                      4    2    2
      + 16*a22 *a33 *b33 *k12*n1*u1*u3 - 16*a22 *a33 *b33 *k12*n1*u2*u3*v2*v3

              4    2    2          2         4    2    2
      + 16*a22 *a33 *b33 *k12*n1*v1  + 16*a22 *a33 *b33 *k16*n1*u1*u2*v1*v2

              4    2    2             2            4    2    2          2
      + 16*a22 *a33 *b33 *k16*n1*u1*v1 *v2 + 16*a22 *a33 *b33 *k16*n1*u2 *v1*v3

              4    2    2
      + 16*a22 *a33 *b33 *k16*n1*u2*u3*v1*v2

              4    2    2             2            4    2    2          3
      + 16*a22 *a33 *b33 *k16*n1*u2*v1 *v3 + 16*a22 *a33 *b33 *k16*n1*v1 *v3

              4    2          2   3         4    2          2
      - 32*a22 *a33 *b33*k1*n1 *u1  - 96*a22 *a33 *b33*k1*n1 *u1

              4    2          2   2            4    2          2
      - 32*a22 *a33 *b33*k1*n1 *u2 *u3 + 32*a22 *a33 *b33*k1*n1 *u2*v1

              4    2           2   3         4    2           2
      + 16*a22 *a33 *b33*k12*n1 *u1  - 16*a22 *a33 *b33*k12*n1 *u1*u2*v1

              4    2           2   2            4    2           2
      + 16*a22 *a33 *b33*k12*n1 *u2 *u3 - 32*a22 *a33 *b33*k12*n1 *u2*v1

              4    2           2   2            4    2           2   2
      - 16*a22 *a33 *b33*k16*n1 *u1 *v2 - 16*a22 *a33 *b33*k16*n1 *u3 *v1

              4    2           2      2         4    2           2   3
      - 16*a22 *a33 *b33*k16*n1 *u3*v1  + 16*a22 *a33 *b33*k38*n1 *u2

              4    2           2                  4    2           2   3
      + 16*a22 *a33 *b33*k38*n1 *u2*u3*v1 + 16*a22 *a33 *b33*k38*n1 *u3

             4        3         3           4        3
      - 8*a22 *a33*b33 *k1*n1*u1 *v3 - 8*a22 *a33*b33 *k1*n1*u1*u2*u3*v1

              4        3            2            4        3            3
      + 32*a22 *a33*b33 *k1*n1*u1*u3 *v2 + 32*a22 *a33*b33 *k1*n1*u1*v1

              4        3                        4        3            3
      - 16*a22 *a33*b33 *k1*n1*u1*v1*v3 - 16*a22 *a33*b33 *k1*n1*u1*v2

              4        3                     4        3         3
      + 12*a22 *a33*b33 *k1*n1*u1*v3 - 16*a22 *a33*b33 *k1*n1*u2 *v1

              4        3         3            4        3
      + 32*a22 *a33*b33 *k1*n1*u3 *v3 + 12*a22 *a33*b33 *k1*n1*v1*v3

             4        3          3           4        3
      + 4*a22 *a33*b33 *k12*n1*u1 *v3 + 4*a22 *a33*b33 *k12*n1*u1*u2*u3*v1

             4        3             2           4        3             3
      - 8*a22 *a33*b33 *k12*n1*u1*u3 *v2 - 8*a22 *a33*b33 *k12*n1*u1*v1

             4        3                        4        3             3
      + 4*a22 *a33*b33 *k12*n1*u1*v1*v3 + 4*a22 *a33*b33 *k12*n1*u1*v2

             4        3                     4        3          3
      - 4*a22 *a33*b33 *k12*n1*u1*v3 + 4*a22 *a33*b33 *k12*n1*u2 *v1

             4        3          3           4        3             2
      - 8*a22 *a33*b33 *k12*n1*u3 *v3 - 4*a22 *a33*b33 *k12*n1*u3*v2

             4        3          3        4        3
      + 4*a22 *a33*b33 *k12*n1*v1  - 4*a22 *a33*b33 *k12*n1*v1*v3

             4        3                2        4        3                2
      - 4*a22 *a33*b33 *k16*n1*u1*v1*v3  - 4*a22 *a33*b33 *k16*n1*u2*v2*v3

             4        3             3        4        3          2   2
      - 4*a22 *a33*b33 *k16*n1*u2*v3  - 8*a22 *a33*b33 *k16*n1*v1 *v3

             4        3             3         4              3   2
      - 4*a22 *a33*b33 *k16*n1*v1*v3  + 16*a22 *a33*b33*k1*n1 *u1

              4               3   2        4    4
      - 16*a22 *a33*b33*k12*n1 *u1  - 8*a22 *b33 *k1*n1*u1*u2*v1*v2

             4    4               2        4    4            2
      + 8*a22 *b33 *k1*n1*u1*u3*v3  - 8*a22 *b33 *k1*n1*u1*v1 *v2

             4    4         2              4    4
      - 8*a22 *b33 *k1*n1*u2 *v1*v3 - 8*a22 *b33 *k1*n1*u2*u3*v1*v2

             4    4            2           4    4         3
      - 8*a22 *b33 *k1*n1*u2*v1 *v3 - 8*a22 *b33 *k1*n1*v1 *v3

             4    4            2           4    4
      + 8*a22 *b33 *k1*n1*v1*v2 *v3 + 4*a22 *b33 *k12*n1*u1*u2*v1*v2

             4    4             2           4    4          2
      + 4*a22 *b33 *k12*n1*u1*v1 *v2 + 4*a22 *b33 *k12*n1*u2 *v1*v3

             4    4                           4    4             2
      + 4*a22 *b33 *k12*n1*u2*u3*v1*v2 + 4*a22 *b33 *k12*n1*u2*v1 *v3

             4    4          3           4    3      2   2
      + 4*a22 *b33 *k12*n1*v1 *v3 + 8*a22 *b33 *k1*n1 *u1 *v2

              4    3      2                  4    3      2   2
      - 16*a22 *b33 *k1*n1 *u2*v1*v3 + 16*a22 *b33 *k1*n1 *u3 *v1

              4    3      2      2        4    3       2
      + 16*a22 *b33 *k1*n1 *u3*v1  - 8*a22 *b33 *k12*n1 *u2*v1*v3

             4    3       2   2           4    3       2      2
      - 4*a22 *b33 *k12*n1 *u3 *v1 - 4*a22 *b33 *k12*n1 *u3*v1

             4    3       2               4    2      3
      + 8*a22 *b33 *k12*n1 *v2*v3 - 16*a22 *b33 *k1*n1 *u1*u3

             3    4    2                 3    4    2
      + 8*a22 *a33 *b33 *r4*u1*v3 + 8*a22 *a33 *b33 *r4*v1*v3

              3    3    2                        3    3    2
      + 32*a22 *a33 *b33 *k1*n1*u1*u3*v1 + 16*a22 *a33 *b33 *k1*n1*u1*u3

              3    3    2                        3    3    2         4
      - 32*a22 *a33 *b33 *k1*n1*u2*v1*v2 + 32*a22 *a33 *b33 *k1*n1*u3

              3    3    2         2   2         3    3    2         2
      + 32*a22 *a33 *b33 *k1*n1*v1 *v2  + 16*a22 *a33 *b33 *k1*n1*v1

              3    3    2                      3    3    2
      - 16*a22 *a33 *b33 *k12*n1*u1*u3 + 16*a22 *a33 *b33 *k12*n1*u2*u3*v2*v3

              3    3    2          2         3    3    2
      - 16*a22 *a33 *b33 *k12*n1*v1  - 16*a22 *a33 *b33 *k16*n1*u1*u2*v1*v2

              3    3    2             2            3    3    2          2
      - 16*a22 *a33 *b33 *k16*n1*u1*v1 *v2 - 16*a22 *a33 *b33 *k16*n1*u2 *v1*v3

              3    3    2
      - 16*a22 *a33 *b33 *k16*n1*u2*u3*v1*v2

              3    3    2             2            3    3    2          3
      - 16*a22 *a33 *b33 *k16*n1*u2*v1 *v3 - 16*a22 *a33 *b33 *k16*n1*v1 *v3

             3    2    3         3           3    2    3
      - 8*a22 *a33 *b33 *k1*n1*u1 *v3 - 8*a22 *a33 *b33 *k1*n1*u1*u2*u3*v1

              3    2    3            2            3    2    3            3
      + 32*a22 *a33 *b33 *k1*n1*u1*u3 *v2 + 32*a22 *a33 *b33 *k1*n1*u1*v1

              3    2    3                        3    2    3            3
      - 16*a22 *a33 *b33 *k1*n1*u1*v1*v3 - 16*a22 *a33 *b33 *k1*n1*u1*v2

              3    2    3                     3    2    3         3
      + 12*a22 *a33 *b33 *k1*n1*u1*v3 - 16*a22 *a33 *b33 *k1*n1*u2 *v1

              3    2    3         3            3    2    3
      + 32*a22 *a33 *b33 *k1*n1*u3 *v3 + 12*a22 *a33 *b33 *k1*n1*v1*v3

             3    2    3          3           3    2    3
      + 4*a22 *a33 *b33 *k12*n1*u1 *v3 + 4*a22 *a33 *b33 *k12*n1*u1*u2*u3*v1

             3    2    3             2           3    2    3             3
      - 8*a22 *a33 *b33 *k12*n1*u1*u3 *v2 - 8*a22 *a33 *b33 *k12*n1*u1*v1

             3    2    3                        3    2    3             3
      + 4*a22 *a33 *b33 *k12*n1*u1*v1*v3 + 4*a22 *a33 *b33 *k12*n1*u1*v2

             3    2    3                     3    2    3          3
      - 4*a22 *a33 *b33 *k12*n1*u1*v3 + 4*a22 *a33 *b33 *k12*n1*u2 *v1

             3    2    3          3           3    2    3             2
      - 8*a22 *a33 *b33 *k12*n1*u3 *v3 - 4*a22 *a33 *b33 *k12*n1*u3*v2

             3    2    3          3        3    2    3
      + 4*a22 *a33 *b33 *k12*n1*v1  - 4*a22 *a33 *b33 *k12*n1*v1*v3

             3    2    3                2        3    2    3                2
      - 4*a22 *a33 *b33 *k16*n1*u1*v1*v3  - 4*a22 *a33 *b33 *k16*n1*u2*v2*v3

             3    2    3             3        3    2    3          2   2
      - 4*a22 *a33 *b33 *k16*n1*u2*v3  - 8*a22 *a33 *b33 *k16*n1*v1 *v3

             3    2    3             3         3    2    2      2
      - 4*a22 *a33 *b33 *k16*n1*v1*v3  + 16*a22 *a33 *b33 *k1*n1 *u1*u3*v2

              3    2    2      2   3         3    2    2      2
      + 16*a22 *a33 *b33 *k1*n1 *u2  - 32*a22 *a33 *b33 *k1*n1 *u2*u3*v1

              3    2    2      2              3    2    2       2
      - 16*a22 *a33 *b33 *k1*n1 *v1*v2 + 8*a22 *a33 *b33 *k12*n1 *u1*u3*v2

              3    2    2       2              3    2    2       2   3
      + 16*a22 *a33 *b33 *k12*n1 *u1*v2 + 8*a22 *a33 *b33 *k12*n1 *u2

             3    2    2       2   2            3    2    2       2
      + 8*a22 *a33 *b33 *k12*n1 *u2 *v3 + 16*a22 *a33 *b33 *k12*n1 *u2*u3*v1

             3    2    2       2                 3    2    2       2      2
      + 8*a22 *a33 *b33 *k16*n1 *u2*v2*v3 + 8*a22 *a33 *b33 *k16*n1 *u2*v3

             3    2    2       2   2           3    2    2       2
      + 8*a22 *a33 *b33 *k16*n1 *v2 *v3 + 8*a22 *a33 *b33 *k38*n1 *u2*v1*v3

             3    2    2       2      2        3    2    2       2   2
      + 8*a22 *a33 *b33 *k38*n1 *u2*v2  + 8*a22 *a33 *b33 *k38*n1 *u3 *v1

             3        4                          3        4               2
      + 8*a22 *a33*b33 *k1*n1*u1*u2*v1*v2 - 8*a22 *a33*b33 *k1*n1*u1*u3*v3

             3        4            2           3        4         2
      + 8*a22 *a33*b33 *k1*n1*u1*v1 *v2 + 8*a22 *a33*b33 *k1*n1*u2 *v1*v3

             3        4                          3        4            2
      + 8*a22 *a33*b33 *k1*n1*u2*u3*v1*v2 + 8*a22 *a33*b33 *k1*n1*u2*v1 *v3

             3        4         3           3        4            2
      + 8*a22 *a33*b33 *k1*n1*v1 *v3 - 8*a22 *a33*b33 *k1*n1*v1*v2 *v3

             3        4                           3        4             2
      - 4*a22 *a33*b33 *k12*n1*u1*u2*v1*v2 - 4*a22 *a33*b33 *k12*n1*u1*v1 *v2

             3        4          2              3        4
      - 4*a22 *a33*b33 *k12*n1*u2 *v1*v3 - 4*a22 *a33*b33 *k12*n1*u2*u3*v1*v2

             3        4             2           3        4          3
      - 4*a22 *a33*b33 *k12*n1*u2*v1 *v3 - 4*a22 *a33*b33 *k12*n1*v1 *v3

              3        2      3              3    5               2
      - 16*a22 *a33*b33 *k1*n1 *u1*u3 + 3*a22 *b33 *k1*n1*u1*v1*v3

             3    5               2        3    5            3
      + 2*a22 *b33 *k1*n1*u2*v2*v3  + 3*a22 *b33 *k1*n1*u2*v3

             3    5         2   2        3    5            3
      + 6*a22 *b33 *k1*n1*v1 *v3  + 2*a22 *b33 *k1*n1*v1*v3

           3    5                2      3    5                2
      - a22 *b33 *k12*n1*u1*v1*v3  - a22 *b33 *k12*n1*u2*v2*v3

           3    5             3        3    5          2   2
      - a22 *b33 *k12*n1*u2*v3  - 2*a22 *b33 *k12*n1*v1 *v3

           3    5             3        3    4      2
      - a22 *b33 *k12*n1*v1*v3  - 8*a22 *b33 *k1*n1 *u2*v2*v3

             3    4      2      2        3    4      2   2
      - 4*a22 *b33 *k1*n1 *u2*v3  - 8*a22 *b33 *k1*n1 *v2 *v3

             3    4       2                 3    4       2      2
      + 2*a22 *b33 *k12*n1 *u2*v2*v3 + 4*a22 *b33 *k12*n1 *u2*v3

             3    4       2   2           3    3      3
      + 2*a22 *b33 *k12*n1 *v2 *v3 + 4*a22 *b33 *k1*n1 *u1*v3

             3    3       3              2    3    3         3
      + 4*a22 *b33 *k12*n1 *u1*v3 + 8*a22 *a33 *b33 *k1*n1*u1 *v3

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

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

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

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

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

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

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

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

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

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

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

             2    3    3                2        2    3    3             3
      + 4*a22 *a33 *b33 *k16*n1*u2*v2*v3  + 4*a22 *a33 *b33 *k16*n1*u2*v3

             2    3    3          2   2        2    3    3             3
      + 8*a22 *a33 *b33 *k16*n1*v1 *v3  + 4*a22 *a33 *b33 *k16*n1*v1*v3

             2    2    4                          2    2    4               2
      + 8*a22 *a33 *b33 *k1*n1*u1*u2*v1*v2 - 8*a22 *a33 *b33 *k1*n1*u1*u3*v3

             2    2    4            2           2    2    4         2
      + 8*a22 *a33 *b33 *k1*n1*u1*v1 *v2 + 8*a22 *a33 *b33 *k1*n1*u2 *v1*v3

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

             2    2    4         3           2    2    4            2
      + 8*a22 *a33 *b33 *k1*n1*v1 *v3 - 8*a22 *a33 *b33 *k1*n1*v1*v2 *v3

             2    2    4                           2    2    4             2
      - 4*a22 *a33 *b33 *k12*n1*u1*u2*v1*v2 - 4*a22 *a33 *b33 *k12*n1*u1*v1 *v2

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

             2    2    4             2           2    2    4          3
      - 4*a22 *a33 *b33 *k12*n1*u2*v1 *v3 - 4*a22 *a33 *b33 *k12*n1*v1 *v3

             2    2    3      2   2            2    2    3      2
      - 8*a22 *a33 *b33 *k1*n1 *u1 *v2 + 16*a22 *a33 *b33 *k1*n1 *u2*v1*v3

              2    2    3      2   2            2    2    3      2      2
      - 16*a22 *a33 *b33 *k1*n1 *u3 *v1 - 16*a22 *a33 *b33 *k1*n1 *u3*v1

             2    2    3       2                 2    2    3       2   2
      + 8*a22 *a33 *b33 *k12*n1 *u2*v1*v3 + 4*a22 *a33 *b33 *k12*n1 *u3 *v1

             2    2    3       2      2        2    2    3       2
      + 4*a22 *a33 *b33 *k12*n1 *u3*v1  - 8*a22 *a33 *b33 *k12*n1 *v2*v3

             2        5               2        2        5               2
      - 3*a22 *a33*b33 *k1*n1*u1*v1*v3  - 2*a22 *a33*b33 *k1*n1*u2*v2*v3

             2        5            3        2        5         2   2
      - 3*a22 *a33*b33 *k1*n1*u2*v3  - 6*a22 *a33*b33 *k1*n1*v1 *v3

             2        5            3      2        5                2
      - 2*a22 *a33*b33 *k1*n1*v1*v3  + a22 *a33*b33 *k12*n1*u1*v1*v3

           2        5                2      2        5             3
      + a22 *a33*b33 *k12*n1*u2*v2*v3  + a22 *a33*b33 *k12*n1*u2*v3

             2        5          2   2      2        5             3
      + 2*a22 *a33*b33 *k12*n1*v1 *v3  + a22 *a33*b33 *k12*n1*v1*v3

             2        3      3              2        3       3
      + 4*a22 *a33*b33 *k1*n1 *u1*v3 + 4*a22 *a33*b33 *k12*n1 *u1*v3

                 3    4                              3    4               2
      - 8*a22*a33 *b33 *k1*n1*u1*u2*v1*v2 + 8*a22*a33 *b33 *k1*n1*u1*u3*v3

                 3    4            2               3    4         2
      - 8*a22*a33 *b33 *k1*n1*u1*v1 *v2 - 8*a22*a33 *b33 *k1*n1*u2 *v1*v3

                 3    4                              3    4            2
      - 8*a22*a33 *b33 *k1*n1*u2*u3*v1*v2 - 8*a22*a33 *b33 *k1*n1*u2*v1 *v3

                 3    4         3               3    4            2
      - 8*a22*a33 *b33 *k1*n1*v1 *v3 + 8*a22*a33 *b33 *k1*n1*v1*v2 *v3

                 3    4                               3    4             2
      + 4*a22*a33 *b33 *k12*n1*u1*u2*v1*v2 + 4*a22*a33 *b33 *k12*n1*u1*v1 *v2

                 3    4          2                  3    4
      + 4*a22*a33 *b33 *k12*n1*u2 *v1*v3 + 4*a22*a33 *b33 *k12*n1*u2*u3*v1*v2

                 3    4             2               3    4          3
      + 4*a22*a33 *b33 *k12*n1*u2*v1 *v3 + 4*a22*a33 *b33 *k12*n1*v1 *v3

                 2    5               2            2    5               2
      - 3*a22*a33 *b33 *k1*n1*u1*v1*v3  - 2*a22*a33 *b33 *k1*n1*u2*v2*v3

                 2    5            3            2    5         2   2
      - 3*a22*a33 *b33 *k1*n1*u2*v3  - 6*a22*a33 *b33 *k1*n1*v1 *v3

                 2    5            3          2    5                2
      - 2*a22*a33 *b33 *k1*n1*v1*v3  + a22*a33 *b33 *k12*n1*u1*v1*v3

               2    5                2          2    5             3
      + a22*a33 *b33 *k12*n1*u2*v2*v3  + a22*a33 *b33 *k12*n1*u2*v3

                 2    5          2   2          2    5             3
      + 2*a22*a33 *b33 *k12*n1*v1 *v3  + a22*a33 *b33 *k12*n1*v1*v3

                 2    4      2                     2    4      2      2
      + 8*a22*a33 *b33 *k1*n1 *u2*v2*v3 + 4*a22*a33 *b33 *k1*n1 *u2*v3

                 2    4      2   2               2    4       2
      + 8*a22*a33 *b33 *k1*n1 *v2 *v3 - 2*a22*a33 *b33 *k12*n1 *u2*v2*v3

                 2    4       2      2            2    4       2   2
      - 4*a22*a33 *b33 *k12*n1 *u2*v3  - 2*a22*a33 *b33 *k12*n1 *v2 *v3

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

             3    5            3        3    5         2   2
      + 3*a33 *b33 *k1*n1*u2*v3  + 6*a33 *b33 *k1*n1*v1 *v3

             3    5            3      3    5                2
      + 2*a33 *b33 *k1*n1*v1*v3  - a33 *b33 *k12*n1*u1*v1*v3

           3    5                2      3    5             3
      - a33 *b33 *k12*n1*u2*v2*v3  - a33 *b33 *k12*n1*u2*v3

             3    5          2   2      3    5             3         4
      - 2*a33 *b33 *k12*n1*v1 *v3  - a33 *b33 *k12*n1*v1*v3 )/(16*a22 *b33*n1

            3      2              2      3
       *(a22  - a22 *a33 - a22*a33  + 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*a33**2*u3**2 -
 4*a22*a33*b33*u3*v3 - 4*a22*a33*n1*u1 - a22*b33**2*v3**2 - 2*a22*b33*n1*v1 + 
a33*b33**2*v3**2 - 2*a33*b33*n1*v1)/(4*a22*(a22 - a33))$

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