Solution 28 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


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


m2=0


     - b33*n3
m3=-----------
    a22 - a33


         2                2
r6=(4*a22 *b33*k1*n3 - a22 *b33*n3*q1 - 4*a22*a33*b33*k1*n3 + a22*a33*b33*n3*q1

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


        1
     - ---*b33*k1*n1*n3
        4
r5=---------------------
         3      2
      a22  - a22 *a33


     11     2              1     2              5
r4=(----*a22 *b33*k1*n1 - ---*a22 *b33*n1*q1 - ---*a22*a33*b33*k1*n1
     4                     2                    2

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

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


       3              2                    2          1           2
r3=(a22 *n3*q1 - 2*a22 *a33*n3*q1 + a22*a33 *n3*q1 - ---*a22*k1*n1 *n3
                                                      2

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

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


      1
     ---*k1*n1*n3
      2
r2=----------------
       2
    a22  - a22*a33


     1     3            3          3     2                  2
r1=(---*a22 *k1*n1 + a22 *n1*q1 - ---*a22 *a33*k1*n1 - 2*a22 *a33*n1*q1
     2                             2

        3         2                2          1           3    1              2
     + ---*a22*a33 *k1*n1 + a22*a33 *n1*q1 - ---*a22*k1*n1  + ---*a22*k1*n1*n3
        2                                     2                2

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

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


      1     2    2          1     2    2
     ---*a22 *b33 *k1*n1 - ---*a33 *b33 *k1*n1
      4                     4
q20=-------------------------------------------
              5        4          3    2
           a22  - 2*a22 *a33 + a22 *a33


             3    2       1     3    2           2        2
q19=( - 2*a22 *b33 *k1 + ---*a22 *b33 *q1 + 4*a22 *a33*b33 *k1
                          4

         1     2        2               2    2       1         2    2
      - ---*a22 *a33*b33 *q1 - 2*a22*a33 *b33 *k1 + ---*a22*a33 *b33 *q1
         2                                           4

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

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


      5     2    2                       2             1     2    2
     ---*a22 *b33 *k1*n1*n3 - a22*a33*b33 *k1*n1*n3 - ---*a33 *b33 *k1*n1*n3
      4                                                4
q18=-------------------------------------------------------------------------
                      6        5            4    2      3    3
                   a22  - 3*a22 *a33 + 3*a22 *a33  - a22 *a33


q17=0


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

         1     2        2               2    2       1         2    2
      - ---*a22 *a33*b33 *q1 - 2*a22*a33 *b33 *k1 + ---*a22*a33 *b33 *q1
         2                                           4

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

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


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


            - b33*k1*n1*n3
q13=------------------------------
        3        2              2
     a22  - 2*a22 *a33 + a22*a33


         1                   1
      - ---*a22*b33*k1*n1 - ---*a33*b33*k1*n1
         2                   2
q11=------------------------------------------
                    3      2
                 a22  - a22 *a33


           4            3                3             3             2    2
q10=( - a22 *k38 + 2*a22 *a33*k38 + 6*a22 *b33*k1 - a22 *b33*q1 - a22 *a33 *k38

              2                   2                       2
      - 12*a22 *a33*b33*k1 + 2*a22 *a33*b33*q1 + 6*a22*a33 *b33*k1

               2           1               2    3               2
      - a22*a33 *b33*q1 + ---*a22*b33*k1*n1  + ---*a22*b33*k1*n3
                           2                    2

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


          2
     - a22 *k38*n3 + a22*a33*k38*n3 + 2*a22*b33*k1*n3 - 2*a33*b33*k1*n3
q9=---------------------------------------------------------------------
                          3        2              2
                       a22  - 2*a22 *a33 + a22*a33


     - k1*n1
q8=----------
      a22


            3         3           2               2                   2
q7=( - 4*a22 *k1 + a22 *q1 + 8*a22 *a33*k1 - 2*a22 *a33*q1 - 4*a22*a33 *k1

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

             2
    + a22*a33 )


        3                      1
     - ---*a22*b33*k1*n1*n3 + ---*a33*b33*k1*n1*n3
        2                      2
q6=------------------------------------------------
               4        3          2    2
            a22  - 2*a22 *a33 + a22 *a33


       2
    a22 *k38*n3 - a22*a33*k38*n3 - 2*a22*b33*k1*n3 + 2*a33*b33*k1*n3
q5=------------------------------------------------------------------
                         3        2              2
                      a22  - 2*a22 *a33 + a22*a33


          4            3                3             3             2    2
q4=( - a22 *k38 + 2*a22 *a33*k38 + 6*a22 *b33*k1 - a22 *b33*q1 - a22 *a33 *k38

             2                   2                       2
     - 12*a22 *a33*b33*k1 + 2*a22 *a33*b33*q1 + 6*a22*a33 *b33*k1

              2           1               2    3               2
     - a22*a33 *b33*q1 - ---*a22*b33*k1*n1  + ---*a22*b33*k1*n3
                          2                    2

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


       k1*n1*n3
q3=----------------
       2
    a22  - a22*a33


q2=0


           3               3         3          1         3
      - a22 *b33*k16*n3 - ---*a22*b33 *k1*n3 + ---*a33*b33 *k1*n3
                           8                    8
p56=--------------------------------------------------------------
                              5      4
                           a22  - a22 *a33


p55=0


           3               7         3          1         3
      - a22 *b33*k16*n3 - ---*a22*b33 *k1*n3 + ---*a33*b33 *k1*n3
                           8                    8
p54=--------------------------------------------------------------
                              5      4
                           a22  - a22 *a33


p53=0


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

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

         4    2
    + a22 *a33 )


p51=0


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

         1             3          1      2    3            6        5
      + ---*a22*a33*b33 *k1*n1 - ----*a33 *b33 *k1*n1)/(a22  - 2*a22 *a33
         2                        16

         4    2
    + a22 *a33 )


           3               7         3          1         3
      - a22 *b33*k16*n3 - ---*a22*b33 *k1*n3 + ---*a33*b33 *k1*n3
                           8                    8
p49=--------------------------------------------------------------
                              5      4
                           a22  - a22 *a33


p48=0


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

         1             3          1      2    3            6        5
      + ---*a22*a33*b33 *k1*n1 - ----*a33 *b33 *k1*n1)/(a22  - 2*a22 *a33
         2                        16

         4    2
    + a22 *a33 )


        4               3                 2    2           3     2    2
p46=(a22 *k16*n3 - 2*a22 *a33*k16*n3 + a22 *a33 *k16*n3 - ---*a22 *b33 *k1*n3
                                                           2

                     2          3     2    2            5        4
      + 3*a22*a33*b33 *k1*n3 - ---*a33 *b33 *k1*n3)/(a22  - 3*a22 *a33
                                2

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


p45=0


        5               4                 3    2           7     3    2
p44=(a22 *k16*n3 - 2*a22 *a33*k16*n3 + a22 *a33 *k16*n3 + ---*a22 *b33 *k1*n3
                                                           8

         15     2        2          9         2    2          1     3    2
      - ----*a22 *a33*b33 *k1*n3 + ---*a22*a33 *b33 *k1*n3 - ---*a33 *b33 *k1*n3
         8                          8                         8

           6        5            4    2      3    3
     )/(a22  - 3*a22 *a33 + 3*a22 *a33  - a22 *a33 )


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

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

   )


      7         2       1         2
     ---*a22*b33 *k1 - ---*a33*b33 *k1
      8                 8
p42=-----------------------------------
                      3
                   a22


        5               4                 3    2           7     3    2
p41=(a22 *k16*n3 - 2*a22 *a33*k16*n3 + a22 *a33 *k16*n3 + ---*a22 *b33 *k1*n3
                                                           8

         15     2        2          9         2    2          1     3    2
      - ----*a22 *a33*b33 *k1*n3 + ---*a22*a33 *b33 *k1*n3 - ---*a33 *b33 *k1*n3
         8                          8                         8

           6        5            4    2      3    3
     )/(a22  - 3*a22 *a33 + 3*a22 *a33  - a22 *a33 )


           3             2                 2
p40=( - a22 *k38*n3 + a22 *a33*k38*n3 + a22 *b33*k1*n3 - 2*a22*a33*b33*k1*n3

           2                4        3          2    2
      + a33 *b33*k1*n3)/(a22  - 2*a22 *a33 + a22 *a33 )


p39=0


         1                   1
      - ---*a22*b33*k1*n1 - ---*a33*b33*k1*n1
         4                   4
p38=------------------------------------------
                    3      2
                 a22  - a22 *a33


      1               1
     ---*a22*k1*n3 + ---*a33*k1*n3
      2               2
p37=-------------------------------
               2
            a22  - a22*a33


p36=0


p35=0


p34=0


      1         2       1         2
     ---*a22*b33 *k1 + ---*a33*b33 *k1
      8                 8
p33=-----------------------------------
                      3
                   a22


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


p31=0


p30=0


      - k38*n3
p29=-----------
     a22 - a33


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


p27=0


      - 2*b33*k1*n3
p26=----------------
        2
     a22  - a22*a33


p25=0


         5                   1
      - ---*a22*b33*k1*n1 - ---*a33*b33*k1*n1
         4                   4
p24=------------------------------------------
                    3      2
                 a22  - a22 *a33


      3               1
     ---*a22*k1*n3 - ---*a33*k1*n3
      2               2
p23=-------------------------------
               2
            a22  - a22*a33


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       1         2
      - ---*a22*b33 *k1 - ---*a33*b33 *k1
         8                 8
p20=--------------------------------------
                        3
                     a22


        4             3               7     2    2                    2
p19=(a22 *k16*n1 - a22 *a33*k16*n1 + ---*a22 *b33 *k1*n1 - a22*a33*b33 *k1*n1
                                      8

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


p18=0


p17=0


p16

     3           1     2               7         2          1         2
  a22 *k16*n1 - ---*a22 *b33*k38*n1 + ---*a22*b33 *k1*n1 - ---*a33*b33 *k1*n1
                 2                     8                    8
=-----------------------------------------------------------------------------
                                   4      3
                                a22  - a22 *a33


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

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


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


      - k38*n3
p13=-----------
     a22 - a33


      1               1
     ---*a22*k1*n1 + ---*a33*k1*n1
      2               2
p12=-------------------------------
               2
            a22  - a22*a33


p11=0


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


p9=0


p8= - 4*k1


     3               1
    ---*a22*k1*n1 - ---*a33*k1*n1
     2               2
p7=-------------------------------
              2
           a22  - a22*a33


     - 2*b33*k1*n3
p6=----------------
       2
    a22  - a22*a33


p5=0


          3             2               5     2
p4=( - a22 *k38*n1 + a22 *a33*k38*n1 - ---*a22 *b33*k1*n1 + a22*a33*b33*k1*n1
                                        4

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


     3               1
    ---*a22*k1*n3 - ---*a33*k1*n3
     2               2
p3=-------------------------------
              2
           a22  - a22*a33


p2=0


     3               1
    ---*a22*k1*n1 - ---*a33*k1*n1
     2               2
p1=-------------------------------
              2
           a22  - a22*a33


k125=0


       1     3    2        3          4       1          4
      ---*a22 *b33 *k16 + ----*a22*b33 *k1 - ----*a33*b33 *k1
       4                   32                 32
k124=---------------------------------------------------------
                                  5
                               a22


k123=0


       1     3    2        5          4       1          4
      ---*a22 *b33 *k16 + ----*a22*b33 *k1 - ----*a33*b33 *k1
       4                   32                 32
k122=---------------------------------------------------------
                                  5
                               a22


k121=0


k120=0


k119=0


k118=0


       1     3    2        3          4       1          4
      ---*a22 *b33 *k16 + ----*a22*b33 *k1 - ----*a33*b33 *k1
       4                   32                 32
k117=---------------------------------------------------------
                                  5
                               a22


k116=0


       1     3    2        5          4       1          4
      ---*a22 *b33 *k16 + ----*a22*b33 *k1 - ----*a33*b33 *k1
       2                   16                 16
k115=---------------------------------------------------------
                                  5
                               a22


k114=0


k113=0


       1     3    2        5          4       1          4
      ---*a22 *b33 *k16 + ----*a22*b33 *k1 - ----*a33*b33 *k1
       4                   32                 32
k112=---------------------------------------------------------
                                  5
                               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


        2        1                1
     a22 *k38 + ---*a22*b33*k1 + ---*a33*b33*k1
                 2                2
k95=--------------------------------------------
                           2
                        a22


k94=0


k93=0


k91=0


           3            3         3       1         3
      - a22 *b33*k16 - ---*a22*b33 *k1 + ---*a33*b33 *k1
                        8                 8
k90=-----------------------------------------------------
                               4
                            a22


k89=0


           3            3         3       1         3
      - a22 *b33*k16 - ---*a22*b33 *k1 + ---*a33*b33 *k1
                        8                 8
k88=-----------------------------------------------------
                               4
                            a22


k87=0


k86=0


k85=0


k84=0


           3            3         3       1         3
      - a22 *b33*k16 - ---*a22*b33 *k1 + ---*a33*b33 *k1
                        8                 8
k83=-----------------------------------------------------
                               4
                            a22


k82=0


k81=0


         7         2       1         2
      - ---*a22*b33 *k1 + ---*a33*b33 *k1
         4                 4
k80=--------------------------------------
                        3
                     a22


k79=0


k78=0


k77=0


k76=0


k75=0


k74=k38


k73=0


k72=0


        3        3         2       1         2
     a22 *k16 + ---*a22*b33 *k1 - ---*a33*b33 *k1
                 8                 8
k71=----------------------------------------------
                            3
                         a22


k70=0


k69=k16


k68=0


k67=0


        3        7         2       1         2
     a22 *k16 + ---*a22*b33 *k1 - ---*a33*b33 *k1
                 8                 8
k66=----------------------------------------------
                            3
                         a22


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


k64=0


k63=0


      1            1
     ---*a22*k1 + ---*a33*k1
      2            2
k62=-------------------------
               a22


k61=0


k59=0


k58=0


k57=k1


k56=0


k55=0


k54=0


k53=0


           3            3         3       1         3
      - a22 *b33*k16 - ---*a22*b33 *k1 + ---*a33*b33 *k1
                        8                 8
k52=-----------------------------------------------------
                               4
                            a22


k51=0


           3            3         3       1         3
      - a22 *b33*k16 - ---*a22*b33 *k1 + ---*a33*b33 *k1
                        8                 8
k50=-----------------------------------------------------
                               4
                            a22


k49=0


k48=0


           3            3         3       1         3
      - a22 *b33*k16 - ---*a22*b33 *k1 + ---*a33*b33 *k1
                        8                 8
k47=-----------------------------------------------------
                               4
                            a22


k46=0


k45=0


k44=0


         7         2       1         2
      - ---*a22*b33 *k1 + ---*a33*b33 *k1
         4                 4
k43=--------------------------------------
                        3
                     a22


k42=0


k41=0


k40=0


k39=0


k37=0


k36=0


k35=0


k34=0


k33=0


         7         2       1         2
      - ---*a22*b33 *k1 + ---*a33*b33 *k1
         4                 4
k32=--------------------------------------
                        3
                     a22


k31=0


k30=0


k29=0


k28=0


k27=0


k26=0


k25=0


k24=0


k23=0


k22=0


        3        3         2       1         2
     a22 *k16 + ---*a22*b33 *k1 - ---*a33*b33 *k1
                 8                 8
k21=----------------------------------------------
                            3
                         a22


k20=0


        3        7         2       1         2
     a22 *k16 + ---*a22*b33 *k1 - ---*a33*b33 *k1
                 8                 8
k19=----------------------------------------------
                            3
                         a22


k18=0


k17=0


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


k14=0


k13=0


      1            1
     ---*a22*k1 + ---*a33*k1
      2            2
k12=-------------------------
               a22


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:
 q1,k38,k16,k1,b33,n3,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   2
      + 4*a22 *b33*u3*v3 + 4*a22 *n1*u1 + 4*a22 *n3*u3 - 4*a22*a33 *u3

                                                                         2   2
      - 4*a22*a33*b33*u3*v3 - 4*a22*a33*n1*u1 - 4*a22*a33*n3*u3 - a22*b33 *v3

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

       *a22*(a22 - a33))

           7      4         7      2               7            2
INT=(32*a22 *k1*u1  + 64*a22 *k1*u1 *u2*u3 + 32*a22 *k1*u1*u2*v2

               7         2         7            2          7
      - 128*a22 *k1*u1*u3  + 16*a22 *k1*u1*v1*v2  - 128*a22 *k1*u1*v1

              7      2   2         7       3            7       3
      + 16*a22 *k1*u2 *u3  + 32*a22 *k16*u1 *v1 + 32*a22 *k16*u1 *v3

              7                         7          3         7       3
      + 32*a22 *k16*u1*u2*u3*v1 + 32*a22 *k16*u1*v2  + 32*a22 *k16*u2 *v1

              7          3         7          2            7
      + 32*a22 *k16*v1*v2  + 32*a22 *k38*u1*u2 *v2 - 64*a22 *k38*u1*u3*v1

              7                   7       3            7
      - 32*a22 *k38*u1*u3 + 32*a22 *k38*u2 *v3 + 32*a22 *k38*u2*u3*v2*v3

              7                      7       2         7      2
      + 64*a22 *k38*u2*v1*v2 - 32*a22 *k38*v1  + 32*a22 *q1*u1

              7                  6          4          6          2
      + 32*a22 *q1*u1*v1 - 64*a22 *a33*k1*u1  - 128*a22 *a33*k1*u1 *u2*u3

              6                2          6             2
      - 64*a22 *a33*k1*u1*u2*v2  + 256*a22 *a33*k1*u1*u3

              6                2          6
      - 16*a22 *a33*k1*u1*v1*v2  + 256*a22 *a33*k1*u1*v1

              6          2   2         6           3            6           3
      - 16*a22 *a33*k1*u2 *u3  - 64*a22 *a33*k16*u1 *v1 - 64*a22 *a33*k16*u1 *v3

              6                             6              3
      - 64*a22 *a33*k16*u1*u2*u3*v1 - 64*a22 *a33*k16*u1*v2

              6           3            6              3
      - 64*a22 *a33*k16*u2 *v1 - 64*a22 *a33*k16*v1*v2

              6              2             6
      - 64*a22 *a33*k38*u1*u2 *v2 + 128*a22 *a33*k38*u1*u3*v1

              6                       6           3
      + 64*a22 *a33*k38*u1*u3 - 64*a22 *a33*k38*u2 *v3

              6                              6
      - 64*a22 *a33*k38*u2*u3*v2*v3 - 128*a22 *a33*k38*u2*v1*v2

              6           2         6          2         6
      + 64*a22 *a33*k38*v1  - 64*a22 *a33*q1*u1  - 64*a22 *a33*q1*u1*v1

              6                          6
      + 64*a22 *b33*k1*u1*u3*v1 + 192*a22 *b33*k1*u1*u3

              6                            6
      + 16*a22 *b33*k1*u2*u3*v2*v3 - 64*a22 *b33*k1*u2*v1*v2

              6          4         6          2   2          6          2
      + 64*a22 *b33*k1*u3  + 64*a22 *b33*k1*v1 *v2  + 192*a22 *b33*k1*v1

              6                             6              2
      - 32*a22 *b33*k16*u1*u2*v1*v2 - 32*a22 *b33*k16*u1*v1 *v2

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

              6              2            6           3
      - 32*a22 *b33*k16*u2*v1 *v3 - 32*a22 *b33*k16*v1 *v3

              6                      6          2         6         3
      - 32*a22 *b33*q1*u1*u3 - 32*a22 *b33*q1*v1  + 48*a22 *k1*n1*u1

              6                        6                  6         2
      + 16*a22 *k1*n1*u1*u2*v1 + 16*a22 *k1*n1*u1 + 48*a22 *k1*n1*u2 *u3

              6                     6            2         6
      - 32*a22 *k1*n1*u2*v1 + 48*a22 *k1*n3*u1*u2  + 16*a22 *k1*n3*u1*u2*v3

              6         2            6          2            6          2
      + 48*a22 *k1*n3*u2 *v2 + 32*a22 *k16*n1*u1 *v2 + 32*a22 *k16*n1*u3 *v1

              6             2         6                         6          2
      + 32*a22 *k16*n1*u3*v1  + 32*a22 *k16*n3*u1*v2*v3 + 32*a22 *k16*n3*u3 *v3

              6                         6          3         6
      + 32*a22 *k16*n3*u3*v1*v3 - 32*a22 *k38*n1*u2  - 32*a22 *k38*n1*u2*u3*v1

              6          3         6          2            6
      - 32*a22 *k38*n1*u3  - 32*a22 *k38*n3*u2 *v1 - 32*a22 *k38*n3*u2*u3*v3

              6                      6                         6
      + 32*a22 *k38*n3*u2*u3 - 32*a22 *k38*n3*u3*v1*v2 - 32*a22 *k38*n3*u3*v1

              6                  6                  5    2      4
      + 32*a22 *n1*q1*u1 + 32*a22 *n3*q1*u3 + 32*a22 *a33 *k1*u1

              5    2      2               5    2            2
      + 64*a22 *a33 *k1*u1 *u2*u3 + 32*a22 *a33 *k1*u1*u2*v2

               5    2         2         5    2            2
      - 128*a22 *a33 *k1*u1*u3  - 16*a22 *a33 *k1*u1*v1*v2

               5    2                  5    2      2   2
      - 128*a22 *a33 *k1*u1*v1 - 16*a22 *a33 *k1*u2 *u3

              5    2       3            5    2       3
      + 32*a22 *a33 *k16*u1 *v1 + 32*a22 *a33 *k16*u1 *v3

              5    2                         5    2          3
      + 32*a22 *a33 *k16*u1*u2*u3*v1 + 32*a22 *a33 *k16*u1*v2

              5    2       3            5    2          3
      + 32*a22 *a33 *k16*u2 *v1 + 32*a22 *a33 *k16*v1*v2

              5    2          2            5    2
      + 32*a22 *a33 *k38*u1*u2 *v2 - 64*a22 *a33 *k38*u1*u3*v1

              5    2                   5    2       3
      - 32*a22 *a33 *k38*u1*u3 + 32*a22 *a33 *k38*u2 *v3

              5    2                         5    2
      + 32*a22 *a33 *k38*u2*u3*v2*v3 + 64*a22 *a33 *k38*u2*v1*v2

              5    2       2         5    2      2         5    2
      - 32*a22 *a33 *k38*v1  + 32*a22 *a33 *q1*u1  + 32*a22 *a33 *q1*u1*v1

               5                              5
      - 128*a22 *a33*b33*k1*u1*u3*v1 - 384*a22 *a33*b33*k1*u1*u3

              5                                 5
      - 16*a22 *a33*b33*k1*u2*u3*v2*v3 + 128*a22 *a33*b33*k1*u2*v1*v2

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

               5              2         5
      - 384*a22 *a33*b33*k1*v1  + 64*a22 *a33*b33*k16*u1*u2*v1*v2

              5                  2            5               2
      + 64*a22 *a33*b33*k16*u1*v1 *v2 + 64*a22 *a33*b33*k16*u2 *v1*v3

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

              5               3            5
      + 64*a22 *a33*b33*k16*v1 *v3 + 64*a22 *a33*b33*q1*u1*u3

              5              2         5             3         5
      + 64*a22 *a33*b33*q1*v1  - 64*a22 *a33*k1*n1*u1  - 32*a22 *a33*k1*n1*u1

              5             2            5
      - 64*a22 *a33*k1*n1*u2 *u3 + 64*a22 *a33*k1*n1*u2*v1

              5                2         5             2
      - 64*a22 *a33*k1*n3*u1*u2  - 64*a22 *a33*k1*n3*u2 *v2

              5              2            5              2
      - 32*a22 *a33*k16*n1*u1 *v2 - 32*a22 *a33*k16*n1*u3 *v1

              5                 2         5
      - 32*a22 *a33*k16*n1*u3*v1  - 32*a22 *a33*k16*n3*u1*v2*v3

              5              2            5
      - 32*a22 *a33*k16*n3*u3 *v3 - 32*a22 *a33*k16*n3*u3*v1*v3

              5              3         5
      + 32*a22 *a33*k38*n1*u2  + 32*a22 *a33*k38*n1*u2*u3*v1

              5              3         5              2
      + 32*a22 *a33*k38*n1*u3  + 32*a22 *a33*k38*n3*u2 *v1

              5                             5
      + 32*a22 *a33*k38*n3*u2*u3*v3 - 32*a22 *a33*k38*n3*u2*u3

              5                             5
      + 32*a22 *a33*k38*n3*u3*v1*v2 + 32*a22 *a33*k38*n3*u3*v1

              5                      5                      5    2      3
      - 32*a22 *a33*n1*q1*u1 - 32*a22 *a33*n3*q1*u3 + 12*a22 *b33 *k1*u1 *v3

              5    2                        5    2         2
      + 12*a22 *b33 *k1*u1*u2*u3*v1 - 56*a22 *b33 *k1*u1*u3 *v2

              5    2         3         5    2
      - 56*a22 *b33 *k1*u1*v1  + 28*a22 *b33 *k1*u1*v1*v3

              5    2         3         5    2                  5    2      3
      + 28*a22 *b33 *k1*u1*v2  - 64*a22 *b33 *k1*u1*v3 + 28*a22 *b33 *k1*u2 *v1

              5    2      3           5    2         2        5    2      3
      - 56*a22 *b33 *k1*u3 *v3 + 4*a22 *b33 *k1*u3*v2  - 4*a22 *b33 *k1*v1

              5    2                 5    2             2
      - 64*a22 *b33 *k1*v1*v3 + 8*a22 *b33 *k16*u1*v1*v3

             5    2             2        5    2          3
      + 8*a22 *b33 *k16*u2*v2*v3  + 8*a22 *b33 *k16*u2*v3

              5    2       2   2        5    2          3        5    2
      + 16*a22 *b33 *k16*v1 *v3  + 8*a22 *b33 *k16*v1*v3  + 8*a22 *b33 *q1*u1*v3

             5    2                  5
      + 8*a22 *b33 *q1*v1*v3 - 40*a22 *b33*k1*n1*u1*u3*v2

              5                         5             3
      - 16*a22 *b33*k1*n1*u1*v2 - 40*a22 *b33*k1*n1*u2

             5             2            5
      - 8*a22 *b33*k1*n1*u2 *v3 + 48*a22 *b33*k1*n1*u2*u3*v1

              5                         5
      + 32*a22 *b33*k1*n1*v1*v2 + 88*a22 *b33*k1*n1*v1

              5                            5
      - 64*a22 *b33*k1*n3*u1*u2*u3 + 32*a22 *b33*k1*n3*u2*u3*v3

              5                         5             2
      - 64*a22 *b33*k1*n3*u2*u3 - 64*a22 *b33*k1*n3*u3 *v2

              5                          5
      + 64*a22 *b33*k1*n3*u3*v1 + 128*a22 *b33*k1*n3*v3

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

              5              2            5
      - 16*a22 *b33*k16*n1*v2 *v3 - 32*a22 *b33*k16*n3*v1*v2*v3

              5                 2         5              3
      - 32*a22 *b33*k16*n3*v1*v3  - 32*a22 *b33*k16*n3*v3

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

              5              2            5                      5
      - 16*a22 *b33*k38*n1*u3 *v1 - 16*a22 *b33*n1*q1*v1 - 32*a22 *b33*n3*q1*v3

              5      2               5            2         5
      - 16*a22 *k1*n1 *u1*v1 + 32*a22 *k1*n1*n3*u2  + 16*a22 *k1*n1*n3*u2

              4    3            2         4    3      2   2
      + 16*a22 *a33 *k1*u1*v1*v2  + 16*a22 *a33 *k1*u2 *u3

              4    2                          4    2
      + 64*a22 *a33 *b33*k1*u1*u3*v1 + 192*a22 *a33 *b33*k1*u1*u3

              4    2                            4    2
      - 16*a22 *a33 *b33*k1*u2*u3*v2*v3 - 64*a22 *a33 *b33*k1*u2*v1*v2

              4    2          4         4    2          2   2
      + 64*a22 *a33 *b33*k1*u3  + 64*a22 *a33 *b33*k1*v1 *v2

               4    2          2         4    2
      + 192*a22 *a33 *b33*k1*v1  - 32*a22 *a33 *b33*k16*u1*u2*v1*v2

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

              4    2                             4    2              2
      - 32*a22 *a33 *b33*k16*u2*u3*v1*v2 - 32*a22 *a33 *b33*k16*u2*v1 *v3

              4    2           3            4    2
      - 32*a22 *a33 *b33*k16*v1 *v3 - 32*a22 *a33 *b33*q1*u1*u3

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

              4    2                        4    2
      - 16*a22 *a33 *k1*n1*u1*u2*v1 + 16*a22 *a33 *k1*n1*u1

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

              4    2            2         4    2
      + 16*a22 *a33 *k1*n3*u1*u2  - 16*a22 *a33 *k1*n3*u1*u2*v3

              4    2         2            4        2      3
      + 16*a22 *a33 *k1*n3*u2 *v2 - 28*a22 *a33*b33 *k1*u1 *v3

              4        2                         4        2         2
      - 28*a22 *a33*b33 *k1*u1*u2*u3*v1 + 120*a22 *a33*b33 *k1*u1*u3 *v2

               4        2         3         4        2
      + 120*a22 *a33*b33 *k1*u1*v1  - 60*a22 *a33*b33 *k1*u1*v1*v3

              4        2         3          4        2
      - 60*a22 *a33*b33 *k1*u1*v2  + 128*a22 *a33*b33 *k1*u1*v3

              4        2      3             4        2      3
      - 60*a22 *a33*b33 *k1*u2 *v1 + 120*a22 *a33*b33 *k1*u3 *v3

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

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

              4        2             2         4        2          3
      - 16*a22 *a33*b33 *k16*u2*v2*v3  - 16*a22 *a33*b33 *k16*u2*v3

              4        2       2   2         4        2          3
      - 32*a22 *a33*b33 *k16*v1 *v3  - 16*a22 *a33*b33 *k16*v1*v3

              4        2                  4        2
      - 16*a22 *a33*b33 *q1*u1*v3 - 16*a22 *a33*b33 *q1*v1*v3

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

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

              4                          4
      - 80*a22 *a33*b33*k1*n1*v1 + 64*a22 *a33*b33*k1*n3*u1*u2*u3

              4                                4
      - 64*a22 *a33*b33*k1*n3*u2*u3*v3 + 64*a22 *a33*b33*k1*n3*u2*u3

              4                 2            4
      + 64*a22 *a33*b33*k1*n3*u3 *v2 - 64*a22 *a33*b33*k1*n3*u3*v1

               4                          4
      - 128*a22 *a33*b33*k1*n3*v3 + 16*a22 *a33*b33*k16*n1*u2*v2*v3

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

              4                                 4                     2
      + 32*a22 *a33*b33*k16*n3*v1*v2*v3 + 32*a22 *a33*b33*k16*n3*v1*v3

              4                  3         4
      + 32*a22 *a33*b33*k16*n3*v3  + 16*a22 *a33*b33*k38*n1*u2*v1*v3

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

              4                          4
      + 16*a22 *a33*b33*n1*q1*v1 + 32*a22 *a33*b33*n3*q1*v3

              4          2               4                2
      + 16*a22 *a33*k1*n1 *u1*v1 - 32*a22 *a33*k1*n1*n3*u2

              4                         4    3
      - 16*a22 *a33*k1*n1*n3*u2 - 12*a22 *b33 *k1*u1*u2*v1*v2

              4    3            2         4    3         2
      + 16*a22 *b33 *k1*u1*u3*v3  - 12*a22 *b33 *k1*u1*v1 *v2

              4    3      2               4    3
      - 12*a22 *b33 *k1*u2 *v1*v3 - 12*a22 *b33 *k1*u2*u3*v1*v2

              4    3         2            4    3      3
      - 12*a22 *b33 *k1*u2*v1 *v3 - 12*a22 *b33 *k1*v1 *v3

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

              4    2                        4    2         2
      - 40*a22 *b33 *k1*n1*u2*v1*v3 + 28*a22 *b33 *k1*n1*u3 *v1

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

              4    2                        4    2         2
      - 48*a22 *b33 *k1*n3*u1*v2*v3 + 28*a22 *b33 *k1*n3*u3 *v3

              4    2                        4          2
      + 28*a22 *b33 *k1*n3*u3*v1*v3 - 16*a22 *b33*k1*n1 *u1*u3

              4          2   2         4                2
      + 16*a22 *b33*k1*n1 *v1  - 48*a22 *b33*k1*n1*n3*u3

              4                           4
      - 32*a22 *b33*k1*n1*n3*u3*v2 - 8*a22 *b33*k1*n1*n3*v2

              4          2               4          2   2         4      3
      + 48*a22 *b33*k1*n3 *u1*u3 + 48*a22 *b33*k1*n3 *v1  - 16*a22 *k1*n1 *u1

              4      2               4         2            4      3
      - 16*a22 *k1*n1 *n3*u3 + 16*a22 *k1*n1*n3 *u1 + 16*a22 *k1*n3 *u3

              3    3                            3    2    2      3
      + 16*a22 *a33 *b33*k1*u2*u3*v2*v3 + 20*a22 *a33 *b33 *k1*u1 *v3

              3    2    2                        3    2    2         2
      + 20*a22 *a33 *b33 *k1*u1*u2*u3*v1 - 72*a22 *a33 *b33 *k1*u1*u3 *v2

              3    2    2         3         3    2    2
      - 72*a22 *a33 *b33 *k1*u1*v1  + 36*a22 *a33 *b33 *k1*u1*v1*v3

              3    2    2         3         3    2    2
      + 36*a22 *a33 *b33 *k1*u1*v2  - 64*a22 *a33 *b33 *k1*u1*v3

              3    2    2      3            3    2    2      3
      + 36*a22 *a33 *b33 *k1*u2 *v1 - 72*a22 *a33 *b33 *k1*u3 *v3

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

              3    2    2                 3    2    2             2
      - 64*a22 *a33 *b33 *k1*v1*v3 + 8*a22 *a33 *b33 *k16*u1*v1*v3

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

              3    2    2       2   2        3    2    2          3
      + 16*a22 *a33 *b33 *k16*v1 *v3  + 8*a22 *a33 *b33 *k16*v1*v3

             3    2    2                 3    2    2
      + 8*a22 *a33 *b33 *q1*u1*v3 + 8*a22 *a33 *b33 *q1*v1*v3

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

             3    2             3        3    2             2
      + 8*a22 *a33 *b33*k1*n1*u2  + 8*a22 *a33 *b33*k1*n1*u2 *v3

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

              3    2                            3        3
      + 32*a22 *a33 *b33*k1*n3*u2*u3*v3 + 28*a22 *a33*b33 *k1*u1*u2*v1*v2

              3        3            2         3        3         2
      - 32*a22 *a33*b33 *k1*u1*u3*v3  + 28*a22 *a33*b33 *k1*u1*v1 *v2

              3        3      2               3        3
      + 28*a22 *a33*b33 *k1*u2 *v1*v3 + 28*a22 *a33*b33 *k1*u2*u3*v1*v2

              3        3         2            3        3      3
      + 28*a22 *a33*b33 *k1*u2*v1 *v3 + 28*a22 *a33*b33 *k1*v1 *v3

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

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

              3        2            2         3        2
      - 32*a22 *a33*b33 *k1*n1*u3*v1  + 48*a22 *a33*b33 *k1*n3*u1*v2*v3

              3        2         2            3        2
      - 32*a22 *a33*b33 *k1*n3*u3 *v3 - 32*a22 *a33*b33 *k1*n3*u3*v1*v3

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

              3                    2        3
      + 16*a22 *a33*b33*k1*n1*n3*u3  + 8*a22 *a33*b33*k1*n1*n3*v2

              3              2               3              2   2
      - 16*a22 *a33*b33*k1*n3 *u1*u3 - 16*a22 *a33*b33*k1*n3 *v1

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

             3    4         3         3    4      2   2        3    4         3
      + 5*a22 *b33 *k1*u2*v3  + 10*a22 *b33 *k1*v1 *v3  + 3*a22 *b33 *k1*v1*v3

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

              3    3         2            3    3
      - 14*a22 *b33 *k1*n1*v2 *v3 - 28*a22 *b33 *k1*n3*v1*v2*v3

              3    3            2         3    3         3
      - 28*a22 *b33 *k1*n3*v1*v3  - 12*a22 *b33 *k1*n3*v3

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

              3    2                        3    2      2
      + 40*a22 *b33 *k1*n1*n3*u3*v3 - 32*a22 *b33 *k1*n3 *u1*v3

              3    2      2              3          3
      - 32*a22 *b33 *k1*n3 *v1*v3 + 8*a22 *b33*k1*n1 *v1

              3          2              3             2
      + 16*a22 *b33*k1*n1 *n3*v3 + 8*a22 *b33*k1*n1*n3 *v1

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

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

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

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

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

              2    2    3                        2    2    3            2
      - 20*a22 *a33 *b33 *k1*u1*u2*v1*v2 + 16*a22 *a33 *b33 *k1*u1*u3*v3

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

              2    2    3                        2    2    3         2
      - 20*a22 *a33 *b33 *k1*u2*u3*v1*v2 - 20*a22 *a33 *b33 *k1*u2*v1 *v3

              2    2    3      3            2    2    3         2
      - 20*a22 *a33 *b33 *k1*v1 *v3 + 16*a22 *a33 *b33 *k1*v1*v2 *v3

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

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

             2    2    2         2           2    2    2
      + 4*a22 *a33 *b33 *k1*n3*u3 *v3 + 4*a22 *a33 *b33 *k1*n3*u3*v1*v3

              2        4            2        2        4            2
      - 11*a22 *a33*b33 *k1*u1*v1*v3  - 7*a22 *a33*b33 *k1*u2*v2*v3

              2        4         3         2        4      2   2
      - 11*a22 *a33*b33 *k1*u2*v3  - 22*a22 *a33*b33 *k1*v1 *v3

             2        4         3         2        3
      - 7*a22 *a33*b33 *k1*v1*v3  + 16*a22 *a33*b33 *k1*n1*u2*v2*v3

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

              2        3                        2        3            2
      + 32*a22 *a33*b33 *k1*n3*v1*v2*v3 + 32*a22 *a33*b33 *k1*n3*v1*v3

              2        3         3        2        2      2
      + 16*a22 *a33*b33 *k1*n3*v3  + 8*a22 *a33*b33 *k1*n1 *u1*v3

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

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

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

                 3    3         2               3    3      3
      + 4*a22*a33 *b33 *k1*u2*v1 *v3 + 4*a22*a33 *b33 *k1*v1 *v3

                 2    4            2            2    4            2
      + 7*a22*a33 *b33 *k1*u1*v1*v3  + 5*a22*a33 *b33 *k1*u2*v2*v3

                 2    4         3             2    4      2   2
      + 7*a22*a33 *b33 *k1*u2*v3  + 14*a22*a33 *b33 *k1*v1 *v3

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

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

                 2    3                           2    3            2
      - 4*a22*a33 *b33 *k1*n3*v1*v2*v3 - 4*a22*a33 *b33 *k1*n3*v1*v3

                 2    3         3      3    4            2
      - 4*a22*a33 *b33 *k1*n3*v3  - a33 *b33 *k1*u1*v1*v3

           3    4            2      3    4         3        3    4      2   2
      - a33 *b33 *k1*u2*v2*v3  - a33 *b33 *k1*u2*v3  - 2*a33 *b33 *k1*v1 *v3

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

INT=(32*a22**7*k1*u1**4 + 64*a22**7*k1*u1**2*u2*u3 + 32*a22**7*k1*u1*u2*v2**2 - 
128*a22**7*k1*u1*u3**2 + 16*a22**7*k1*u1*v1*v2**2 - 128*a22**7*k1*u1*v1 + 16*a22
**7*k1*u2**2*u3**2 + 32*a22**7*k16*u1**3*v1 + 32*a22**7*k16*u1**3*v3 + 32*a22**7
*k16*u1*u2*u3*v1 + 32*a22**7*k16*u1*v2**3 + 32*a22**7*k16*u2**3*v1 + 32*a22**7*
k16*v1*v2**3 + 32*a22**7*k38*u1*u2**2*v2 - 64*a22**7*k38*u1*u3*v1 - 32*a22**7*
k38*u1*u3 + 32*a22**7*k38*u2**3*v3 + 32*a22**7*k38*u2*u3*v2*v3 + 64*a22**7*k38*
u2*v1*v2 - 32*a22**7*k38*v1**2 + 32*a22**7*q1*u1**2 + 32*a22**7*q1*u1*v1 - 64*
a22**6*a33*k1*u1**4 - 128*a22**6*a33*k1*u1**2*u2*u3 - 64*a22**6*a33*k1*u1*u2*v2
**2 + 256*a22**6*a33*k1*u1*u3**2 - 16*a22**6*a33*k1*u1*v1*v2**2 + 256*a22**6*a33
*k1*u1*v1 - 16*a22**6*a33*k1*u2**2*u3**2 - 64*a22**6*a33*k16*u1**3*v1 - 64*a22**
6*a33*k16*u1**3*v3 - 64*a22**6*a33*k16*u1*u2*u3*v1 - 64*a22**6*a33*k16*u1*v2**3 
- 64*a22**6*a33*k16*u2**3*v1 - 64*a22**6*a33*k16*v1*v2**3 - 64*a22**6*a33*k38*u1
*u2**2*v2 + 128*a22**6*a33*k38*u1*u3*v1 + 64*a22**6*a33*k38*u1*u3 - 64*a22**6*
a33*k38*u2**3*v3 - 64*a22**6*a33*k38*u2*u3*v2*v3 - 128*a22**6*a33*k38*u2*v1*v2 +
 64*a22**6*a33*k38*v1**2 - 64*a22**6*a33*q1*u1**2 - 64*a22**6*a33*q1*u1*v1 + 64*
a22**6*b33*k1*u1*u3*v1 + 192*a22**6*b33*k1*u1*u3 + 16*a22**6*b33*k1*u2*u3*v2*v3 
- 64*a22**6*b33*k1*u2*v1*v2 + 64*a22**6*b33*k1*u3**4 + 64*a22**6*b33*k1*v1**2*v2
**2 + 192*a22**6*b33*k1*v1**2 - 32*a22**6*b33*k16*u1*u2*v1*v2 - 32*a22**6*b33*
k16*u1*v1**2*v2 - 32*a22**6*b33*k16*u2**2*v1*v3 - 32*a22**6*b33*k16*u2*u3*v1*v2 
- 32*a22**6*b33*k16*u2*v1**2*v3 - 32*a22**6*b33*k16*v1**3*v3 - 32*a22**6*b33*q1*
u1*u3 - 32*a22**6*b33*q1*v1**2 + 48*a22**6*k1*n1*u1**3 + 16*a22**6*k1*n1*u1*u2*
v1 + 16*a22**6*k1*n1*u1 + 48*a22**6*k1*n1*u2**2*u3 - 32*a22**6*k1*n1*u2*v1 + 48*
a22**6*k1*n3*u1*u2**2 + 16*a22**6*k1*n3*u1*u2*v3 + 48*a22**6*k1*n3*u2**2*v2 + 32
*a22**6*k16*n1*u1**2*v2 + 32*a22**6*k16*n1*u3**2*v1 + 32*a22**6*k16*n1*u3*v1**2 
+ 32*a22**6*k16*n3*u1*v2*v3 + 32*a22**6*k16*n3*u3**2*v3 + 32*a22**6*k16*n3*u3*v1
*v3 - 32*a22**6*k38*n1*u2**3 - 32*a22**6*k38*n1*u2*u3*v1 - 32*a22**6*k38*n1*u3**
3 - 32*a22**6*k38*n3*u2**2*v1 - 32*a22**6*k38*n3*u2*u3*v3 + 32*a22**6*k38*n3*u2*
u3 - 32*a22**6*k38*n3*u3*v1*v2 - 32*a22**6*k38*n3*u3*v1 + 32*a22**6*n1*q1*u1 + 
32*a22**6*n3*q1*u3 + 32*a22**5*a33**2*k1*u1**4 + 64*a22**5*a33**2*k1*u1**2*u2*u3
 + 32*a22**5*a33**2*k1*u1*u2*v2**2 - 128*a22**5*a33**2*k1*u1*u3**2 - 16*a22**5*
a33**2*k1*u1*v1*v2**2 - 128*a22**5*a33**2*k1*u1*v1 - 16*a22**5*a33**2*k1*u2**2*
u3**2 + 32*a22**5*a33**2*k16*u1**3*v1 + 32*a22**5*a33**2*k16*u1**3*v3 + 32*a22**
5*a33**2*k16*u1*u2*u3*v1 + 32*a22**5*a33**2*k16*u1*v2**3 + 32*a22**5*a33**2*k16*
u2**3*v1 + 32*a22**5*a33**2*k16*v1*v2**3 + 32*a22**5*a33**2*k38*u1*u2**2*v2 - 64
*a22**5*a33**2*k38*u1*u3*v1 - 32*a22**5*a33**2*k38*u1*u3 + 32*a22**5*a33**2*k38*
u2**3*v3 + 32*a22**5*a33**2*k38*u2*u3*v2*v3 + 64*a22**5*a33**2*k38*u2*v1*v2 - 32
*a22**5*a33**2*k38*v1**2 + 32*a22**5*a33**2*q1*u1**2 + 32*a22**5*a33**2*q1*u1*v1
 - 128*a22**5*a33*b33*k1*u1*u3*v1 - 384*a22**5*a33*b33*k1*u1*u3 - 16*a22**5*a33*
b33*k1*u2*u3*v2*v3 + 128*a22**5*a33*b33*k1*u2*v1*v2 - 128*a22**5*a33*b33*k1*u3**
4 - 128*a22**5*a33*b33*k1*v1**2*v2**2 - 384*a22**5*a33*b33*k1*v1**2 + 64*a22**5*
a33*b33*k16*u1*u2*v1*v2 + 64*a22**5*a33*b33*k16*u1*v1**2*v2 + 64*a22**5*a33*b33*
k16*u2**2*v1*v3 + 64*a22**5*a33*b33*k16*u2*u3*v1*v2 + 64*a22**5*a33*b33*k16*u2*
v1**2*v3 + 64*a22**5*a33*b33*k16*v1**3*v3 + 64*a22**5*a33*b33*q1*u1*u3 + 64*a22
**5*a33*b33*q1*v1**2 - 64*a22**5*a33*k1*n1*u1**3 - 32*a22**5*a33*k1*n1*u1 - 64*
a22**5*a33*k1*n1*u2**2*u3 + 64*a22**5*a33*k1*n1*u2*v1 - 64*a22**5*a33*k1*n3*u1*
u2**2 - 64*a22**5*a33*k1*n3*u2**2*v2 - 32*a22**5*a33*k16*n1*u1**2*v2 - 32*a22**5
*a33*k16*n1*u3**2*v1 - 32*a22**5*a33*k16*n1*u3*v1**2 - 32*a22**5*a33*k16*n3*u1*
v2*v3 - 32*a22**5*a33*k16*n3*u3**2*v3 - 32*a22**5*a33*k16*n3*u3*v1*v3 + 32*a22**
5*a33*k38*n1*u2**3 + 32*a22**5*a33*k38*n1*u2*u3*v1 + 32*a22**5*a33*k38*n1*u3**3 
+ 32*a22**5*a33*k38*n3*u2**2*v1 + 32*a22**5*a33*k38*n3*u2*u3*v3 - 32*a22**5*a33*
k38*n3*u2*u3 + 32*a22**5*a33*k38*n3*u3*v1*v2 + 32*a22**5*a33*k38*n3*u3*v1 - 32*
a22**5*a33*n1*q1*u1 - 32*a22**5*a33*n3*q1*u3 + 12*a22**5*b33**2*k1*u1**3*v3 + 12
*a22**5*b33**2*k1*u1*u2*u3*v1 - 56*a22**5*b33**2*k1*u1*u3**2*v2 - 56*a22**5*b33
**2*k1*u1*v1**3 + 28*a22**5*b33**2*k1*u1*v1*v3 + 28*a22**5*b33**2*k1*u1*v2**3 - 
64*a22**5*b33**2*k1*u1*v3 + 28*a22**5*b33**2*k1*u2**3*v1 - 56*a22**5*b33**2*k1*
u3**3*v3 + 4*a22**5*b33**2*k1*u3*v2**2 - 4*a22**5*b33**2*k1*v1**3 - 64*a22**5*
b33**2*k1*v1*v3 + 8*a22**5*b33**2*k16*u1*v1*v3**2 + 8*a22**5*b33**2*k16*u2*v2*v3
**2 + 8*a22**5*b33**2*k16*u2*v3**3 + 16*a22**5*b33**2*k16*v1**2*v3**2 + 8*a22**5
*b33**2*k16*v1*v3**3 + 8*a22**5*b33**2*q1*u1*v3 + 8*a22**5*b33**2*q1*v1*v3 - 40*
a22**5*b33*k1*n1*u1*u3*v2 - 16*a22**5*b33*k1*n1*u1*v2 - 40*a22**5*b33*k1*n1*u2**
3 - 8*a22**5*b33*k1*n1*u2**2*v3 + 48*a22**5*b33*k1*n1*u2*u3*v1 + 32*a22**5*b33*
k1*n1*v1*v2 + 88*a22**5*b33*k1*n1*v1 - 64*a22**5*b33*k1*n3*u1*u2*u3 + 32*a22**5*
b33*k1*n3*u2*u3*v3 - 64*a22**5*b33*k1*n3*u2*u3 - 64*a22**5*b33*k1*n3*u3**2*v2 + 
64*a22**5*b33*k1*n3*u3*v1 + 128*a22**5*b33*k1*n3*v3 - 16*a22**5*b33*k16*n1*u2*v2
*v3 - 16*a22**5*b33*k16*n1*u2*v3**2 - 16*a22**5*b33*k16*n1*v2**2*v3 - 32*a22**5*
b33*k16*n3*v1*v2*v3 - 32*a22**5*b33*k16*n3*v1*v3**2 - 32*a22**5*b33*k16*n3*v3**3
 - 16*a22**5*b33*k38*n1*u2*v1*v3 - 16*a22**5*b33*k38*n1*u2*v2**2 - 16*a22**5*b33
*k38*n1*u3**2*v1 - 16*a22**5*b33*n1*q1*v1 - 32*a22**5*b33*n3*q1*v3 - 16*a22**5*
k1*n1**2*u1*v1 + 32*a22**5*k1*n1*n3*u2**2 + 16*a22**5*k1*n1*n3*u2 + 16*a22**4*
a33**3*k1*u1*v1*v2**2 + 16*a22**4*a33**3*k1*u2**2*u3**2 + 64*a22**4*a33**2*b33*
k1*u1*u3*v1 + 192*a22**4*a33**2*b33*k1*u1*u3 - 16*a22**4*a33**2*b33*k1*u2*u3*v2*
v3 - 64*a22**4*a33**2*b33*k1*u2*v1*v2 + 64*a22**4*a33**2*b33*k1*u3**4 + 64*a22**
4*a33**2*b33*k1*v1**2*v2**2 + 192*a22**4*a33**2*b33*k1*v1**2 - 32*a22**4*a33**2*
b33*k16*u1*u2*v1*v2 - 32*a22**4*a33**2*b33*k16*u1*v1**2*v2 - 32*a22**4*a33**2*
b33*k16*u2**2*v1*v3 - 32*a22**4*a33**2*b33*k16*u2*u3*v1*v2 - 32*a22**4*a33**2*
b33*k16*u2*v1**2*v3 - 32*a22**4*a33**2*b33*k16*v1**3*v3 - 32*a22**4*a33**2*b33*
q1*u1*u3 - 32*a22**4*a33**2*b33*q1*v1**2 + 16*a22**4*a33**2*k1*n1*u1**3 - 16*a22
**4*a33**2*k1*n1*u1*u2*v1 + 16*a22**4*a33**2*k1*n1*u1 + 16*a22**4*a33**2*k1*n1*
u2**2*u3 - 32*a22**4*a33**2*k1*n1*u2*v1 + 16*a22**4*a33**2*k1*n3*u1*u2**2 - 16*
a22**4*a33**2*k1*n3*u1*u2*v3 + 16*a22**4*a33**2*k1*n3*u2**2*v2 - 28*a22**4*a33*
b33**2*k1*u1**3*v3 - 28*a22**4*a33*b33**2*k1*u1*u2*u3*v1 + 120*a22**4*a33*b33**2
*k1*u1*u3**2*v2 + 120*a22**4*a33*b33**2*k1*u1*v1**3 - 60*a22**4*a33*b33**2*k1*u1
*v1*v3 - 60*a22**4*a33*b33**2*k1*u1*v2**3 + 128*a22**4*a33*b33**2*k1*u1*v3 - 60*
a22**4*a33*b33**2*k1*u2**3*v1 + 120*a22**4*a33*b33**2*k1*u3**3*v3 - 4*a22**4*a33
*b33**2*k1*u3*v2**2 + 4*a22**4*a33*b33**2*k1*v1**3 + 128*a22**4*a33*b33**2*k1*v1
*v3 - 16*a22**4*a33*b33**2*k16*u1*v1*v3**2 - 16*a22**4*a33*b33**2*k16*u2*v2*v3**
2 - 16*a22**4*a33*b33**2*k16*u2*v3**3 - 32*a22**4*a33*b33**2*k16*v1**2*v3**2 - 
16*a22**4*a33*b33**2*k16*v1*v3**3 - 16*a22**4*a33*b33**2*q1*u1*v3 - 16*a22**4*
a33*b33**2*q1*v1*v3 + 32*a22**4*a33*b33*k1*n1*u1*u3*v2 + 32*a22**4*a33*b33*k1*n1
*u2**3 - 64*a22**4*a33*b33*k1*n1*u2*u3*v1 - 32*a22**4*a33*b33*k1*n1*v1*v2 - 80*
a22**4*a33*b33*k1*n1*v1 + 64*a22**4*a33*b33*k1*n3*u1*u2*u3 - 64*a22**4*a33*b33*
k1*n3*u2*u3*v3 + 64*a22**4*a33*b33*k1*n3*u2*u3 + 64*a22**4*a33*b33*k1*n3*u3**2*
v2 - 64*a22**4*a33*b33*k1*n3*u3*v1 - 128*a22**4*a33*b33*k1*n3*v3 + 16*a22**4*a33
*b33*k16*n1*u2*v2*v3 + 16*a22**4*a33*b33*k16*n1*u2*v3**2 + 16*a22**4*a33*b33*k16
*n1*v2**2*v3 + 32*a22**4*a33*b33*k16*n3*v1*v2*v3 + 32*a22**4*a33*b33*k16*n3*v1*
v3**2 + 32*a22**4*a33*b33*k16*n3*v3**3 + 16*a22**4*a33*b33*k38*n1*u2*v1*v3 + 16*
a22**4*a33*b33*k38*n1*u2*v2**2 + 16*a22**4*a33*b33*k38*n1*u3**2*v1 + 16*a22**4*
a33*b33*n1*q1*v1 + 32*a22**4*a33*b33*n3*q1*v3 + 16*a22**4*a33*k1*n1**2*u1*v1 - 
32*a22**4*a33*k1*n1*n3*u2**2 - 16*a22**4*a33*k1*n1*n3*u2 - 12*a22**4*b33**3*k1*
u1*u2*v1*v2 + 16*a22**4*b33**3*k1*u1*u3*v3**2 - 12*a22**4*b33**3*k1*u1*v1**2*v2 
- 12*a22**4*b33**3*k1*u2**2*v1*v3 - 12*a22**4*b33**3*k1*u2*u3*v1*v2 - 12*a22**4*
b33**3*k1*u2*v1**2*v3 - 12*a22**4*b33**3*k1*v1**3*v3 + 16*a22**4*b33**3*k1*v1*v2
**2*v3 + 16*a22**4*b33**2*k1*n1*u1**2*v2 - 40*a22**4*b33**2*k1*n1*u2*v1*v3 + 28*
a22**4*b33**2*k1*n1*u3**2*v1 + 28*a22**4*b33**2*k1*n1*u3*v1**2 + 8*a22**4*b33**2
*k1*n1*v2*v3 - 48*a22**4*b33**2*k1*n3*u1*v2*v3 + 28*a22**4*b33**2*k1*n3*u3**2*v3
 + 28*a22**4*b33**2*k1*n3*u3*v1*v3 - 16*a22**4*b33*k1*n1**2*u1*u3 + 16*a22**4*
b33*k1*n1**2*v1**2 - 48*a22**4*b33*k1*n1*n3*u3**2 - 32*a22**4*b33*k1*n1*n3*u3*v2
 - 8*a22**4*b33*k1*n1*n3*v2 + 48*a22**4*b33*k1*n3**2*u1*u3 + 48*a22**4*b33*k1*n3
**2*v1**2 - 16*a22**4*k1*n1**3*u1 - 16*a22**4*k1*n1**2*n3*u3 + 16*a22**4*k1*n1*
n3**2*u1 + 16*a22**4*k1*n3**3*u3 + 16*a22**3*a33**3*b33*k1*u2*u3*v2*v3 + 20*a22
**3*a33**2*b33**2*k1*u1**3*v3 + 20*a22**3*a33**2*b33**2*k1*u1*u2*u3*v1 - 72*a22
**3*a33**2*b33**2*k1*u1*u3**2*v2 - 72*a22**3*a33**2*b33**2*k1*u1*v1**3 + 36*a22
**3*a33**2*b33**2*k1*u1*v1*v3 + 36*a22**3*a33**2*b33**2*k1*u1*v2**3 - 64*a22**3*
a33**2*b33**2*k1*u1*v3 + 36*a22**3*a33**2*b33**2*k1*u2**3*v1 - 72*a22**3*a33**2*
b33**2*k1*u3**3*v3 - 4*a22**3*a33**2*b33**2*k1*u3*v2**2 + 4*a22**3*a33**2*b33**2
*k1*v1**3 - 64*a22**3*a33**2*b33**2*k1*v1*v3 + 8*a22**3*a33**2*b33**2*k16*u1*v1*
v3**2 + 8*a22**3*a33**2*b33**2*k16*u2*v2*v3**2 + 8*a22**3*a33**2*b33**2*k16*u2*
v3**3 + 16*a22**3*a33**2*b33**2*k16*v1**2*v3**2 + 8*a22**3*a33**2*b33**2*k16*v1*
v3**3 + 8*a22**3*a33**2*b33**2*q1*u1*v3 + 8*a22**3*a33**2*b33**2*q1*v1*v3 + 8*
a22**3*a33**2*b33*k1*n1*u1*u3*v2 + 16*a22**3*a33**2*b33*k1*n1*u1*v2 + 8*a22**3*
a33**2*b33*k1*n1*u2**3 + 8*a22**3*a33**2*b33*k1*n1*u2**2*v3 + 16*a22**3*a33**2*
b33*k1*n1*u2*u3*v1 - 8*a22**3*a33**2*b33*k1*n1*v1 + 32*a22**3*a33**2*b33*k1*n3*
u2*u3*v3 + 28*a22**3*a33*b33**3*k1*u1*u2*v1*v2 - 32*a22**3*a33*b33**3*k1*u1*u3*
v3**2 + 28*a22**3*a33*b33**3*k1*u1*v1**2*v2 + 28*a22**3*a33*b33**3*k1*u2**2*v1*
v3 + 28*a22**3*a33*b33**3*k1*u2*u3*v1*v2 + 28*a22**3*a33*b33**3*k1*u2*v1**2*v3 +
 28*a22**3*a33*b33**3*k1*v1**3*v3 - 32*a22**3*a33*b33**3*k1*v1*v2**2*v3 - 16*a22
**3*a33*b33**2*k1*n1*u1**2*v2 + 32*a22**3*a33*b33**2*k1*n1*u2*v1*v3 - 32*a22**3*
a33*b33**2*k1*n1*u3**2*v1 - 32*a22**3*a33*b33**2*k1*n1*u3*v1**2 + 48*a22**3*a33*
b33**2*k1*n3*u1*v2*v3 - 32*a22**3*a33*b33**2*k1*n3*u3**2*v3 - 32*a22**3*a33*b33
**2*k1*n3*u3*v1*v3 - 16*a22**3*a33*b33*k1*n1**2*u1*u3 - 16*a22**3*a33*b33*k1*n1
**2*v1**2 + 16*a22**3*a33*b33*k1*n1*n3*u3**2 + 8*a22**3*a33*b33*k1*n1*n3*v2 - 16
*a22**3*a33*b33*k1*n3**2*u1*u3 - 16*a22**3*a33*b33*k1*n3**2*v1**2 + 5*a22**3*b33
**4*k1*u1*v1*v3**2 + 3*a22**3*b33**4*k1*u2*v2*v3**2 + 5*a22**3*b33**4*k1*u2*v3**
3 + 10*a22**3*b33**4*k1*v1**2*v3**2 + 3*a22**3*b33**4*k1*v1*v3**3 - 14*a22**3*
b33**3*k1*n1*u2*v2*v3 - 4*a22**3*b33**3*k1*n1*u2*v3**2 - 14*a22**3*b33**3*k1*n1*
v2**2*v3 - 28*a22**3*b33**3*k1*n3*v1*v2*v3 - 28*a22**3*b33**3*k1*n3*v1*v3**2 - 
12*a22**3*b33**3*k1*n3*v3**3 + 8*a22**3*b33**2*k1*n1**2*u1*v3 - 4*a22**3*b33**2*
k1*n1**2*v1*v3 + 40*a22**3*b33**2*k1*n1*n3*u3*v3 - 32*a22**3*b33**2*k1*n3**2*u1*
v3 - 32*a22**3*b33**2*k1*n3**2*v1*v3 + 8*a22**3*b33*k1*n1**3*v1 + 16*a22**3*b33*
k1*n1**2*n3*v3 + 8*a22**3*b33*k1*n1*n3**2*v1 - 4*a22**2*a33**3*b33**2*k1*u1**3*
v3 - 4*a22**2*a33**3*b33**2*k1*u1*u2*u3*v1 + 8*a22**2*a33**3*b33**2*k1*u1*u3**2*
v2 + 8*a22**2*a33**3*b33**2*k1*u1*v1**3 - 4*a22**2*a33**3*b33**2*k1*u1*v1*v3 - 4
*a22**2*a33**3*b33**2*k1*u1*v2**3 - 4*a22**2*a33**3*b33**2*k1*u2**3*v1 + 8*a22**
2*a33**3*b33**2*k1*u3**3*v3 + 4*a22**2*a33**3*b33**2*k1*u3*v2**2 - 4*a22**2*a33
**3*b33**2*k1*v1**3 - 20*a22**2*a33**2*b33**3*k1*u1*u2*v1*v2 + 16*a22**2*a33**2*
b33**3*k1*u1*u3*v3**2 - 20*a22**2*a33**2*b33**3*k1*u1*v1**2*v2 - 20*a22**2*a33**
2*b33**3*k1*u2**2*v1*v3 - 20*a22**2*a33**2*b33**3*k1*u2*u3*v1*v2 - 20*a22**2*a33
**2*b33**3*k1*u2*v1**2*v3 - 20*a22**2*a33**2*b33**3*k1*v1**3*v3 + 16*a22**2*a33
**2*b33**3*k1*v1*v2**2*v3 + 8*a22**2*a33**2*b33**2*k1*n1*u2*v1*v3 + 4*a22**2*a33
**2*b33**2*k1*n1*u3**2*v1 + 4*a22**2*a33**2*b33**2*k1*n1*u3*v1**2 - 8*a22**2*a33
**2*b33**2*k1*n1*v2*v3 + 4*a22**2*a33**2*b33**2*k1*n3*u3**2*v3 + 4*a22**2*a33**2
*b33**2*k1*n3*u3*v1*v3 - 11*a22**2*a33*b33**4*k1*u1*v1*v3**2 - 7*a22**2*a33*b33
**4*k1*u2*v2*v3**2 - 11*a22**2*a33*b33**4*k1*u2*v3**3 - 22*a22**2*a33*b33**4*k1*
v1**2*v3**2 - 7*a22**2*a33*b33**4*k1*v1*v3**3 + 16*a22**2*a33*b33**3*k1*n1*u2*v2
*v3 + 8*a22**2*a33*b33**3*k1*n1*u2*v3**2 + 16*a22**2*a33*b33**3*k1*n1*v2**2*v3 +
 32*a22**2*a33*b33**3*k1*n3*v1*v2*v3 + 32*a22**2*a33*b33**3*k1*n3*v1*v3**2 + 16*
a22**2*a33*b33**3*k1*n3*v3**3 + 8*a22**2*a33*b33**2*k1*n1**2*u1*v3 + 4*a22**2*
a33*b33**2*k1*n1**2*v1*v3 + 8*a22**2*a33*b33**2*k1*n1*n3*u3*v3 + 4*a22*a33**3*
b33**3*k1*u1*u2*v1*v2 + 4*a22*a33**3*b33**3*k1*u1*v1**2*v2 + 4*a22*a33**3*b33**3
*k1*u2**2*v1*v3 + 4*a22*a33**3*b33**3*k1*u2*u3*v1*v2 + 4*a22*a33**3*b33**3*k1*u2
*v1**2*v3 + 4*a22*a33**3*b33**3*k1*v1**3*v3 + 7*a22*a33**2*b33**4*k1*u1*v1*v3**2
 + 5*a22*a33**2*b33**4*k1*u2*v2*v3**2 + 7*a22*a33**2*b33**4*k1*u2*v3**3 + 14*a22
*a33**2*b33**4*k1*v1**2*v3**2 + 5*a22*a33**2*b33**4*k1*v1*v3**3 - 2*a22*a33**2*
b33**3*k1*n1*u2*v2*v3 - 4*a22*a33**2*b33**3*k1*n1*u2*v3**2 - 2*a22*a33**2*b33**3
*k1*n1*v2**2*v3 - 4*a22*a33**2*b33**3*k1*n3*v1*v2*v3 - 4*a22*a33**2*b33**3*k1*n3
*v1*v3**2 - 4*a22*a33**2*b33**3*k1*n3*v3**3 - a33**3*b33**4*k1*u1*v1*v3**2 - a33
**3*b33**4*k1*u2*v2*v3**2 - a33**3*b33**4*k1*u2*v3**3 - 2*a33**3*b33**4*k1*v1**2
*v3**2 - a33**3*b33**4*k1*v1*v3**3)/(32*a22**5*(a22**2 - 2*a22*a33 + a33**2))$