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


n3=0


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


        1                1
     - ---*a22*b33*n2 - ---*a33*b33*n2
        2                2
m2=------------------------------------
                 2
              a22  - a22*a33


m3=0


r6=0


     5     2              1     2               1     2
r5=(---*a22 *b33*k1*n2 + ---*a22 *b33*k12*n2 - ---*a22 *b33*n2*q1
     2                    2                     2

     - 5*a22*a33*b33*k1*n2 - a22*a33*b33*k12*n2 + a22*a33*b33*n2*q1

        5     2              1     2               1     2
     + ---*a33 *b33*k1*n2 + ---*a33 *b33*k12*n2 - ---*a33 *b33*n2*q1
        2                    2                     2

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

             3
    - a22*a33 )


     5     2              1     2               1     2
r4=(---*a22 *b33*k1*n1 + ---*a22 *b33*k12*n1 - ---*a22 *b33*n1*q1
     2                    2                     2

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

        5     2              1     2               1     2
     + ---*a33 *b33*k1*n1 + ---*a33 *b33*k12*n1 - ---*a33 *b33*n1*q1
        2                    2                     2

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

             3
    - a22*a33 )


     - k1*n1*n2 + k12*n1*n2
r3=-------------------------
       2                  2
    a22  - 2*a22*a33 + a33


          2            2                                                2
r2=( - a22 *k1*n2 + a22 *n2*q1 + 2*a22*a33*k1*n2 - 2*a22*a33*n2*q1 - a33 *k1*n2

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

         3
    - a33 )


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

                            2            2             2              3
     - 2*a22*a33*n1*q1 + a33 *k1*n1 - a33 *k12*n1 + a33 *n1*q1 - k1*n1

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


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


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

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

         1     2      2    1     2      2    1     2       2    1     2       2
      - ---*b33 *k1*n1  + ---*b33 *k1*n2  + ---*b33 *k12*n1  + ---*b33 *k12*n2 )
         4                 4                 4                  4

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


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


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


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

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

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


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


        3               2                 2                    2
q13=(a22 *k38*n2 - 2*a22 *a33*k38*n2 - a22 *b33*k1*n2 + a22*a33 *k38*n2

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

             3
    - a22*a33 )


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


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

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

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


           - b33*k1*n1*n2
q9=------------------------------
       3        2              2
    a22  - 2*a22 *a33 + a22*a33


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


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

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


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


           - b33*k1*n1*n2
q5=------------------------------
       3        2              2
    a22  - 2*a22 *a33 + a22*a33


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

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

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


     - 2*k1*n2
q3=------------
    a22 - a33


    2*k1*n1*n2 - 2*k12*n1*n2
q2=--------------------------
       2                  2
    a22  - 2*a22*a33 + a33


p56=0


         1     4                 3                   1     2    2
p55=( - ---*a22 *b33*k16*n2 + a22 *a33*b33*k16*n2 - ---*a22 *a33 *b33*k16*n2
         2                                           2

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

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

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


p54=0


         1     4                 3                   1     2    2
p53=( - ---*a22 *b33*k16*n2 + a22 *a33*b33*k16*n2 - ---*a22 *a33 *b33*k16*n2
         2                                           2

         1     2    3          1     2    3                     3
      - ---*a22 *b33 *k1*n2 + ---*a22 *b33 *k12*n2 + a22*a33*b33 *k1*n2
         2                     8

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

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


         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


         1     4                 3                   1     2    2
p48=( - ---*a22 *b33*k16*n2 + a22 *a33*b33*k16*n2 - ---*a22 *a33 *b33*k16*n2
         2                                           2

         1     2    3          1     2    3                     3
      - ---*a22 *b33 *k1*n2 + ---*a22 *b33 *k12*n2 + a22*a33*b33 *k1*n2
         2                     8

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

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


         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


         1     2               1                              2
p45=( - ---*a22 *b33*k38*n2 + ---*a22*a33*b33*k38*n2 - a22*b33 *k1*n2
         2                     2

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

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


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


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


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


p37=0


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

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

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


p35=0


        4               3               1     3                 2    2
p34=(a22 *k16*n2 - 2*a22 *a33*k16*n2 - ---*a22 *b33*k38*n2 + a22 *a33 *k16*n2
                                        2

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

         1         2                           2          1             2
      - ---*a22*a33 *b33*k38*n2 - 2*a22*a33*b33 *k1*n2 + ---*a22*a33*b33 *k12*n2
         2                                                2

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

         2    3
    - a22 *a33 )


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


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


        4               3                 2    2             2    2
p31=(a22 *k16*n2 - 2*a22 *a33*k16*n2 + a22 *a33 *k16*n2 + a22 *b33 *k1*n2

         1     2    2                       2          1             2
      - ---*a22 *b33 *k12*n2 - 2*a22*a33*b33 *k1*n2 + ---*a22*a33*b33 *k12*n2
         4                                             2

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

         2    3
    - a22 *a33 )


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

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


p29=0


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


      k12*n2
p27=-----------
     a22 - a33


p26=0


           2                                            1
p25=( - a22 *k38*n2 + a22*a33*k38*n2 - a22*b33*k1*n2 - ---*a22*b33*k12*n2
                                                        2

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


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


p23=0


     2*k1*n2 - k12*n2
p22=------------------
        a22 - a33


        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


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


        2           1                      2          1     2
     a22 *k16*n1 - ---*a22*b33*k38*n1 + b33 *k1*n1 - ---*b33 *k12*n1
                    2                                 4
p16=-----------------------------------------------------------------
                                3      2
                             a22  - a22 *a33


           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


     - k38*n2
p9=-----------
    a22 - a33


p8= - 4*k1


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


p6=0


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


          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


    2*k1*n2 - k12*n2
p2=------------------
       a22 - a33


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

Relevance for the application:

The following expression INT is a first

integral for the Hamiltonian HAM:

          3   2        3   2        2       2        2       2        2       2
HAM=(4*a22 *u1  + 4*a22 *u2  - 4*a22 *a33*u1  - 4*a22 *a33*u2  + 4*a22 *a33*u3

             2                  2              2                  2   2
      + 4*a22 *b33*u3*v3 + 4*a22 *n1*u1 + 4*a22 *n2*u2 - 4*a22*a33 *u3

                                                                         2   2
      - 4*a22*a33*b33*u3*v3 - 4*a22*a33*n1*u1 - 4*a22*a33*n2*u2 - a22*b33 *v3

                                                   2   2
      - 2*a22*b33*n1*v1 - 2*a22*b33*n2*v2 + a33*b33 *v3  - 2*a33*b33*n1*v1

      - 2*a33*b33*n2*v2)/(4*a22*(a22 - a33))

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

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

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

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

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

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

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

              7                  6          4         6          2
      + 16*a22 *q1*u1*v1 - 48*a22 *a33*k1*u1  - 96*a22 *a33*k1*u1 *u2*u3

              6                2          6             2
      - 48*a22 *a33*k1*u1*u2*v2  + 192*a22 *a33*k1*u1*u3

               6                      6                 2
      + 192*a22 *a33*k1*u1*v1 - 48*a22 *a33*k12*u1*v1*v2

              6           2   2         6           3
      - 48*a22 *a33*k12*u2 *u3  - 48*a22 *a33*k16*u1 *v1

              6           3            6
      - 48*a22 *a33*k16*u1 *v3 - 48*a22 *a33*k16*u1*u2*u3*v1

              6              3         6           3            6              3
      - 48*a22 *a33*k16*u1*v2  - 48*a22 *a33*k16*u2 *v1 - 48*a22 *a33*k16*v1*v2

              6              2            6
      - 48*a22 *a33*k38*u1*u2 *v2 + 96*a22 *a33*k38*u1*u3*v1

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

              6                             6
      - 48*a22 *a33*k38*u2*u3*v2*v3 - 96*a22 *a33*k38*u2*v1*v2

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

              6                         6                      6
      + 32*a22 *b33*k1*u1*u3*v1 + 96*a22 *b33*k1*u1*u3 - 32*a22 *b33*k1*u2*v1*v2

              6          4         6          2   2         6          2
      + 32*a22 *b33*k1*u3  + 32*a22 *b33*k1*v1 *v2  + 96*a22 *b33*k1*v1

              6                             6
      + 16*a22 *b33*k12*u2*u3*v2*v3 - 16*a22 *b33*k16*u1*u2*v1*v2

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

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

              6           3            6                      6          2
      - 16*a22 *b33*k16*v1 *v3 - 16*a22 *b33*q1*u1*u3 - 16*a22 *b33*q1*v1

              6         3         6                  6         2
      + 32*a22 *k1*n1*u1  + 16*a22 *k1*n1*u1 + 32*a22 *k1*n1*u2 *u3

              6                     6         2            6
      - 32*a22 *k1*n1*u2*v1 + 32*a22 *k1*n2*u1 *u2 + 32*a22 *k1*n2*u1*u2*v2

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

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

              6                      6          2            6
      + 32*a22 *k12*n1*u2*v1 - 16*a22 *k12*n2*u1 *u2 - 16*a22 *k12*n2*u1*u2*v2

              6                         6          2            6          2
      + 16*a22 *k12*n2*u1*v1*v2 + 16*a22 *k16*n1*u1 *v2 + 16*a22 *k16*n1*u3 *v1

              6             2         6          2            6             2
      + 16*a22 *k16*n1*u3*v1  + 16*a22 *k16*n2*u1 *v3 + 16*a22 *k16*n2*u1*v2

              6             2         6          3         6
      + 16*a22 *k16*n2*v1*v2  - 16*a22 *k38*n1*u2  - 16*a22 *k38*n1*u2*u3*v1

              6          3         6             2         6
      - 16*a22 *k38*n1*u3  - 16*a22 *k38*n2*u2*u3  - 16*a22 *k38*n2*u2*u3*v2

              6                      6          2            6
      + 16*a22 *k38*n2*u3*v2 - 16*a22 *k38*n2*v1 *v2 + 16*a22 *n1*q1*u1

              6                  5    2      4         5    2      2
      + 16*a22 *n2*q1*u2 + 48*a22 *a33 *k1*u1  + 96*a22 *a33 *k1*u1 *u2*u3

              5    2            2          5    2         2
      + 48*a22 *a33 *k1*u1*u2*v2  - 192*a22 *a33 *k1*u1*u3

               5    2                  5    2             2
      - 192*a22 *a33 *k1*u1*v1 + 48*a22 *a33 *k12*u1*v1*v2

              5    2       2   2         5    2       3
      + 48*a22 *a33 *k12*u2 *u3  + 48*a22 *a33 *k16*u1 *v1

              5    2       3            5    2
      + 48*a22 *a33 *k16*u1 *v3 + 48*a22 *a33 *k16*u1*u2*u3*v1

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

              5    2          3         5    2          2
      + 48*a22 *a33 *k16*v1*v2  + 48*a22 *a33 *k38*u1*u2 *v2

              5    2                      5    2
      - 96*a22 *a33 *k38*u1*u3*v1 - 48*a22 *a33 *k38*u1*u3

              5    2       3            5    2
      + 48*a22 *a33 *k38*u2 *v3 + 48*a22 *a33 *k38*u2*u3*v2*v3

              5    2                      5    2       2         5    2      2
      + 96*a22 *a33 *k38*u2*v1*v2 - 48*a22 *a33 *k38*v1  + 48*a22 *a33 *q1*u1

              5    2                  5
      + 48*a22 *a33 *q1*u1*v1 - 96*a22 *a33*b33*k1*u1*u3*v1

               5                          5
      - 288*a22 *a33*b33*k1*u1*u3 + 96*a22 *a33*b33*k1*u2*v1*v2

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

               5              2         5
      - 288*a22 *a33*b33*k1*v1  - 48*a22 *a33*b33*k12*u2*u3*v2*v3

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

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

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

              5                          5              2
      + 48*a22 *a33*b33*q1*u1*u3 + 48*a22 *a33*b33*q1*v1

              5             3         5                      5             2
      - 64*a22 *a33*k1*n1*u1  - 32*a22 *a33*k1*n1*u1 - 64*a22 *a33*k1*n1*u2 *u3

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

              5                            5             2
      - 64*a22 *a33*k1*n2*u1*u2*v2 + 64*a22 *a33*k1*n2*u2

              5                      5              3
      + 32*a22 *a33*k1*n2*u2 + 32*a22 *a33*k12*n1*u1

              5                             5
      - 32*a22 *a33*k12*n1*u1*u2*v1 + 32*a22 *a33*k12*n1*u1

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

              5              2            5
      + 32*a22 *a33*k12*n2*u1 *u2 + 32*a22 *a33*k12*n2*u1*u2*v2

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

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

              5              2            5                 2
      - 32*a22 *a33*k16*n2*u1 *v3 - 32*a22 *a33*k16*n2*u1*v2

              5                 2         5              3
      - 32*a22 *a33*k16*n2*v1*v2  + 32*a22 *a33*k38*n1*u2

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

              5                 2         5
      + 32*a22 *a33*k38*n2*u2*u3  + 32*a22 *a33*k38*n2*u2*u3*v2

              5                          5              2
      - 32*a22 *a33*k38*n2*u3*v2 + 32*a22 *a33*k38*n2*v1 *v2

              5                      5                     5    2      3
      - 32*a22 *a33*n1*q1*u1 - 32*a22 *a33*n2*q1*u2 + 8*a22 *b33 *k1*u1 *v3

             5    2                        5    2         2
      + 8*a22 *b33 *k1*u1*u2*u3*v1 - 32*a22 *b33 *k1*u1*u3 *v2

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

              5    2         3         5    2                  5    2      3
      + 16*a22 *b33 *k1*u1*v2  - 32*a22 *b33 *k1*u1*v3 + 16*a22 *b33 *k1*u2 *v1

              5    2      3            5    2                 5    2       3
      - 32*a22 *b33 *k1*u3 *v3 - 32*a22 *b33 *k1*v1*v3 - 4*a22 *b33 *k12*u1 *v3

             5    2                        5    2          2
      - 4*a22 *b33 *k12*u1*u2*u3*v1 + 8*a22 *b33 *k12*u1*u3 *v2

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

             5    2          3        5    2       3           5    2       3
      - 4*a22 *b33 *k12*u1*v2  - 4*a22 *b33 *k12*u2 *v1 + 8*a22 *b33 *k12*u3 *v3

             5    2          2        5    2       3        5    2             2
      + 4*a22 *b33 *k12*u3*v2  - 4*a22 *b33 *k12*v1  + 4*a22 *b33 *k16*u1*v1*v3

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

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

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

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

              5                         5
      + 16*a22 *b33*k1*n1*v1*v2 + 40*a22 *b33*k1*n1*v1

              5             2            5
      - 16*a22 *b33*k1*n2*u1 *u3 - 16*a22 *b33*k1*n2*u2*u3*v2

              5                         5             2
      - 16*a22 *b33*k1*n2*u3*v2 + 32*a22 *b33*k1*n2*v1 *v2

              5                     5
      + 40*a22 *b33*k1*n2*v2 - 8*a22 *b33*k12*n1*u1*u3*v2

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

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

             5                      5              2
      + 8*a22 *b33*k12*n1*v1 - 8*a22 *b33*k12*n2*u1 *u3

             5                            5
      - 8*a22 *b33*k12*n2*u1*u3*v3 - 8*a22 *b33*k12*n2*u2*u3*v2

              5              2         5              2
      + 16*a22 *b33*k12*n2*u3  - 16*a22 *b33*k12*n2*v1 *v2

             5                      5
      + 8*a22 *b33*k12*n2*v2 - 8*a22 *b33*k16*n1*u2*v2*v3

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

             5                            5                 2
      - 8*a22 *b33*k16*n2*u3*v2*v3 - 8*a22 *b33*k16*n2*u3*v3

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

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

             5                 2        5              2
      - 8*a22 *b33*k38*n2*u1*v1  - 8*a22 *b33*k38*n2*v1 *v3

             5                 2        5                     5
      - 8*a22 *b33*k38*n2*v1*v2  - 8*a22 *b33*n1*q1*v1 - 8*a22 *b33*n2*q1*v2

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

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

              5                      5       2               4    3      4
      + 16*a22 *k12*n1*n2*u3 - 16*a22 *k12*n2 *u1*v1 - 16*a22 *a33 *k1*u1

              4    3      2               4    3            2
      - 32*a22 *a33 *k1*u1 *u2*u3 - 16*a22 *a33 *k1*u1*u2*v2

              4    3         2         4    3
      + 64*a22 *a33 *k1*u1*u3  + 64*a22 *a33 *k1*u1*v1

              4    3             2         4    3       2   2
      - 16*a22 *a33 *k12*u1*v1*v2  - 16*a22 *a33 *k12*u2 *u3

              4    3       3            4    3       3
      - 16*a22 *a33 *k16*u1 *v1 - 16*a22 *a33 *k16*u1 *v3

              4    3                         4    3          3
      - 16*a22 *a33 *k16*u1*u2*u3*v1 - 16*a22 *a33 *k16*u1*v2

              4    3       3            4    3          3
      - 16*a22 *a33 *k16*u2 *v1 - 16*a22 *a33 *k16*v1*v2

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

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

              4    3                         4    3
      - 16*a22 *a33 *k38*u2*u3*v2*v3 - 32*a22 *a33 *k38*u2*v1*v2

              4    3       2         4    3      2         4    3
      + 16*a22 *a33 *k38*v1  - 16*a22 *a33 *q1*u1  - 16*a22 *a33 *q1*u1*v1

              4    2                          4    2
      + 96*a22 *a33 *b33*k1*u1*u3*v1 + 288*a22 *a33 *b33*k1*u1*u3

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

              4    2          2   2          4    2          2
      + 96*a22 *a33 *b33*k1*v1 *v2  + 288*a22 *a33 *b33*k1*v1

              4    2                             4    2
      + 48*a22 *a33 *b33*k12*u2*u3*v2*v3 - 48*a22 *a33 *b33*k16*u1*u2*v1*v2

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

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

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

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

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

              4    2         2            4    2
      + 32*a22 *a33 *k1*n2*u1 *u2 + 32*a22 *a33 *k1*n2*u1*u2*v2

              4    2         2         4    2                  4    2          3
      - 32*a22 *a33 *k1*n2*u2  - 16*a22 *a33 *k1*n2*u2 - 16*a22 *a33 *k12*n1*u1

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

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

              4    2          2            4    2
      - 16*a22 *a33 *k12*n2*u1 *u2 - 16*a22 *a33 *k12*n2*u1*u2*v2

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

              4    2          2            4    2             2
      + 16*a22 *a33 *k16*n1*u3 *v1 + 16*a22 *a33 *k16*n1*u3*v1

              4    2          2            4    2             2
      + 16*a22 *a33 *k16*n2*u1 *v3 + 16*a22 *a33 *k16*n2*u1*v2

              4    2             2         4    2          3
      + 16*a22 *a33 *k16*n2*v1*v2  - 16*a22 *a33 *k38*n1*u2

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

              4    2             2         4    2
      - 16*a22 *a33 *k38*n2*u2*u3  - 16*a22 *a33 *k38*n2*u2*u3*v2

              4    2                      4    2          2
      + 16*a22 *a33 *k38*n2*u3*v2 - 16*a22 *a33 *k38*n2*v1 *v2

              4    2                  4    2
      + 16*a22 *a33 *n1*q1*u1 + 16*a22 *a33 *n2*q1*u2

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

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

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

              4        2                  4        2      3
      + 96*a22 *a33*b33 *k1*u1*v3 - 48*a22 *a33*b33 *k1*u2 *v1

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

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

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

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

              4        2       3            4        2       3
      + 12*a22 *a33*b33 *k12*u2 *v1 - 24*a22 *a33*b33 *k12*u3 *v3

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

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

              4        2          3         4        2       2   2
      - 12*a22 *a33*b33 *k16*u2*v3  - 24*a22 *a33*b33 *k16*v1 *v3

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

              4        2                  4
      - 12*a22 *a33*b33 *q1*v1*v3 + 32*a22 *a33*b33*k1*n1*u1*u3*v2

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

              4                             4
      - 32*a22 *a33*b33*k1*n1*v1*v2 - 80*a22 *a33*b33*k1*n1*v1

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

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

              4                          4
      - 80*a22 *a33*b33*k1*n2*v2 + 16*a22 *a33*b33*k12*n1*u1*u3*v2

              4                              4                  3
      + 32*a22 *a33*b33*k12*n1*u1*v2 + 16*a22 *a33*b33*k12*n1*u2

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

              4                           4                  2
      - 16*a22 *a33*b33*k12*n1*v1 + 16*a22 *a33*b33*k12*n2*u1 *u3

              4                                 4
      + 16*a22 *a33*b33*k12*n2*u1*u3*v3 + 16*a22 *a33*b33*k12*n2*u2*u3*v2

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

              4                           4
      - 16*a22 *a33*b33*k12*n2*v2 + 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
      + 16*a22 *a33*b33*k16*n2*u3*v2*v3 + 16*a22 *a33*b33*k16*n2*u3*v3

              4                     2         4
      + 16*a22 *a33*b33*k16*n2*v2*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                     2         4                  2
      + 16*a22 *a33*b33*k38*n2*u1*v1  + 16*a22 *a33*b33*k38*n2*v1 *v3

              4                     2         4
      + 16*a22 *a33*b33*k38*n2*v1*v2  + 16*a22 *a33*b33*n1*q1*v1

              4                          4          2
      + 16*a22 *a33*b33*n2*q1*v2 + 16*a22 *a33*k1*n1 *u1*v1

              4                            4
      - 32*a22 *a33*k1*n1*n2*u1*u2 + 16*a22 *a33*k1*n1*n2*u3

              4          2               4           2
      - 16*a22 *a33*k1*n2 *u1*v1 - 16*a22 *a33*k12*n1 *u1*v1

              4                             4
      + 32*a22 *a33*k12*n1*n2*u1*u2 - 16*a22 *a33*k12*n1*n2*u3

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

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

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

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

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

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

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

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

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

              4    2            2        4    2         2
      + 16*a22 *b33 *k1*n1*u3*v1  + 8*a22 *b33 *k1*n2*u1 *v3

              4    2            2         4    2         2
      + 16*a22 *b33 *k1*n2*u1*v2  - 16*a22 *b33 *k1*n2*v1 *v3

              4    2            2        4    2
      + 16*a22 *b33 *k1*n2*v1*v2  - 8*a22 *b33 *k12*n1*u2*v1*v3

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

             4    2                     4    2             2
      + 8*a22 *b33 *k12*n1*v2*v3 - 4*a22 *b33 *k12*n2*u1*v2

             4    2                     4    2          2
      - 8*a22 *b33 *k12*n2*u3*v3 - 8*a22 *b33 *k12*n2*v1 *v3

             4    2             2         4          2   2
      - 4*a22 *b33 *k12*n2*v1*v2  + 16*a22 *b33*k1*n1 *v1

              4                            4
      - 16*a22 *b33*k1*n1*n2*u2*u3 - 16*a22 *b33*k1*n1*n2*u3*v1

              4          2   2         4           2
      - 16*a22 *b33*k1*n2 *v1  - 16*a22 *b33*k12*n1 *u1*u3

              4           2   2         4      3            4      2
      - 16*a22 *b33*k12*n1 *v1  - 16*a22 *k1*n1 *u1 - 16*a22 *k1*n1 *n2*u2

              4       3            4       2
      + 16*a22 *k12*n1 *u1 + 16*a22 *k12*n1 *n2*u2

              3    3                         3    3
      - 32*a22 *a33 *b33*k1*u1*u3*v1 - 96*a22 *a33 *b33*k1*u1*u3

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

              3    3          2   2         3    3          2
      - 32*a22 *a33 *b33*k1*v1 *v2  - 96*a22 *a33 *b33*k1*v1

              3    3                             3    3
      - 16*a22 *a33 *b33*k12*u2*u3*v2*v3 + 16*a22 *a33 *b33*k16*u1*u2*v1*v2

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

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

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

              3    3          2         3    2    2      3
      + 16*a22 *a33 *b33*q1*v1  + 24*a22 *a33 *b33 *k1*u1 *v3

              3    2    2                        3    2    2         2
      + 24*a22 *a33 *b33 *k1*u1*u2*u3*v1 - 96*a22 *a33 *b33 *k1*u1*u3 *v2

              3    2    2         3         3    2    2
      - 96*a22 *a33 *b33 *k1*u1*v1  + 48*a22 *a33 *b33 *k1*u1*v1*v3

              3    2    2         3         3    2    2
      + 48*a22 *a33 *b33 *k1*u1*v2  - 96*a22 *a33 *b33 *k1*u1*v3

              3    2    2      3            3    2    2      3
      + 48*a22 *a33 *b33 *k1*u2 *v1 - 96*a22 *a33 *b33 *k1*u3 *v3

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

              3    2    2                         3    2    2          2
      - 12*a22 *a33 *b33 *k12*u1*u2*u3*v1 + 24*a22 *a33 *b33 *k12*u1*u3 *v2

              3    2    2          3         3    2    2
      + 24*a22 *a33 *b33 *k12*u1*v1  - 12*a22 *a33 *b33 *k12*u1*v1*v3

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

              3    2    2       3            3    2    2          2
      + 24*a22 *a33 *b33 *k12*u3 *v3 + 12*a22 *a33 *b33 *k12*u3*v2

              3    2    2       3         3    2    2             2
      - 12*a22 *a33 *b33 *k12*v1  + 12*a22 *a33 *b33 *k16*u1*v1*v3

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

              3    2    2       2   2         3    2    2          3
      + 24*a22 *a33 *b33 *k16*v1 *v3  + 12*a22 *a33 *b33 *k16*v1*v3

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

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

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

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

              3    2                            3    2
      - 16*a22 *a33 *b33*k1*n2*u2*u3*v2 - 16*a22 *a33 *b33*k1*n2*u3*v2

              3    2             2            3    2
      + 32*a22 *a33 *b33*k1*n2*v1 *v2 + 40*a22 *a33 *b33*k1*n2*v2

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

             3    2              3        3    2              2
      - 8*a22 *a33 *b33*k12*n1*u2  - 8*a22 *a33 *b33*k12*n1*u2 *v3

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

             3    2              2           3    2
      - 8*a22 *a33 *b33*k12*n2*u1 *u3 - 8*a22 *a33 *b33*k12*n2*u1*u3*v3

             3    2                             3    2              2
      - 8*a22 *a33 *b33*k12*n2*u2*u3*v2 + 16*a22 *a33 *b33*k12*n2*u3

              3    2              2           3    2
      - 16*a22 *a33 *b33*k12*n2*v1 *v2 + 8*a22 *a33 *b33*k12*n2*v2

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

             3    2              2           3    2
      - 8*a22 *a33 *b33*k16*n1*v2 *v3 - 8*a22 *a33 *b33*k16*n2*u3*v2*v3

             3    2                 2        3    2                 2
      - 8*a22 *a33 *b33*k16*n2*u3*v3  - 8*a22 *a33 *b33*k16*n2*v2*v3

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

             3    2              2           3    2                 2
      - 8*a22 *a33 *b33*k38*n1*u3 *v1 - 8*a22 *a33 *b33*k38*n2*u1*v1

             3    2              2           3    2                 2
      - 8*a22 *a33 *b33*k38*n2*v1 *v3 - 8*a22 *a33 *b33*k38*n2*v1*v2

             3    2                     3    2
      - 8*a22 *a33 *b33*n1*q1*v1 - 8*a22 *a33 *b33*n2*q1*v2

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

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

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

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

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

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

              3        3          2            3        3       3
      - 12*a22 *a33*b33 *k12*u2*v1 *v3 - 12*a22 *a33*b33 *k12*v1 *v3

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

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

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

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

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

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

             3        2             2         3        2
      + 8*a22 *a33*b33 *k12*n2*u1*v2  + 16*a22 *a33*b33 *k12*n2*u3*v3

              3        2          2           3        2             2
      + 16*a22 *a33*b33 *k12*n2*v1 *v3 + 8*a22 *a33*b33 *k12*n2*v1*v2

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

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

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

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

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

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

             3    4       2   2      3    4          3
      - 2*a22 *b33 *k12*v1 *v3  - a22 *b33 *k12*v1*v3

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

             3    3         2           3    3
      - 8*a22 *b33 *k1*n1*v2 *v3 - 8*a22 *b33 *k1*n2*u3*v2*v3

             3    3            2        3    3            2
      - 8*a22 *b33 *k1*n2*u3*v3  - 4*a22 *b33 *k1*n2*v2*v3

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

             3    3          2           3    3
      + 2*a22 *b33 *k12*n1*v2 *v3 + 2*a22 *b33 *k12*n2*u3*v2*v3

             3    3             2        3    3             2
      + 2*a22 *b33 *k12*n2*u3*v3  + 4*a22 *b33 *k12*n2*v2*v3

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

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

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

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

             3          2              3           3
      + 8*a22 *b33*k1*n1 *n2*v2 - 8*a22 *b33*k12*n1 *v1

             3           2              2    3    2      3
      - 8*a22 *b33*k12*n1 *n2*v2 - 8*a22 *a33 *b33 *k1*u1 *v3

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

              2    2    3          2            2    2    3       3
      + 12*a22 *a33 *b33 *k12*u2*v1 *v3 + 12*a22 *a33 *b33 *k12*v1 *v3

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

              2    2    2         2            2    2    2            2
      + 16*a22 *a33 *b33 *k1*n1*u3 *v1 + 16*a22 *a33 *b33 *k1*n1*u3*v1

             2    2    2         2            2    2    2            2
      + 8*a22 *a33 *b33 *k1*n2*u1 *v3 + 16*a22 *a33 *b33 *k1*n2*u1*v2

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

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

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

             2    2    2             2        2    2    2
      - 4*a22 *a33 *b33 *k12*n2*u1*v2  - 8*a22 *a33 *b33 *k12*n2*u3*v3

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

             2        4            2        2        4            2
      - 9*a22 *a33*b33 *k1*u1*v1*v3  - 6*a22 *a33*b33 *k1*u2*v2*v3

             2        4         3         2        4      2   2
      - 9*a22 *a33*b33 *k1*u2*v3  - 18*a22 *a33*b33 *k1*v1 *v3

             2        4         3        2        4             2
      - 6*a22 *a33*b33 *k1*v1*v3  + 3*a22 *a33*b33 *k12*u1*v1*v3

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

             2        4       2   2        2        4          3
      + 6*a22 *a33*b33 *k12*v1 *v3  + 3*a22 *a33*b33 *k12*v1*v3

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

              2        3         2            2        3
      + 16*a22 *a33*b33 *k1*n1*v2 *v3 + 16*a22 *a33*b33 *k1*n2*u3*v2*v3

              2        3            2        2        3            2
      + 16*a22 *a33*b33 *k1*n2*u3*v3  + 8*a22 *a33*b33 *k1*n2*v2*v3

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

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

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

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

             2        2      2              2        2       2
      - 4*a22 *a33*b33 *k1*n2 *v1*v3 - 8*a22 *a33*b33 *k12*n1 *u1*v3

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

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

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

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

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

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

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

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

                 3    3       3               2    4            2
      - 4*a22*a33 *b33 *k12*v1 *v3 + 9*a22*a33 *b33 *k1*u1*v1*v3

                 2    4            2            2    4         3
      + 6*a22*a33 *b33 *k1*u2*v2*v3  + 9*a22*a33 *b33 *k1*u2*v3

                  2    4      2   2            2    4         3
      + 18*a22*a33 *b33 *k1*v1 *v3  + 6*a22*a33 *b33 *k1*v1*v3

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

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

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

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

                 2    3                           2    3            2
      - 8*a22*a33 *b33 *k1*n2*u3*v2*v3 - 8*a22*a33 *b33 *k1*n2*u3*v3

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

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

                 2    3                            2    3             2
      + 2*a22*a33 *b33 *k12*n2*u3*v2*v3 + 2*a22*a33 *b33 *k12*n2*u3*v3

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

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

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

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

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

INT=(16*a22**7*k1*u1**4 + 32*a22**7*k1*u1**2*u2*u3 + 16*a22**7*k1*u1*u2*v2**2 - 
64*a22**7*k1*u1*u3**2 - 64*a22**7*k1*u1*v1 + 16*a22**7*k12*u1*v1*v2**2 + 16*a22
**7*k12*u2**2*u3**2 + 16*a22**7*k16*u1**3*v1 + 16*a22**7*k16*u1**3*v3 + 16*a22**
7*k16*u1*u2*u3*v1 + 16*a22**7*k16*u1*v2**3 + 16*a22**7*k16*u2**3*v1 + 16*a22**7*
k16*v1*v2**3 + 16*a22**7*k38*u1*u2**2*v2 - 32*a22**7*k38*u1*u3*v1 - 16*a22**7*
k38*u1*u3 + 16*a22**7*k38*u2**3*v3 + 16*a22**7*k38*u2*u3*v2*v3 + 32*a22**7*k38*
u2*v1*v2 - 16*a22**7*k38*v1**2 + 16*a22**7*q1*u1**2 + 16*a22**7*q1*u1*v1 - 48*
a22**6*a33*k1*u1**4 - 96*a22**6*a33*k1*u1**2*u2*u3 - 48*a22**6*a33*k1*u1*u2*v2**
2 + 192*a22**6*a33*k1*u1*u3**2 + 192*a22**6*a33*k1*u1*v1 - 48*a22**6*a33*k12*u1*
v1*v2**2 - 48*a22**6*a33*k12*u2**2*u3**2 - 48*a22**6*a33*k16*u1**3*v1 - 48*a22**
6*a33*k16*u1**3*v3 - 48*a22**6*a33*k16*u1*u2*u3*v1 - 48*a22**6*a33*k16*u1*v2**3 
- 48*a22**6*a33*k16*u2**3*v1 - 48*a22**6*a33*k16*v1*v2**3 - 48*a22**6*a33*k38*u1
*u2**2*v2 + 96*a22**6*a33*k38*u1*u3*v1 + 48*a22**6*a33*k38*u1*u3 - 48*a22**6*a33
*k38*u2**3*v3 - 48*a22**6*a33*k38*u2*u3*v2*v3 - 96*a22**6*a33*k38*u2*v1*v2 + 48*
a22**6*a33*k38*v1**2 - 48*a22**6*a33*q1*u1**2 - 48*a22**6*a33*q1*u1*v1 + 32*a22
**6*b33*k1*u1*u3*v1 + 96*a22**6*b33*k1*u1*u3 - 32*a22**6*b33*k1*u2*v1*v2 + 32*
a22**6*b33*k1*u3**4 + 32*a22**6*b33*k1*v1**2*v2**2 + 96*a22**6*b33*k1*v1**2 + 16
*a22**6*b33*k12*u2*u3*v2*v3 - 16*a22**6*b33*k16*u1*u2*v1*v2 - 16*a22**6*b33*k16*
u1*v1**2*v2 - 16*a22**6*b33*k16*u2**2*v1*v3 - 16*a22**6*b33*k16*u2*u3*v1*v2 - 16
*a22**6*b33*k16*u2*v1**2*v3 - 16*a22**6*b33*k16*v1**3*v3 - 16*a22**6*b33*q1*u1*
u3 - 16*a22**6*b33*q1*v1**2 + 32*a22**6*k1*n1*u1**3 + 16*a22**6*k1*n1*u1 + 32*
a22**6*k1*n1*u2**2*u3 - 32*a22**6*k1*n1*u2*v1 + 32*a22**6*k1*n2*u1**2*u2 + 32*
a22**6*k1*n2*u1*u2*v2 - 32*a22**6*k1*n2*u2**2 - 16*a22**6*k1*n2*u2 - 16*a22**6*
k12*n1*u1**3 + 16*a22**6*k12*n1*u1*u2*v1 - 16*a22**6*k12*n1*u1 - 16*a22**6*k12*
n1*u2**2*u3 + 32*a22**6*k12*n1*u2*v1 - 16*a22**6*k12*n2*u1**2*u2 - 16*a22**6*k12
*n2*u1*u2*v2 + 16*a22**6*k12*n2*u1*v1*v2 + 16*a22**6*k16*n1*u1**2*v2 + 16*a22**6
*k16*n1*u3**2*v1 + 16*a22**6*k16*n1*u3*v1**2 + 16*a22**6*k16*n2*u1**2*v3 + 16*
a22**6*k16*n2*u1*v2**2 + 16*a22**6*k16*n2*v1*v2**2 - 16*a22**6*k38*n1*u2**3 - 16
*a22**6*k38*n1*u2*u3*v1 - 16*a22**6*k38*n1*u3**3 - 16*a22**6*k38*n2*u2*u3**2 - 
16*a22**6*k38*n2*u2*u3*v2 + 16*a22**6*k38*n2*u3*v2 - 16*a22**6*k38*n2*v1**2*v2 +
 16*a22**6*n1*q1*u1 + 16*a22**6*n2*q1*u2 + 48*a22**5*a33**2*k1*u1**4 + 96*a22**5
*a33**2*k1*u1**2*u2*u3 + 48*a22**5*a33**2*k1*u1*u2*v2**2 - 192*a22**5*a33**2*k1*
u1*u3**2 - 192*a22**5*a33**2*k1*u1*v1 + 48*a22**5*a33**2*k12*u1*v1*v2**2 + 48*
a22**5*a33**2*k12*u2**2*u3**2 + 48*a22**5*a33**2*k16*u1**3*v1 + 48*a22**5*a33**2
*k16*u1**3*v3 + 48*a22**5*a33**2*k16*u1*u2*u3*v1 + 48*a22**5*a33**2*k16*u1*v2**3
 + 48*a22**5*a33**2*k16*u2**3*v1 + 48*a22**5*a33**2*k16*v1*v2**3 + 48*a22**5*a33
**2*k38*u1*u2**2*v2 - 96*a22**5*a33**2*k38*u1*u3*v1 - 48*a22**5*a33**2*k38*u1*u3
 + 48*a22**5*a33**2*k38*u2**3*v3 + 48*a22**5*a33**2*k38*u2*u3*v2*v3 + 96*a22**5*
a33**2*k38*u2*v1*v2 - 48*a22**5*a33**2*k38*v1**2 + 48*a22**5*a33**2*q1*u1**2 + 
48*a22**5*a33**2*q1*u1*v1 - 96*a22**5*a33*b33*k1*u1*u3*v1 - 288*a22**5*a33*b33*
k1*u1*u3 + 96*a22**5*a33*b33*k1*u2*v1*v2 - 96*a22**5*a33*b33*k1*u3**4 - 96*a22**
5*a33*b33*k1*v1**2*v2**2 - 288*a22**5*a33*b33*k1*v1**2 - 48*a22**5*a33*b33*k12*
u2*u3*v2*v3 + 48*a22**5*a33*b33*k16*u1*u2*v1*v2 + 48*a22**5*a33*b33*k16*u1*v1**2
*v2 + 48*a22**5*a33*b33*k16*u2**2*v1*v3 + 48*a22**5*a33*b33*k16*u2*u3*v1*v2 + 48
*a22**5*a33*b33*k16*u2*v1**2*v3 + 48*a22**5*a33*b33*k16*v1**3*v3 + 48*a22**5*a33
*b33*q1*u1*u3 + 48*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*n2*u1**2*u2 - 64*a22**5*a33*k1*n2*u1*u2*v2 + 64*a22**5*a33*k1*n2*
u2**2 + 32*a22**5*a33*k1*n2*u2 + 32*a22**5*a33*k12*n1*u1**3 - 32*a22**5*a33*k12*
n1*u1*u2*v1 + 32*a22**5*a33*k12*n1*u1 + 32*a22**5*a33*k12*n1*u2**2*u3 - 64*a22**
5*a33*k12*n1*u2*v1 + 32*a22**5*a33*k12*n2*u1**2*u2 + 32*a22**5*a33*k12*n2*u1*u2*
v2 - 32*a22**5*a33*k12*n2*u1*v1*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*n2*u1**2
*v3 - 32*a22**5*a33*k16*n2*u1*v2**2 - 32*a22**5*a33*k16*n2*v1*v2**2 + 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*n2*u2*u3**2 + 32*a22**5*a33*k38*n2*u2*u3*v2 - 32*a22**5*a33*
k38*n2*u3*v2 + 32*a22**5*a33*k38*n2*v1**2*v2 - 32*a22**5*a33*n1*q1*u1 - 32*a22**
5*a33*n2*q1*u2 + 8*a22**5*b33**2*k1*u1**3*v3 + 8*a22**5*b33**2*k1*u1*u2*u3*v1 - 
32*a22**5*b33**2*k1*u1*u3**2*v2 - 32*a22**5*b33**2*k1*u1*v1**3 + 16*a22**5*b33**
2*k1*u1*v1*v3 + 16*a22**5*b33**2*k1*u1*v2**3 - 32*a22**5*b33**2*k1*u1*v3 + 16*
a22**5*b33**2*k1*u2**3*v1 - 32*a22**5*b33**2*k1*u3**3*v3 - 32*a22**5*b33**2*k1*
v1*v3 - 4*a22**5*b33**2*k12*u1**3*v3 - 4*a22**5*b33**2*k12*u1*u2*u3*v1 + 8*a22**
5*b33**2*k12*u1*u3**2*v2 + 8*a22**5*b33**2*k12*u1*v1**3 - 4*a22**5*b33**2*k12*u1
*v1*v3 - 4*a22**5*b33**2*k12*u1*v2**3 - 4*a22**5*b33**2*k12*u2**3*v1 + 8*a22**5*
b33**2*k12*u3**3*v3 + 4*a22**5*b33**2*k12*u3*v2**2 - 4*a22**5*b33**2*k12*v1**3 +
 4*a22**5*b33**2*k16*u1*v1*v3**2 + 4*a22**5*b33**2*k16*u2*v2*v3**2 + 4*a22**5*
b33**2*k16*u2*v3**3 + 8*a22**5*b33**2*k16*v1**2*v3**2 + 4*a22**5*b33**2*k16*v1*
v3**3 + 4*a22**5*b33**2*q1*u1*v3 + 4*a22**5*b33**2*q1*v1*v3 - 16*a22**5*b33*k1*
n1*u1*u3*v2 - 16*a22**5*b33*k1*n1*u2**3 + 32*a22**5*b33*k1*n1*u2*u3*v1 + 16*a22
**5*b33*k1*n1*v1*v2 + 40*a22**5*b33*k1*n1*v1 - 16*a22**5*b33*k1*n2*u1**2*u3 - 16
*a22**5*b33*k1*n2*u2*u3*v2 - 16*a22**5*b33*k1*n2*u3*v2 + 32*a22**5*b33*k1*n2*v1
**2*v2 + 40*a22**5*b33*k1*n2*v2 - 8*a22**5*b33*k12*n1*u1*u3*v2 - 16*a22**5*b33*
k12*n1*u1*v2 - 8*a22**5*b33*k12*n1*u2**3 - 8*a22**5*b33*k12*n1*u2**2*v3 - 16*a22
**5*b33*k12*n1*u2*u3*v1 + 8*a22**5*b33*k12*n1*v1 - 8*a22**5*b33*k12*n2*u1**2*u3 
- 8*a22**5*b33*k12*n2*u1*u3*v3 - 8*a22**5*b33*k12*n2*u2*u3*v2 + 16*a22**5*b33*
k12*n2*u3**2 - 16*a22**5*b33*k12*n2*v1**2*v2 + 8*a22**5*b33*k12*n2*v2 - 8*a22**5
*b33*k16*n1*u2*v2*v3 - 8*a22**5*b33*k16*n1*u2*v3**2 - 8*a22**5*b33*k16*n1*v2**2*
v3 - 8*a22**5*b33*k16*n2*u3*v2*v3 - 8*a22**5*b33*k16*n2*u3*v3**2 - 8*a22**5*b33*
k16*n2*v2*v3**2 - 8*a22**5*b33*k38*n1*u2*v1*v3 - 8*a22**5*b33*k38*n1*u2*v2**2 - 
8*a22**5*b33*k38*n1*u3**2*v1 - 8*a22**5*b33*k38*n2*u1*v1**2 - 8*a22**5*b33*k38*
n2*v1**2*v3 - 8*a22**5*b33*k38*n2*v1*v2**2 - 8*a22**5*b33*n1*q1*v1 - 8*a22**5*
b33*n2*q1*v2 - 16*a22**5*k1*n1**2*u1*v1 + 32*a22**5*k1*n1*n2*u1*u2 - 16*a22**5*
k1*n1*n2*u3 + 16*a22**5*k1*n2**2*u1*v1 + 16*a22**5*k12*n1**2*u1*v1 - 32*a22**5*
k12*n1*n2*u1*u2 + 16*a22**5*k12*n1*n2*u3 - 16*a22**5*k12*n2**2*u1*v1 - 16*a22**4
*a33**3*k1*u1**4 - 32*a22**4*a33**3*k1*u1**2*u2*u3 - 16*a22**4*a33**3*k1*u1*u2*
v2**2 + 64*a22**4*a33**3*k1*u1*u3**2 + 64*a22**4*a33**3*k1*u1*v1 - 16*a22**4*a33
**3*k12*u1*v1*v2**2 - 16*a22**4*a33**3*k12*u2**2*u3**2 - 16*a22**4*a33**3*k16*u1
**3*v1 - 16*a22**4*a33**3*k16*u1**3*v3 - 16*a22**4*a33**3*k16*u1*u2*u3*v1 - 16*
a22**4*a33**3*k16*u1*v2**3 - 16*a22**4*a33**3*k16*u2**3*v1 - 16*a22**4*a33**3*
k16*v1*v2**3 - 16*a22**4*a33**3*k38*u1*u2**2*v2 + 32*a22**4*a33**3*k38*u1*u3*v1 
+ 16*a22**4*a33**3*k38*u1*u3 - 16*a22**4*a33**3*k38*u2**3*v3 - 16*a22**4*a33**3*
k38*u2*u3*v2*v3 - 32*a22**4*a33**3*k38*u2*v1*v2 + 16*a22**4*a33**3*k38*v1**2 - 
16*a22**4*a33**3*q1*u1**2 - 16*a22**4*a33**3*q1*u1*v1 + 96*a22**4*a33**2*b33*k1*
u1*u3*v1 + 288*a22**4*a33**2*b33*k1*u1*u3 - 96*a22**4*a33**2*b33*k1*u2*v1*v2 + 
96*a22**4*a33**2*b33*k1*u3**4 + 96*a22**4*a33**2*b33*k1*v1**2*v2**2 + 288*a22**4
*a33**2*b33*k1*v1**2 + 48*a22**4*a33**2*b33*k12*u2*u3*v2*v3 - 48*a22**4*a33**2*
b33*k16*u1*u2*v1*v2 - 48*a22**4*a33**2*b33*k16*u1*v1**2*v2 - 48*a22**4*a33**2*
b33*k16*u2**2*v1*v3 - 48*a22**4*a33**2*b33*k16*u2*u3*v1*v2 - 48*a22**4*a33**2*
b33*k16*u2*v1**2*v3 - 48*a22**4*a33**2*b33*k16*v1**3*v3 - 48*a22**4*a33**2*b33*
q1*u1*u3 - 48*a22**4*a33**2*b33*q1*v1**2 + 32*a22**4*a33**2*k1*n1*u1**3 + 16*a22
**4*a33**2*k1*n1*u1 + 32*a22**4*a33**2*k1*n1*u2**2*u3 - 32*a22**4*a33**2*k1*n1*
u2*v1 + 32*a22**4*a33**2*k1*n2*u1**2*u2 + 32*a22**4*a33**2*k1*n2*u1*u2*v2 - 32*
a22**4*a33**2*k1*n2*u2**2 - 16*a22**4*a33**2*k1*n2*u2 - 16*a22**4*a33**2*k12*n1*
u1**3 + 16*a22**4*a33**2*k12*n1*u1*u2*v1 - 16*a22**4*a33**2*k12*n1*u1 - 16*a22**
4*a33**2*k12*n1*u2**2*u3 + 32*a22**4*a33**2*k12*n1*u2*v1 - 16*a22**4*a33**2*k12*
n2*u1**2*u2 - 16*a22**4*a33**2*k12*n2*u1*u2*v2 + 16*a22**4*a33**2*k12*n2*u1*v1*
v2 + 16*a22**4*a33**2*k16*n1*u1**2*v2 + 16*a22**4*a33**2*k16*n1*u3**2*v1 + 16*
a22**4*a33**2*k16*n1*u3*v1**2 + 16*a22**4*a33**2*k16*n2*u1**2*v3 + 16*a22**4*a33
**2*k16*n2*u1*v2**2 + 16*a22**4*a33**2*k16*n2*v1*v2**2 - 16*a22**4*a33**2*k38*n1
*u2**3 - 16*a22**4*a33**2*k38*n1*u2*u3*v1 - 16*a22**4*a33**2*k38*n1*u3**3 - 16*
a22**4*a33**2*k38*n2*u2*u3**2 - 16*a22**4*a33**2*k38*n2*u2*u3*v2 + 16*a22**4*a33
**2*k38*n2*u3*v2 - 16*a22**4*a33**2*k38*n2*v1**2*v2 + 16*a22**4*a33**2*n1*q1*u1 
+ 16*a22**4*a33**2*n2*q1*u2 - 24*a22**4*a33*b33**2*k1*u1**3*v3 - 24*a22**4*a33*
b33**2*k1*u1*u2*u3*v1 + 96*a22**4*a33*b33**2*k1*u1*u3**2*v2 + 96*a22**4*a33*b33
**2*k1*u1*v1**3 - 48*a22**4*a33*b33**2*k1*u1*v1*v3 - 48*a22**4*a33*b33**2*k1*u1*
v2**3 + 96*a22**4*a33*b33**2*k1*u1*v3 - 48*a22**4*a33*b33**2*k1*u2**3*v1 + 96*
a22**4*a33*b33**2*k1*u3**3*v3 + 96*a22**4*a33*b33**2*k1*v1*v3 + 12*a22**4*a33*
b33**2*k12*u1**3*v3 + 12*a22**4*a33*b33**2*k12*u1*u2*u3*v1 - 24*a22**4*a33*b33**
2*k12*u1*u3**2*v2 - 24*a22**4*a33*b33**2*k12*u1*v1**3 + 12*a22**4*a33*b33**2*k12
*u1*v1*v3 + 12*a22**4*a33*b33**2*k12*u1*v2**3 + 12*a22**4*a33*b33**2*k12*u2**3*
v1 - 24*a22**4*a33*b33**2*k12*u3**3*v3 - 12*a22**4*a33*b33**2*k12*u3*v2**2 + 12*
a22**4*a33*b33**2*k12*v1**3 - 12*a22**4*a33*b33**2*k16*u1*v1*v3**2 - 12*a22**4*
a33*b33**2*k16*u2*v2*v3**2 - 12*a22**4*a33*b33**2*k16*u2*v3**3 - 24*a22**4*a33*
b33**2*k16*v1**2*v3**2 - 12*a22**4*a33*b33**2*k16*v1*v3**3 - 12*a22**4*a33*b33**
2*q1*u1*v3 - 12*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 + 32*a22**4*a33*b33*k1*n2*u1**2*u3
 + 32*a22**4*a33*b33*k1*n2*u2*u3*v2 + 32*a22**4*a33*b33*k1*n2*u3*v2 - 64*a22**4*
a33*b33*k1*n2*v1**2*v2 - 80*a22**4*a33*b33*k1*n2*v2 + 16*a22**4*a33*b33*k12*n1*
u1*u3*v2 + 32*a22**4*a33*b33*k12*n1*u1*v2 + 16*a22**4*a33*b33*k12*n1*u2**3 + 16*
a22**4*a33*b33*k12*n1*u2**2*v3 + 32*a22**4*a33*b33*k12*n1*u2*u3*v1 - 16*a22**4*
a33*b33*k12*n1*v1 + 16*a22**4*a33*b33*k12*n2*u1**2*u3 + 16*a22**4*a33*b33*k12*n2
*u1*u3*v3 + 16*a22**4*a33*b33*k12*n2*u2*u3*v2 - 32*a22**4*a33*b33*k12*n2*u3**2 +
 32*a22**4*a33*b33*k12*n2*v1**2*v2 - 16*a22**4*a33*b33*k12*n2*v2 + 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 + 16*a22**4*a33*b33*k16*n2*u3*v2*v3 + 16*a22**4*a33*b33*k16*n2*u3*
v3**2 + 16*a22**4*a33*b33*k16*n2*v2*v3**2 + 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*k38*n2*u1*v1**2 + 16*a22**4*a33*b33*k38*n2*v1**2*v3 + 16*a22**4*a33*
b33*k38*n2*v1*v2**2 + 16*a22**4*a33*b33*n1*q1*v1 + 16*a22**4*a33*b33*n2*q1*v2 + 
16*a22**4*a33*k1*n1**2*u1*v1 - 32*a22**4*a33*k1*n1*n2*u1*u2 + 16*a22**4*a33*k1*
n1*n2*u3 - 16*a22**4*a33*k1*n2**2*u1*v1 - 16*a22**4*a33*k12*n1**2*u1*v1 + 32*a22
**4*a33*k12*n1*n2*u1*u2 - 16*a22**4*a33*k12*n1*n2*u3 + 16*a22**4*a33*k12*n2**2*
u1*v1 - 8*a22**4*b33**3*k1*u1*u2*v1*v2 + 8*a22**4*b33**3*k1*u1*u3*v3**2 - 8*a22
**4*b33**3*k1*u1*v1**2*v2 - 8*a22**4*b33**3*k1*u2**2*v1*v3 - 8*a22**4*b33**3*k1*
u2*u3*v1*v2 - 8*a22**4*b33**3*k1*u2*v1**2*v3 - 8*a22**4*b33**3*k1*v1**3*v3 + 8*
a22**4*b33**3*k1*v1*v2**2*v3 + 4*a22**4*b33**3*k12*u1*u2*v1*v2 + 4*a22**4*b33**3
*k12*u1*v1**2*v2 + 4*a22**4*b33**3*k12*u2**2*v1*v3 + 4*a22**4*b33**3*k12*u2*u3*
v1*v2 + 4*a22**4*b33**3*k12*u2*v1**2*v3 + 4*a22**4*b33**3*k12*v1**3*v3 + 8*a22**
4*b33**2*k1*n1*u1**2*v2 - 16*a22**4*b33**2*k1*n1*u2*v1*v3 + 16*a22**4*b33**2*k1*
n1*u3**2*v1 + 16*a22**4*b33**2*k1*n1*u3*v1**2 + 8*a22**4*b33**2*k1*n2*u1**2*v3 +
 16*a22**4*b33**2*k1*n2*u1*v2**2 - 16*a22**4*b33**2*k1*n2*v1**2*v3 + 16*a22**4*
b33**2*k1*n2*v1*v2**2 - 8*a22**4*b33**2*k12*n1*u2*v1*v3 - 4*a22**4*b33**2*k12*n1
*u3**2*v1 - 4*a22**4*b33**2*k12*n1*u3*v1**2 + 8*a22**4*b33**2*k12*n1*v2*v3 - 4*
a22**4*b33**2*k12*n2*u1*v2**2 - 8*a22**4*b33**2*k12*n2*u3*v3 - 8*a22**4*b33**2*
k12*n2*v1**2*v3 - 4*a22**4*b33**2*k12*n2*v1*v2**2 + 16*a22**4*b33*k1*n1**2*v1**2
 - 16*a22**4*b33*k1*n1*n2*u2*u3 - 16*a22**4*b33*k1*n1*n2*u3*v1 - 16*a22**4*b33*
k1*n2**2*v1**2 - 16*a22**4*b33*k12*n1**2*u1*u3 - 16*a22**4*b33*k12*n1**2*v1**2 -
 16*a22**4*k1*n1**3*u1 - 16*a22**4*k1*n1**2*n2*u2 + 16*a22**4*k12*n1**3*u1 + 16*
a22**4*k12*n1**2*n2*u2 - 32*a22**3*a33**3*b33*k1*u1*u3*v1 - 96*a22**3*a33**3*b33
*k1*u1*u3 + 32*a22**3*a33**3*b33*k1*u2*v1*v2 - 32*a22**3*a33**3*b33*k1*u3**4 - 
32*a22**3*a33**3*b33*k1*v1**2*v2**2 - 96*a22**3*a33**3*b33*k1*v1**2 - 16*a22**3*
a33**3*b33*k12*u2*u3*v2*v3 + 16*a22**3*a33**3*b33*k16*u1*u2*v1*v2 + 16*a22**3*
a33**3*b33*k16*u1*v1**2*v2 + 16*a22**3*a33**3*b33*k16*u2**2*v1*v3 + 16*a22**3*
a33**3*b33*k16*u2*u3*v1*v2 + 16*a22**3*a33**3*b33*k16*u2*v1**2*v3 + 16*a22**3*
a33**3*b33*k16*v1**3*v3 + 16*a22**3*a33**3*b33*q1*u1*u3 + 16*a22**3*a33**3*b33*
q1*v1**2 + 24*a22**3*a33**2*b33**2*k1*u1**3*v3 + 24*a22**3*a33**2*b33**2*k1*u1*
u2*u3*v1 - 96*a22**3*a33**2*b33**2*k1*u1*u3**2*v2 - 96*a22**3*a33**2*b33**2*k1*
u1*v1**3 + 48*a22**3*a33**2*b33**2*k1*u1*v1*v3 + 48*a22**3*a33**2*b33**2*k1*u1*
v2**3 - 96*a22**3*a33**2*b33**2*k1*u1*v3 + 48*a22**3*a33**2*b33**2*k1*u2**3*v1 -
 96*a22**3*a33**2*b33**2*k1*u3**3*v3 - 96*a22**3*a33**2*b33**2*k1*v1*v3 - 12*a22
**3*a33**2*b33**2*k12*u1**3*v3 - 12*a22**3*a33**2*b33**2*k12*u1*u2*u3*v1 + 24*
a22**3*a33**2*b33**2*k12*u1*u3**2*v2 + 24*a22**3*a33**2*b33**2*k12*u1*v1**3 - 12
*a22**3*a33**2*b33**2*k12*u1*v1*v3 - 12*a22**3*a33**2*b33**2*k12*u1*v2**3 - 12*
a22**3*a33**2*b33**2*k12*u2**3*v1 + 24*a22**3*a33**2*b33**2*k12*u3**3*v3 + 12*
a22**3*a33**2*b33**2*k12*u3*v2**2 - 12*a22**3*a33**2*b33**2*k12*v1**3 + 12*a22**
3*a33**2*b33**2*k16*u1*v1*v3**2 + 12*a22**3*a33**2*b33**2*k16*u2*v2*v3**2 + 12*
a22**3*a33**2*b33**2*k16*u2*v3**3 + 24*a22**3*a33**2*b33**2*k16*v1**2*v3**2 + 12
*a22**3*a33**2*b33**2*k16*v1*v3**3 + 12*a22**3*a33**2*b33**2*q1*u1*v3 + 12*a22**
3*a33**2*b33**2*q1*v1*v3 - 16*a22**3*a33**2*b33*k1*n1*u1*u3*v2 - 16*a22**3*a33**
2*b33*k1*n1*u2**3 + 32*a22**3*a33**2*b33*k1*n1*u2*u3*v1 + 16*a22**3*a33**2*b33*
k1*n1*v1*v2 + 40*a22**3*a33**2*b33*k1*n1*v1 - 16*a22**3*a33**2*b33*k1*n2*u1**2*
u3 - 16*a22**3*a33**2*b33*k1*n2*u2*u3*v2 - 16*a22**3*a33**2*b33*k1*n2*u3*v2 + 32
*a22**3*a33**2*b33*k1*n2*v1**2*v2 + 40*a22**3*a33**2*b33*k1*n2*v2 - 8*a22**3*a33
**2*b33*k12*n1*u1*u3*v2 - 16*a22**3*a33**2*b33*k12*n1*u1*v2 - 8*a22**3*a33**2*
b33*k12*n1*u2**3 - 8*a22**3*a33**2*b33*k12*n1*u2**2*v3 - 16*a22**3*a33**2*b33*
k12*n1*u2*u3*v1 + 8*a22**3*a33**2*b33*k12*n1*v1 - 8*a22**3*a33**2*b33*k12*n2*u1
**2*u3 - 8*a22**3*a33**2*b33*k12*n2*u1*u3*v3 - 8*a22**3*a33**2*b33*k12*n2*u2*u3*
v2 + 16*a22**3*a33**2*b33*k12*n2*u3**2 - 16*a22**3*a33**2*b33*k12*n2*v1**2*v2 + 
8*a22**3*a33**2*b33*k12*n2*v2 - 8*a22**3*a33**2*b33*k16*n1*u2*v2*v3 - 8*a22**3*
a33**2*b33*k16*n1*u2*v3**2 - 8*a22**3*a33**2*b33*k16*n1*v2**2*v3 - 8*a22**3*a33
**2*b33*k16*n2*u3*v2*v3 - 8*a22**3*a33**2*b33*k16*n2*u3*v3**2 - 8*a22**3*a33**2*
b33*k16*n2*v2*v3**2 - 8*a22**3*a33**2*b33*k38*n1*u2*v1*v3 - 8*a22**3*a33**2*b33*
k38*n1*u2*v2**2 - 8*a22**3*a33**2*b33*k38*n1*u3**2*v1 - 8*a22**3*a33**2*b33*k38*
n2*u1*v1**2 - 8*a22**3*a33**2*b33*k38*n2*v1**2*v3 - 8*a22**3*a33**2*b33*k38*n2*
v1*v2**2 - 8*a22**3*a33**2*b33*n1*q1*v1 - 8*a22**3*a33**2*b33*n2*q1*v2 + 24*a22
**3*a33*b33**3*k1*u1*u2*v1*v2 - 24*a22**3*a33*b33**3*k1*u1*u3*v3**2 + 24*a22**3*
a33*b33**3*k1*u1*v1**2*v2 + 24*a22**3*a33*b33**3*k1*u2**2*v1*v3 + 24*a22**3*a33*
b33**3*k1*u2*u3*v1*v2 + 24*a22**3*a33*b33**3*k1*u2*v1**2*v3 + 24*a22**3*a33*b33
**3*k1*v1**3*v3 - 24*a22**3*a33*b33**3*k1*v1*v2**2*v3 - 12*a22**3*a33*b33**3*k12
*u1*u2*v1*v2 - 12*a22**3*a33*b33**3*k12*u1*v1**2*v2 - 12*a22**3*a33*b33**3*k12*
u2**2*v1*v3 - 12*a22**3*a33*b33**3*k12*u2*u3*v1*v2 - 12*a22**3*a33*b33**3*k12*u2
*v1**2*v3 - 12*a22**3*a33*b33**3*k12*v1**3*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 - 16*a22**3*a33*b33**2*k1*n2*u1**2*v3 - 32
*a22**3*a33*b33**2*k1*n2*u1*v2**2 + 32*a22**3*a33*b33**2*k1*n2*v1**2*v3 - 32*a22
**3*a33*b33**2*k1*n2*v1*v2**2 + 16*a22**3*a33*b33**2*k12*n1*u2*v1*v3 + 8*a22**3*
a33*b33**2*k12*n1*u3**2*v1 + 8*a22**3*a33*b33**2*k12*n1*u3*v1**2 - 16*a22**3*a33
*b33**2*k12*n1*v2*v3 + 8*a22**3*a33*b33**2*k12*n2*u1*v2**2 + 16*a22**3*a33*b33**
2*k12*n2*u3*v3 + 16*a22**3*a33*b33**2*k12*n2*v1**2*v3 + 8*a22**3*a33*b33**2*k12*
n2*v1*v2**2 - 16*a22**3*a33*b33*k1*n1**2*v1**2 + 16*a22**3*a33*b33*k1*n1*n2*u2*
u3 + 16*a22**3*a33*b33*k1*n1*n2*u3*v1 + 16*a22**3*a33*b33*k1*n2**2*v1**2 + 16*
a22**3*a33*b33*k12*n1**2*u1*u3 + 16*a22**3*a33*b33*k12*n1**2*v1**2 + 3*a22**3*
b33**4*k1*u1*v1*v3**2 + 2*a22**3*b33**4*k1*u2*v2*v3**2 + 3*a22**3*b33**4*k1*u2*
v3**3 + 6*a22**3*b33**4*k1*v1**2*v3**2 + 2*a22**3*b33**4*k1*v1*v3**3 - a22**3*
b33**4*k12*u1*v1*v3**2 - a22**3*b33**4*k12*u2*v2*v3**2 - a22**3*b33**4*k12*u2*v3
**3 - 2*a22**3*b33**4*k12*v1**2*v3**2 - a22**3*b33**4*k12*v1*v3**3 - 8*a22**3*
b33**3*k1*n1*u2*v2*v3 - 4*a22**3*b33**3*k1*n1*u2*v3**2 - 8*a22**3*b33**3*k1*n1*
v2**2*v3 - 8*a22**3*b33**3*k1*n2*u3*v2*v3 - 8*a22**3*b33**3*k1*n2*u3*v3**2 - 4*
a22**3*b33**3*k1*n2*v2*v3**2 + 2*a22**3*b33**3*k12*n1*u2*v2*v3 + 4*a22**3*b33**3
*k12*n1*u2*v3**2 + 2*a22**3*b33**3*k12*n1*v2**2*v3 + 2*a22**3*b33**3*k12*n2*u3*
v2*v3 + 2*a22**3*b33**3*k12*n2*u3*v3**2 + 4*a22**3*b33**3*k12*n2*v2*v3**2 - 4*
a22**3*b33**2*k1*n1**2*v1*v3 + 8*a22**3*b33**2*k1*n1*n2*u2*v3 + 4*a22**3*b33**2*
k1*n2**2*v1*v3 + 8*a22**3*b33**2*k12*n1**2*u1*v3 + 4*a22**3*b33**2*k12*n1**2*v1*
v3 + 8*a22**3*b33**2*k12*n1*n2*u2*v3 + 4*a22**3*b33**2*k12*n2**2*v1*v3 + 8*a22**
3*b33*k1*n1**3*v1 + 8*a22**3*b33*k1*n1**2*n2*v2 - 8*a22**3*b33*k12*n1**3*v1 - 8*
a22**3*b33*k12*n1**2*n2*v2 - 8*a22**2*a33**3*b33**2*k1*u1**3*v3 - 8*a22**2*a33**
3*b33**2*k1*u1*u2*u3*v1 + 32*a22**2*a33**3*b33**2*k1*u1*u3**2*v2 + 32*a22**2*a33
**3*b33**2*k1*u1*v1**3 - 16*a22**2*a33**3*b33**2*k1*u1*v1*v3 - 16*a22**2*a33**3*
b33**2*k1*u1*v2**3 + 32*a22**2*a33**3*b33**2*k1*u1*v3 - 16*a22**2*a33**3*b33**2*
k1*u2**3*v1 + 32*a22**2*a33**3*b33**2*k1*u3**3*v3 + 32*a22**2*a33**3*b33**2*k1*
v1*v3 + 4*a22**2*a33**3*b33**2*k12*u1**3*v3 + 4*a22**2*a33**3*b33**2*k12*u1*u2*
u3*v1 - 8*a22**2*a33**3*b33**2*k12*u1*u3**2*v2 - 8*a22**2*a33**3*b33**2*k12*u1*
v1**3 + 4*a22**2*a33**3*b33**2*k12*u1*v1*v3 + 4*a22**2*a33**3*b33**2*k12*u1*v2**
3 + 4*a22**2*a33**3*b33**2*k12*u2**3*v1 - 8*a22**2*a33**3*b33**2*k12*u3**3*v3 - 
4*a22**2*a33**3*b33**2*k12*u3*v2**2 + 4*a22**2*a33**3*b33**2*k12*v1**3 - 4*a22**
2*a33**3*b33**2*k16*u1*v1*v3**2 - 4*a22**2*a33**3*b33**2*k16*u2*v2*v3**2 - 4*a22
**2*a33**3*b33**2*k16*u2*v3**3 - 8*a22**2*a33**3*b33**2*k16*v1**2*v3**2 - 4*a22
**2*a33**3*b33**2*k16*v1*v3**3 - 4*a22**2*a33**3*b33**2*q1*u1*v3 - 4*a22**2*a33
**3*b33**2*q1*v1*v3 - 24*a22**2*a33**2*b33**3*k1*u1*u2*v1*v2 + 24*a22**2*a33**2*
b33**3*k1*u1*u3*v3**2 - 24*a22**2*a33**2*b33**3*k1*u1*v1**2*v2 - 24*a22**2*a33**
2*b33**3*k1*u2**2*v1*v3 - 24*a22**2*a33**2*b33**3*k1*u2*u3*v1*v2 - 24*a22**2*a33
**2*b33**3*k1*u2*v1**2*v3 - 24*a22**2*a33**2*b33**3*k1*v1**3*v3 + 24*a22**2*a33
**2*b33**3*k1*v1*v2**2*v3 + 12*a22**2*a33**2*b33**3*k12*u1*u2*v1*v2 + 12*a22**2*
a33**2*b33**3*k12*u1*v1**2*v2 + 12*a22**2*a33**2*b33**3*k12*u2**2*v1*v3 + 12*a22
**2*a33**2*b33**3*k12*u2*u3*v1*v2 + 12*a22**2*a33**2*b33**3*k12*u2*v1**2*v3 + 12
*a22**2*a33**2*b33**3*k12*v1**3*v3 + 8*a22**2*a33**2*b33**2*k1*n1*u1**2*v2 - 16*
a22**2*a33**2*b33**2*k1*n1*u2*v1*v3 + 16*a22**2*a33**2*b33**2*k1*n1*u3**2*v1 + 
16*a22**2*a33**2*b33**2*k1*n1*u3*v1**2 + 8*a22**2*a33**2*b33**2*k1*n2*u1**2*v3 +
 16*a22**2*a33**2*b33**2*k1*n2*u1*v2**2 - 16*a22**2*a33**2*b33**2*k1*n2*v1**2*v3
 + 16*a22**2*a33**2*b33**2*k1*n2*v1*v2**2 - 8*a22**2*a33**2*b33**2*k12*n1*u2*v1*
v3 - 4*a22**2*a33**2*b33**2*k12*n1*u3**2*v1 - 4*a22**2*a33**2*b33**2*k12*n1*u3*
v1**2 + 8*a22**2*a33**2*b33**2*k12*n1*v2*v3 - 4*a22**2*a33**2*b33**2*k12*n2*u1*
v2**2 - 8*a22**2*a33**2*b33**2*k12*n2*u3*v3 - 8*a22**2*a33**2*b33**2*k12*n2*v1**
2*v3 - 4*a22**2*a33**2*b33**2*k12*n2*v1*v2**2 - 9*a22**2*a33*b33**4*k1*u1*v1*v3
**2 - 6*a22**2*a33*b33**4*k1*u2*v2*v3**2 - 9*a22**2*a33*b33**4*k1*u2*v3**3 - 18*
a22**2*a33*b33**4*k1*v1**2*v3**2 - 6*a22**2*a33*b33**4*k1*v1*v3**3 + 3*a22**2*
a33*b33**4*k12*u1*v1*v3**2 + 3*a22**2*a33*b33**4*k12*u2*v2*v3**2 + 3*a22**2*a33*
b33**4*k12*u2*v3**3 + 6*a22**2*a33*b33**4*k12*v1**2*v3**2 + 3*a22**2*a33*b33**4*
k12*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 + 16*a22**2*a33*b33**3*k1*n2*u3*
v2*v3 + 16*a22**2*a33*b33**3*k1*n2*u3*v3**2 + 8*a22**2*a33*b33**3*k1*n2*v2*v3**2
 - 4*a22**2*a33*b33**3*k12*n1*u2*v2*v3 - 8*a22**2*a33*b33**3*k12*n1*u2*v3**2 - 4
*a22**2*a33*b33**3*k12*n1*v2**2*v3 - 4*a22**2*a33*b33**3*k12*n2*u3*v2*v3 - 4*a22
**2*a33*b33**3*k12*n2*u3*v3**2 - 8*a22**2*a33*b33**3*k12*n2*v2*v3**2 + 4*a22**2*
a33*b33**2*k1*n1**2*v1*v3 - 8*a22**2*a33*b33**2*k1*n1*n2*u2*v3 - 4*a22**2*a33*
b33**2*k1*n2**2*v1*v3 - 8*a22**2*a33*b33**2*k12*n1**2*u1*v3 - 4*a22**2*a33*b33**
2*k12*n1**2*v1*v3 - 8*a22**2*a33*b33**2*k12*n1*n2*u2*v3 - 4*a22**2*a33*b33**2*
k12*n2**2*v1*v3 + 8*a22*a33**3*b33**3*k1*u1*u2*v1*v2 - 8*a22*a33**3*b33**3*k1*u1
*u3*v3**2 + 8*a22*a33**3*b33**3*k1*u1*v1**2*v2 + 8*a22*a33**3*b33**3*k1*u2**2*v1
*v3 + 8*a22*a33**3*b33**3*k1*u2*u3*v1*v2 + 8*a22*a33**3*b33**3*k1*u2*v1**2*v3 + 
8*a22*a33**3*b33**3*k1*v1**3*v3 - 8*a22*a33**3*b33**3*k1*v1*v2**2*v3 - 4*a22*a33
**3*b33**3*k12*u1*u2*v1*v2 - 4*a22*a33**3*b33**3*k12*u1*v1**2*v2 - 4*a22*a33**3*
b33**3*k12*u2**2*v1*v3 - 4*a22*a33**3*b33**3*k12*u2*u3*v1*v2 - 4*a22*a33**3*b33
**3*k12*u2*v1**2*v3 - 4*a22*a33**3*b33**3*k12*v1**3*v3 + 9*a22*a33**2*b33**4*k1*
u1*v1*v3**2 + 6*a22*a33**2*b33**4*k1*u2*v2*v3**2 + 9*a22*a33**2*b33**4*k1*u2*v3
**3 + 18*a22*a33**2*b33**4*k1*v1**2*v3**2 + 6*a22*a33**2*b33**4*k1*v1*v3**3 - 3*
a22*a33**2*b33**4*k12*u1*v1*v3**2 - 3*a22*a33**2*b33**4*k12*u2*v2*v3**2 - 3*a22*
a33**2*b33**4*k12*u2*v3**3 - 6*a22*a33**2*b33**4*k12*v1**2*v3**2 - 3*a22*a33**2*
b33**4*k12*v1*v3**3 - 8*a22*a33**2*b33**3*k1*n1*u2*v2*v3 - 4*a22*a33**2*b33**3*
k1*n1*u2*v3**2 - 8*a22*a33**2*b33**3*k1*n1*v2**2*v3 - 8*a22*a33**2*b33**3*k1*n2*
u3*v2*v3 - 8*a22*a33**2*b33**3*k1*n2*u3*v3**2 - 4*a22*a33**2*b33**3*k1*n2*v2*v3
**2 + 2*a22*a33**2*b33**3*k12*n1*u2*v2*v3 + 4*a22*a33**2*b33**3*k12*n1*u2*v3**2 
+ 2*a22*a33**2*b33**3*k12*n1*v2**2*v3 + 2*a22*a33**2*b33**3*k12*n2*u3*v2*v3 + 2*
a22*a33**2*b33**3*k12*n2*u3*v3**2 + 4*a22*a33**2*b33**3*k12*n2*v2*v3**2 - 3*a33
**3*b33**4*k1*u1*v1*v3**2 - 2*a33**3*b33**4*k1*u2*v2*v3**2 - 3*a33**3*b33**4*k1*
u2*v3**3 - 6*a33**3*b33**4*k1*v1**2*v3**2 - 2*a33**3*b33**4*k1*v1*v3**3 + a33**3
*b33**4*k12*u1*v1*v3**2 + a33**3*b33**4*k12*u2*v2*v3**2 + a33**3*b33**4*k12*u2*
v3**3 + 2*a33**3*b33**4*k12*v1**2*v3**2 + a33**3*b33**4*k12*v1*v3**3)/(16*a22**4
*(a22**3 - 3*a22**2*a33 + 3*a22*a33**2 - a33**3))$