Solution 8 to problem over


Expressions | Parameters | Inequalities | Relevance | Back to problem over

Expressions

The solution is given through the following expressions:

                                                       2
r10=( - 4*a33*b13*c13*c23*m2*n1*n3*r446 - 4*a33*b13*c23 *m2*n2*n3*r446

                     2   2                    3
      + 8*a33*c13*c23 *n1 *n3*r446 + 8*a33*c23 *n1*n2*n3*r446

             2    2      2           2           2
      - 2*b13 *c13 *n1*n3 *r446 + b13 *c13*c23*n1 *n2*r446

           2           3             2              2           2    2   3
      + b13 *c13*c23*n2 *r446 - 2*b13 *c13*c23*n2*n3 *r446 - b13 *c23 *n1 *r446

           2    2      2             2    3    2        2    3    2
      - b13 *c23 *n1*n2 *r446)/(4*a33 *b13 *c13  + 4*a33 *b13 *c23 )


                                                     2      2
r11=( - 2*a33*b13*c13*c23*m2*n1*n2*r446 + a33*b13*c23 *m2*n1 *r446

                   2      2                     2   2
      - a33*b13*c23 *m2*n2 *r446 + 4*a33*c13*c23 *n1 *n2*r446

                 3      2           2    2
      + 4*a33*c23 *n1*n2 *r446 - b13 *c13 *n1*n2*n3*r446

           2           2                2    3    2        2    3    2
      - b13 *c13*c23*n2 *n3*r446)/(2*a33 *b13 *c13  + 2*a33 *b13 *c23 )


                               2                               2
r12=( - 2*a33*b13*c13*c23*m2*n1 *r446 + 2*a33*b13*c13*c23*m2*n2 *r446

                     2                              2   3
      - 4*a33*b13*c23 *m2*n1*n2*r446 + 4*a33*c13*c23 *n1 *r446

                     2      2                 3   2
      - 4*a33*c13*c23 *n1*n2 *r446 + 4*a33*c23 *n1 *n2*r446

                 3   3           2    2   2              2    2   2
      - 4*a33*c23 *n2 *r446 - b13 *c13 *n1 *n3*r446 + b13 *c13 *n2 *n3*r446

             2                            2    2   2
      - 2*b13 *c13*c23*n1*n2*n3*r446 + b13 *c23 *n1 *n3*r446

           2    2   2                2    3    2        2    3    2
      + b13 *c23 *n2 *n3*r446)/(4*a33 *b13 *c13  + 4*a33 *b13 *c23 )


        2   2              2   2
     c23 *n1 *n3*r446 + c23 *n2 *n3*r446
r13=-------------------------------------
               2    2            2    2
      2*a33*b13 *c13  + 2*a33*b13 *c23


        2   2              2   3
     c23 *n1 *n2*r446 + c23 *n2 *r446
r14=-----------------------------------
              2    2            2    2
     2*a33*b13 *c13  + 2*a33*b13 *c23


        2   3           2      2
     c23 *n1 *r446 + c23 *n1*n2 *r446
r15=-----------------------------------
              2    2            2    2
     2*a33*b13 *c13  + 2*a33*b13 *c23


          2    2    2    2   2             2    2    4   2
r20=(4*a33 *b13 *c13 *c23 *m2 *r446 + 4*a33 *b13 *c23 *m2 *r446

              2        2    3                    2            4
      - 16*a33 *b13*c13 *c23 *m2*n1*r446 - 16*a33 *b13*c13*c23 *m2*n2*r446

              2    2    4   2              2        5
      + 16*a33 *c13 *c23 *n1 *r446 + 32*a33 *c13*c23 *n1*n2*r446

              2    6   2                 3    3
      + 16*a33 *c23 *n2 *r446 + 4*a33*b13 *c13 *c23*m2*n3*r446

                 3        3                       2    3    2
      + 4*a33*b13 *c13*c23 *m2*n3*r446 - 8*a33*b13 *c13 *c23 *n1*n3*r446

                 2    2    3                        2        4
      - 8*a33*b13 *c13 *c23 *n2*n3*r446 - 16*a33*b13 *c13*c23 *n1*n3*r446

                  2    5                 4    4   2
      - 16*a33*b13 *c23 *n2*n3*r446 + b13 *c13 *n3 *r446

           4    2    2   2             4    2    2   2
      + b13 *c13 *c23 *n1 *r446 - 3*b13 *c13 *c23 *n2 *r446

           4    2    2   2             4        3
      + b13 *c13 *c23 *n3 *r446 + 4*b13 *c13*c23 *n1*n2*r446

           4    4   2           4    4   2             2    4    2    2
      - b13 *c23 *n1 *r446 - b13 *c23 *n2 *r446)/(8*a33 *b13 *c13 *c23

           2    4    4
    + 8*a33 *b13 *c23 )


                  2                               3
r21=(2*a33*b13*c13 *c23*m2*n1*r446 + 2*a33*b13*c23 *m2*n1*r446

                 2    2   2                     3
      - 4*a33*c13 *c23 *n1 *r446 - 8*a33*c13*c23 *n1*n2*r446

                 4   2           2    3                 2        2
      + 4*a33*c23 *n1 *r446 + b13 *c13 *n1*n3*r446 - b13 *c13*c23 *n1*n3*r446

             2    3                   2    2    2            2    2    3
      - 2*b13 *c23 *n2*n3*r446)/(4*a33 *b13 *c13 *c23 + 4*a33 *b13 *c23 )


       2          2
     n1 *r446 - n2 *r446
r22=---------------------
                2
           4*a33


                     2                               3
r23=( - 2*a33*b13*c13 *c23*m2*n2*r446 - 2*a33*b13*c23 *m2*n2*r446

                 2    2                           3   2
      + 4*a33*c13 *c23 *n1*n2*r446 + 8*a33*c13*c23 *n2 *r446

                 4                 2    3
      - 4*a33*c23 *n1*n2*r446 - b13 *c13 *n2*n3*r446

             2    2                       2        2
      - 2*b13 *c13 *c23*n1*n3*r446 - 3*b13 *c13*c23 *n2*n3*r446)/(

        2    2    2            2    2    3
   4*a33 *b13 *c13 *c23 + 4*a33 *b13 *c23 )


      - n1*n2*r446
r24=---------------
             2
        2*a33


                                                    2
r26=( - 2*a33*b13*c13*c23*m2*n1*r446 - 2*a33*b13*c23 *m2*n2*r446

                     2   2                 3                 2    2
      + 4*a33*c13*c23 *n1 *r446 + 4*a33*c23 *n1*n2*r446 - b13 *c13 *n1*n3*r446

           2                               3    2            3    2
      - b13 *c13*c23*n2*n3*r446)/(2*a33*b13 *c13  + 2*a33*b13 *c23 )


                  2                  2             2
      - c13*c23*n1 *r446 + c13*c23*n2 *r446 - 2*c23 *n1*n2*r446
r27=------------------------------------------------------------
                               2                2
                  2*a33*b13*c13  + 2*a33*b13*c23


                               2   2           2   2
     2*c13*c23*n1*n2*r446 - c23 *n1 *r446 + c23 *n2 *r446
r28=------------------------------------------------------
                            2                2
               2*a33*b13*c13  + 2*a33*b13*c23


                   2           2
       - c13*c23*n2 *r446 + c23 *n1*n2*r446
r210=---------------------------------------
                      2              2
           a33*b13*c13  + a33*b13*c23


                                                     2
r212=( - 2*a33*b13*c13*c23*m2*n2*r446 + 2*a33*b13*c23 *m2*n1*r446

                      2                       3   2           2    2
       + 4*a33*c13*c23 *n1*n2*r446 + 4*a33*c23 *n2 *r446 - b13 *c13 *n2*n3*r446

            2                           2    2                       3    2
       - b13 *c13*c23*n1*n3*r446 - 2*b13 *c23 *n2*n3*r446)/(2*a33*b13 *c13

               3    2
    + 2*a33*b13 *c23 )


r213=0


            2   2           2   2
       - c23 *n1 *r446 - c23 *n2 *r446
r214=----------------------------------
              2    2        2    2
         2*b13 *c13  + 2*b13 *c23


                                 2   2
       - c13*c23*n1*n2*r446 + c23 *n1 *r446
r215=---------------------------------------
                      2              2
           a33*b13*c13  + a33*b13*c23


                                                  2
r216=(2*a33*b13*c13*c23*m2*n2*r446 - 2*a33*b13*c23 *m2*n1*r446

                      2                       3   2           2    2
       - 4*a33*c13*c23 *n1*n2*r446 - 4*a33*c23 *n2 *r446 + b13 *c13 *n2*n3*r446

            2                           2    2                       3    2
       + b13 *c13*c23*n1*n3*r446 + 2*b13 *c23 *n2*n3*r446)/(2*a33*b13 *c13

               3    2
    + 2*a33*b13 *c23 )


r217=0


r218=0


r219=0


            2   2           2   2
       - c23 *n1 *r446 - c23 *n2 *r446
r220=----------------------------------
              2    2        2    2
         2*b13 *c13  + 2*b13 *c23


     b13*c13*n2*r446 - b13*c23*n1*r446
r30=-----------------------------------
                     2
                4*a33 *c23


                                      2              2
     2*a33*b13*c23*m2*r446 + 4*a33*c23 *n1*r446 + b13 *c13*n3*r446
r31=---------------------------------------------------------------
                                 2
                            4*a33 *b13*c23


             2        2    2                2        4
r32=( - 8*a33 *b13*c13 *c23 *m2*r446 - 8*a33 *b13*c23 *m2*r446

              2    2    3                 2        4
      + 16*a33 *c13 *c23 *n1*r446 + 16*a33 *c13*c23 *n2*r446

                 2    3                        2        3
      - 4*a33*b13 *c13 *c23*n3*r446 - 4*a33*b13 *c13*c23 *n3*r446

           4    3              4        2                2    3    2
      + b13 *c13 *n2*r446 + b13 *c13*c23 *n2*r446)/(4*a33 *b13 *c13 *c23

           2    3    3
    + 4*a33 *b13 *c23 )


        2                      2
     c13 *c23*n1*r446 + c13*c23 *n2*r446
r33=-------------------------------------
                    2              2
         a33*b13*c13  + a33*b13*c23


                                              2                    3
r34=(2*a33*b13*c13*c23*m2*r446 - 4*a33*c13*c23 *n1*r446 - 8*a33*c23 *n2*r446

           2    2                2        2
      + b13 *c13 *n3*r446)/(4*a33 *b13*c23 )


     b13*c13*n1*r446 - b13*c23*n2*r446
r35=-----------------------------------
                     2
                4*a33 *c23


               2              3
      - c13*c23 *n1*r446 - c23 *n2*r446
r36=------------------------------------
                   2              2
        a33*b13*c13  + a33*b13*c23


             2        2                    2        3
r37=( - 8*a33 *b13*c13 *c23*m2*r446 - 8*a33 *b13*c23 *m2*r446

              2    2    2                 2        3
      + 16*a33 *c13 *c23 *n1*r446 + 16*a33 *c13*c23 *n2*r446

                 2    3                    2        2
      - 4*a33*b13 *c13 *n3*r446 - 4*a33*b13 *c13*c23 *n3*r446

           4    2              4    2                2    3    2
      - b13 *c13 *n1*r446 - b13 *c23 *n1*r446)/(4*a33 *b13 *c13

           2    3    2
    + 4*a33 *b13 *c23 )


        2                      2
     c13 *c23*n1*r446 + c13*c23 *n2*r446
r38=-------------------------------------
                    2              2
         a33*b13*c13  + a33*b13*c23


               2              3
      - c13*c23 *n1*r446 - c23 *n2*r446
r39=------------------------------------
                   2              2
        a33*b13*c13  + a33*b13*c23


r310=0


       - n2*r446
r311=------------
        2*a33


               2                3
      2*c13*c23 *n1*r446 + 2*c23 *n2*r446
r312=-------------------------------------
                2    2      2    2
             b13 *c13  + b13 *c23


       - n1*r446
r313=------------
        2*a33


r314=0


               2                3
      2*c13*c23 *n1*r446 + 2*c23 *n2*r446
r315=-------------------------------------
                2    2      2    2
             b13 *c13  + b13 *c23


r316=0


r317=0


r318=0


r319=0


      n2*r446
r320=---------
        a33


r321=0


         2                                  2
      c13 *n2*r446 + c13*c23*n1*r446 + 2*c23 *n2*r446
r322=-------------------------------------------------
                           2            2
                  2*a33*c13  + 2*a33*c23


                   2                            3
r323=(2*a33*b13*c13 *c23*m2*r446 + 2*a33*b13*c23 *m2*r446

                  2    2                        3                    4
       - 4*a33*c13 *c23 *n1*r446 - 8*a33*c13*c23 *n2*r446 + 4*a33*c23 *n1*r446

            2    3              2        2                    2    2
       + b13 *c13 *n3*r446 + b13 *c13*c23 *n3*r446)/(2*a33*b13 *c13 *c23

               2    3
    + 2*a33*b13 *c23 )


      n1*r446
r324=---------
       2*a33


                           2
      c13*c23*n1*r446 + c23 *n2*r446
r325=--------------------------------
                  2            2
         2*a33*c13  + 2*a33*c23


r326=0


                              2
       - c13*c23*n2*r446 + c23 *n1*r446
r327=-----------------------------------
                    2          2
             b13*c13  + b13*c23


                              2
       - c13*c23*n1*r446 - c23 *n2*r446
r328=-----------------------------------
                    2          2
             b13*c13  + b13*c23


r329=0


                           2
      c13*c23*n2*r446 - c23 *n1*r446
r330=--------------------------------
                  2          2
           b13*c13  + b13*c23


r331=0


r332=0


r333=0


r334=0


      n1*r446
r335=---------
        a33


                      2                            3
r336=( - 2*a33*b13*c13 *c23*m2*r446 - 2*a33*b13*c23 *m2*r446

                  2    2                        3                    4
       + 4*a33*c13 *c23 *n1*r446 + 8*a33*c13*c23 *n2*r446 - 4*a33*c23 *n1*r446

            2    3              2        2                    2    2
       - b13 *c13 *n3*r446 - b13 *c13*c23 *n3*r446)/(2*a33*b13 *c13 *c23

               2    3
    + 2*a33*b13 *c23 )


         2
      c13 *n1*r446 + c13*c23*n2*r446
r337=--------------------------------
                  2            2
         2*a33*c13  + 2*a33*c23


r338=0


      n2*r446
r339=---------
       2*a33


           2                                2
      2*c13 *n1*r446 + c13*c23*n2*r446 + c23 *n1*r446
r340=-------------------------------------------------
                           2            2
                  2*a33*c13  + 2*a33*c23


r341=0


                           2
      c13*c23*n1*r446 + c23 *n2*r446
r342=--------------------------------
                  2          2
           b13*c13  + b13*c23


                              2
       - c13*c23*n2*r446 + c23 *n1*r446
r343=-----------------------------------
                    2          2
             b13*c13  + b13*c23


r344=0


r345=0


r346=0


r347=0


r348=0


r349=0


                           2
      c13*c23*n2*r446 - c23 *n1*r446
r350=--------------------------------
                  2          2
           b13*c13  + b13*c23


r351=0


r352=0


r353=0


r354=0


r355=0


           2    2           2    2
      - b13 *c13 *r446 - b13 *c23 *r446
r40=------------------------------------
                     2    2
                8*a33 *c23


     c23*r446
r41=----------
       a33


              2    4           4    2
      - 16*a33 *c23 *r446 - b13 *c13 *r446
r42=---------------------------------------
                    2    2    2
               8*a33 *b13 *c23


r43=0


r44=0


     c13*r446
r45=----------
       a33


        2
     b13 *c13*r446
r46=---------------
           2
      4*a33 *c23


r47=0


r48=0


              2    2           4
      - 16*a33 *c23 *r446 - b13 *r446
r49=----------------------------------
                    2    2
               8*a33 *b13


r410=0


r411=0


r412=0


r413=0


r415=0


r416=0


r417=0


r418=0


r419=0


r420=0


r421=0


r422=0


r423=0


r424=0


      r446
r425=------
       2


r426=0


r427=0


r428=0


r429=0


r431=0


r432=0


r433=0


       - b13*c13*r446
r435=-----------------
         2*a33*c23


r436=0


       - b13*c13*r446
r437=-----------------
         2*a33*c23


r438=0


      2*c23*r446
r439=------------
         b13


r441=0


r442=0


r444=0


r445=0


r447=0


r448=0


r449=0


r450=0


r451=0


r453=0


r454=0


       - r446
r455=---------
         2


r456=0


r458=0


r459=0


r460=0


r461=0


r462=0


r463=0


r464=0


r465=0


r466=0


r467=0


r468=0


r469=0


      b13*r446
r470=----------
       2*a33


       - 2*c23*r446
r471=---------------
           b13


r472=0


r473=0


r474=0


       - b13*c13*r446
r475=-----------------
         2*a33*c23


r476=0


      b13*r446
r477=----------
       2*a33


r478=0


r479=0


r480=0


r481=0


r482=0


r483=r446


r484=0


r485=0


r486=0


r487=0


r488=0


r489=0


r490=0


r491=0


r492=0


r493=0


r494=0


r495=0


r496=0


r497=0


r498=0


r499=0


r4100=0


r4101=0


r4102=0


r4103=0


r4104=0


        - r446
r4105=---------
          2


r4106=0


r4107=0


r4108=0


r4109=0


r4110=0


r4111=0


r4112=0


r4113=0


r4114=0


r4115=0


r4116=0


r4117=0


r4118=0


r4119=0


r4120=0


r4121=0


r4122=0


r4123=0


r4124=0


r4125=0


    b13*c13*n2 - b13*c23*n1
m3=-------------------------
           2*a33*c23


                                        2               3         2    2
m1=(2*a33*b13*c13*c23*m2 - 4*a33*c13*c23 *n1 - 4*a33*c23 *n2 + b13 *c13 *n3

          2    2                   2
     + b13 *c23 *n3)/(2*a33*b13*c23 )


             2    2      2    2
      - 2*b13 *c13  - b13 *c23
c33=----------------------------
                      2
             4*a33*c23


           2    2      2    2
      - b13 *c13  + b13 *c23
c22=--------------------------
                     2
            4*a33*c23


        2
     b13 *c13
c12=-----------
     4*a33*c23


c11=0


b33=0


b32=0


b31=0


      - b13*c13
b23=------------
        c23


b22=0


b21=0


b12=0


b11=0


a23=0


a22=0


a13=0


a12=0


a11=0


      b13*r446
r440=----------
       2*a33


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:
 r446, m2, n3, n2, n1, b13, c23, c13, a33

Inequalities

In the following not identically vanishing expressions are shown. Any auxiliary variables g00?? are used to express that at least one of their coefficients must not vanish, e.g. g0019*p4 + g0020*p3 means that either p4 or p3 or both are non-vanishing.
 
{n1,b13,a33,r440,n3}


Relevance for the application:

The new Hamiltonian in form of a list of vanishing expressions: 

{a11,
a12,
a13,
a22,
a23,
b11,
b12,
b21,
b22,
b13*c13 + b23*c23,
b31,
b32,
b33,
c11,
4*a33*c12*c23 - b13**2*c13,
4*a33*c22*c23**2 + b13**2*c13**2 - b13**2*c23**2,
4*a33*c23**2*c33 + 2*b13**2*c13**2 + b13**2*c23**2,
 - 2*a33*b13*c13*c23*m2 + 2*a33*b13*c23**2*m1 + 4*a33*c13*c23**2*n1 + 4*a33*c23
**3*n2 - b13**2*c13**2*n3 - b13**2*c23**2*n3,
2*a33*c23*m3 - b13*c13*n2 + b13*c23*n1}$

The system of equations related to the Hamiltonian HAM:

                          - u2*v3*b13*c13              2
HAM=u1*v3*b13 + u1*n1 + ------------------ + u2*n2 + u3 *a33 + u3*n3
                               c23

                 2
        v1*v2*b13 *c13
     + ---------------- + 2*v1*v3*c13 + v1*(2*a33*b13*c13*c23*m2
          2*a33*c23

                       2               3         2    2         2    2
        - 4*a33*c13*c23 *n1 - 4*a33*c23 *n2 + b13 *c13 *n3 + b13 *c23 *n3)/(2

                                   2    2      2    2
                   2      2   - b13 *c13  + b13 *c23
       *a33*b13*c23 ) + v2 *-------------------------- + 2*v2*v3*c23 + v2*m2
                                             2
                                    4*a33*c23

                    2    2      2    2
         2   - 2*b13 *c13  - b13 *c23         b13*c13*n2 - b13*c23*n1
     + v3 *---------------------------- + v3*-------------------------
                             2                       2*a33*c23
                    4*a33*c23

has apart from the Hamiltonian and Casimirs only the following first integral: 

     2   2          2    4    2    2        2    4    4
FI=u1 *v3 *( - 4*a33 *b13 *c13 *c23  - 4*a33 *b13 *c23 )

        2          2    3        3           2    3    4
    + u1 *v3*(8*a33 *b13 *c13*c23 *n2 - 8*a33 *b13 *c23 *n1)

        2          2    2    4   2        2    2    4   2
    + u1 *( - 4*a33 *b13 *c23 *n1  - 4*a33 *b13 *c23 *n2 )

                        2    4    2    2        2    4    4
    + u1*u3*v1*v3*(8*a33 *b13 *c13 *c23  + 8*a33 *b13 *c23 )

                        2    3        3           2    3    4
    + u1*u3*v1*( - 8*a33 *b13 *c13*c23 *n2 + 8*a33 *b13 *c23 *n1)

                     2    3        3           2    3    4
    + u1*u3*v2*(8*a33 *b13 *c13*c23 *n1 + 8*a33 *b13 *c23 *n2)

           2              5    2    2            5    4         2
    + u1*v1 *v3*(4*a33*b13 *c13 *c23  + 4*a33*b13 *c23 ) + u1*v1

              4    2    2               4        3               4    4
   *(8*a33*b13 *c13 *c23 *n1 + 4*a33*b13 *c13*c23 *n2 + 4*a33*b13 *c23 *n1)

                               5    3                5        3
    + u1*v1*v2*v3*( - 4*a33*b13 *c13 *c23 - 4*a33*b13 *c13*c23 )

                         4    2    2               4    4
    + u1*v1*v2*(4*a33*b13 *c13 *c23 *n2 + 4*a33*b13 *c23 *n2)

           2           4    2    2               4        3
    + u1*v2 *(4*a33*b13 *c13 *c23 *n1 + 4*a33*b13 *c13*c23 *n2)

              2           2    3    2    3         2    3    5
    + u1*v2*v3 *( - 16*a33 *b13 *c13 *c23  - 16*a33 *b13 *c23 ) + u1*v2*v3*(

              2    3    2    2           2    3    4
       - 8*a33 *b13 *c13 *c23 *m2 - 8*a33 *b13 *c23 *m2

               2    2    2    3            2    2        4
       + 16*a33 *b13 *c13 *c23 *n1 + 32*a33 *b13 *c13*c23 *n2

               2    2    5               4    3                   4        3
       - 16*a33 *b13 *c23 *n1 - 4*a33*b13 *c13 *c23*n3 - 4*a33*b13 *c13*c23 *n3)

                  2    2        3              2    2    4
    + u1*v2*(8*a33 *b13 *c13*c23 *m2*n2 - 8*a33 *b13 *c23 *m2*n1

                      2            4               2        5   2
              - 16*a33 *b13*c13*c23 *n1*n2 - 16*a33 *b13*c23 *n2

                         3    2    2                  3        3
              + 4*a33*b13 *c13 *c23 *n2*n3 + 4*a33*b13 *c13*c23 *n1*n3

                         3    4
              + 8*a33*b13 *c23 *n2*n3)

           3           5    2    2            5    4
    + u1*v3 *(4*a33*b13 *c13 *c23  + 4*a33*b13 *c23 )

           2           4    2    2               4    4
    + u1*v3 *(8*a33*b13 *c13 *c23 *n1 + 8*a33*b13 *c23 *n1)

                         3        3                  3    4   2
    + u1*v3*( - 8*a33*b13 *c13*c23 *n1*n2 + 8*a33*b13 *c23 *n1 )

                   2    4   3            2    4      2
    + u1*(4*a33*b13 *c23 *n1  + 4*a33*b13 *c23 *n1*n2 )

        2   2          2    4    2    2        2    4    4
    + u2 *v3 *( - 4*a33 *b13 *c13 *c23  - 4*a33 *b13 *c23 )

        2          2    3        3           2    3    4
    + u2 *v3*(8*a33 *b13 *c13*c23 *n2 - 8*a33 *b13 *c23 *n1)

        2          2    2    4   2        2    2    4   2
    + u2 *( - 4*a33 *b13 *c23 *n1  - 4*a33 *b13 *c23 *n2 )

                        2    3        3           2    3    4
    + u2*u3*v1*( - 8*a33 *b13 *c13*c23 *n1 - 8*a33 *b13 *c23 *n2)

                        2    4    2    2        2    4    4
    + u2*u3*v2*v3*(8*a33 *b13 *c13 *c23  + 8*a33 *b13 *c23 )

                        2    3        3           2    3    4
    + u2*u3*v2*( - 8*a33 *b13 *c13*c23 *n2 + 8*a33 *b13 *c23 *n1)

           2           4        3               4    4
    + u2*v1 *(4*a33*b13 *c13*c23 *n1 + 4*a33*b13 *c23 *n2)

                            5    2    2            5    4
    + u2*v1*v2*v3*(4*a33*b13 *c13 *c23  + 4*a33*b13 *c23 )

                         4    2    2               4    4
    + u2*v1*v2*(4*a33*b13 *c13 *c23 *n1 + 4*a33*b13 *c23 *n1)

              2        2    3    2    3         2    3    5
    + u2*v1*v3 *(16*a33 *b13 *c13 *c23  + 16*a33 *b13 *c23 ) + u2*v1*v3*(

           2    3    2    2           2    3    4            2    2    2    3
      8*a33 *b13 *c13 *c23 *m2 + 8*a33 *b13 *c23 *m2 - 16*a33 *b13 *c13 *c23 *n1

               2    2        4            2    2    5
       - 32*a33 *b13 *c13*c23 *n2 + 16*a33 *b13 *c23 *n1

                  4    3                   4        3
       + 4*a33*b13 *c13 *c23*n3 + 4*a33*b13 *c13*c23 *n3) + u2*v1*(

              2    2        3              2    2    4
       - 8*a33 *b13 *c13*c23 *m2*n2 + 8*a33 *b13 *c23 *m2*n1

               2            4               2        5   2
       + 16*a33 *b13*c13*c23 *n1*n2 + 16*a33 *b13*c23 *n2

                  3    2    2                  3        3
       - 4*a33*b13 *c13 *c23 *n2*n3 - 4*a33*b13 *c13*c23 *n1*n3

                  3    4
       - 8*a33*b13 *c23 *n2*n3)

           2                 5    3                5        3         2
    + u2*v2 *v3*( - 4*a33*b13 *c13 *c23 - 4*a33*b13 *c13*c23 ) + u2*v2

              4    2    2               4        3               4    4
   *(4*a33*b13 *c13 *c23 *n2 + 4*a33*b13 *c13*c23 *n1 + 8*a33*b13 *c23 *n2)

           3              5    3                5        3
    + u2*v3 *( - 4*a33*b13 *c13 *c23 - 4*a33*b13 *c13*c23 )

           2           4    2    2               4    4
    + u2*v3 *(8*a33*b13 *c13 *c23 *n2 + 8*a33*b13 *c23 *n2)

                         3        3   2            3    4
    + u2*v3*( - 8*a33*b13 *c13*c23 *n2  + 8*a33*b13 *c23 *n1*n2)

                   2    4   2               2    4   3
    + u2*(4*a33*b13 *c23 *n1 *n2 + 4*a33*b13 *c23 *n2 )

        2   2       2    4    2    2        2    4    4
    + u3 *v3 *(4*a33 *b13 *c13 *c23  + 4*a33 *b13 *c23 )

           2        2    2        4            2    2    5
    + u3*v1 *(16*a33 *b13 *c13*c23 *n1 + 16*a33 *b13 *c23 *n2)

                            4    2    2               4    4
    + u3*v1*v3*( - 4*a33*b13 *c13 *c23 *n1 - 4*a33*b13 *c23 *n1) + u3*v1

              3        3                  3    4   2            3    4   2
   *(8*a33*b13 *c13*c23 *n1*n2 - 4*a33*b13 *c23 *n1  + 4*a33*b13 *c23 *n2 )

           2        2    2        4            2    2    5
    + u3*v2 *(16*a33 *b13 *c13*c23 *n1 + 16*a33 *b13 *c23 *n2)

                            4    2    2               4    4
    + u3*v2*v3*( - 4*a33*b13 *c13 *c23 *n2 - 4*a33*b13 *c23 *n2) + u3*v2*(

                  3        3   2            3        3   2
       - 4*a33*b13 *c13*c23 *n1  + 4*a33*b13 *c13*c23 *n2

                  3    4                         2    2        3
       - 8*a33*b13 *c23 *n1*n2) + u3*v3*( - 8*a33 *b13 *c13*c23 *m2*n1

              2    2    4               2            4   2
       - 8*a33 *b13 *c23 *m2*n2 + 16*a33 *b13*c13*c23 *n1

               2        5                  3    2    2
       + 16*a33 *b13*c23 *n1*n2 - 4*a33*b13 *c13 *c23 *n1*n3

                  3        3
       - 4*a33*b13 *c13*c23 *n2*n3)

                   2    4   2               2    4   2
    + u3*(4*a33*b13 *c23 *n1 *n3 + 4*a33*b13 *c23 *n2 *n3)

        3              3        4               3    5
    + v1 *( - 8*a33*b13 *c13*c23 *n1 - 8*a33*b13 *c23 *n2)

        2              3    2    3               3        4         2   2
    + v1 *v2*(8*a33*b13 *c13 *c23 *n1 + 8*a33*b13 *c13*c23 *n2) + v1 *v3

              2    2    2    4         2    2    6      6    2    2      6    4
   *( - 16*a33 *b13 *c13 *c23  - 16*a33 *b13 *c23  - b13 *c13 *c23  - b13 *c23 )

        2              2    2    2    3            2    2    5
    + v1 *v3*( - 16*a33 *b13 *c13 *c23 *m2 - 16*a33 *b13 *c23 *m2

                       2        2    4            2            5
               + 32*a33 *b13*c13 *c23 *n1 + 32*a33 *b13*c13*c23 *n2

                          3    3    2               3        4
               - 8*a33*b13 *c13 *c23 *n3 - 8*a33*b13 *c13*c23 *n3

                      5    2    2           5    4
               - 2*b13 *c13 *c23 *n1 - 2*b13 *c23 *n1)

           2              3        4               3    5
    + v1*v2 *( - 8*a33*b13 *c13*c23 *n1 - 8*a33*b13 *c23 *n2)

              2       6    3            6        3
    + v1*v2*v3 *(2*b13 *c13 *c23 + 2*b13 *c13*c23 ) + v1*v2*v3*(

           5    3               5    2    2           5        3
      2*b13 *c13 *c23*n1 - 2*b13 *c13 *c23 *n2 + 2*b13 *c13*c23 *n1

              5    4                      4    2    2              4    4
       - 2*b13 *c23 *n2) + v1*v2*( - 4*b13 *c13 *c23 *n1*n2 - 4*b13 *c23 *n1*n2)

           3           4    3    2            4        4         2
    + v1*v3 *(8*a33*b13 *c13 *c23  + 8*a33*b13 *c13*c23 ) + v1*v3 *(

               4    3                   4        3               3    3    2
      4*a33*b13 *c13 *c23*m2 + 4*a33*b13 *c13*c23 *m2 - 8*a33*b13 *c13 *c23 *n1

                   3    2    3               3        4                3    5
       - 16*a33*b13 *c13 *c23 *n2 - 8*a33*b13 *c13*c23 *n1 - 16*a33*b13 *c23 *n2

              5    4           5    2    2
       + 2*b13 *c13 *n3 + 2*b13 *c13 *c23 *n3) + v1*v3*(

                  3    2    2                  3    4
       - 4*a33*b13 *c13 *c23 *m2*n2 - 4*a33*b13 *c23 *m2*n2

                  2    2    3                   2        4   2
       + 8*a33*b13 *c13 *c23 *n1*n2 + 16*a33*b13 *c13*c23 *n2

                  2    5              4    3                  4    2    2
       - 8*a33*b13 *c23 *n1*n2 - 2*b13 *c13 *c23*n2*n3 - 4*b13 *c13 *c23 *n1*n3

              4        3                          2        3      2
       - 6*b13 *c13*c23 *n2*n3) + v1*( - 4*a33*b13 *c13*c23 *m2*n1

                  2        3      2            2    4
       + 4*a33*b13 *c13*c23 *m2*n2  - 8*a33*b13 *c23 *m2*n1*n2

                          4   3                    4      2
       + 8*a33*b13*c13*c23 *n1  - 8*a33*b13*c13*c23 *n1*n2

                      5   2                   5   3        3    2    2   2
       + 8*a33*b13*c23 *n1 *n2 - 8*a33*b13*c23 *n2  - 2*b13 *c13 *c23 *n1 *n3

              3    2    2   2           3        3                 3    4   2
       + 2*b13 *c13 *c23 *n2 *n3 - 4*b13 *c13*c23 *n1*n2*n3 + 2*b13 *c23 *n1 *n3

              3    4   2
       + 2*b13 *c23 *n2 *n3)

        3           3    2    3               3        4         2   2
    + v2 *(8*a33*b13 *c13 *c23 *n1 + 8*a33*b13 *c13*c23 *n2) + v2 *v3

              2    2    2    4         2    2    6      6    4      6    2    2
   *( - 16*a33 *b13 *c13 *c23  - 16*a33 *b13 *c23  - b13 *c13  - b13 *c13 *c23 )

        2              2    2    2    3            2    2    5
    + v2 *v3*( - 16*a33 *b13 *c13 *c23 *m2 - 16*a33 *b13 *c23 *m2

                       2        2    4            2            5
               + 32*a33 *b13*c13 *c23 *n1 + 32*a33 *b13*c13*c23 *n2

                          3    3    2               3        4
               - 8*a33*b13 *c13 *c23 *n3 - 8*a33*b13 *c13*c23 *n3

                      5    3               5        3         2
               + 2*b13 *c13 *c23*n2 + 2*b13 *c13*c23 *n2) + v2 *(

           4    2    2   2        4    2    2   2        4    4   2
      2*b13 *c13 *c23 *n1  - 2*b13 *c13 *c23 *n2  + 2*b13 *c23 *n1

              4    4   2         3           4    2    3            4    5
       - 2*b13 *c23 *n2 ) + v2*v3 *(8*a33*b13 *c13 *c23  + 8*a33*b13 *c23 ) + v2

      2           4    2    2               4    4               3    2    3
   *v3 *(4*a33*b13 *c13 *c23 *m2 + 4*a33*b13 *c23 *m2 + 8*a33*b13 *c13 *c23 *n1

                     3    5           5    3               5        3
          + 8*a33*b13 *c23 *n1 + 2*b13 *c13 *c23*n3 + 2*b13 *c13*c23 *n3) + v2

                 3    2    2                  3    4
   *v3*(4*a33*b13 *c13 *c23 *m2*n1 + 4*a33*b13 *c23 *m2*n1

                    2    2    3   2             2        4
         - 8*a33*b13 *c13 *c23 *n1  - 16*a33*b13 *c13*c23 *n1*n2

                    2    5   2        4    3                  4        3
         + 8*a33*b13 *c23 *n1  + 2*b13 *c13 *c23*n1*n3 - 2*b13 *c13*c23 *n1*n3

                4    4                          2        3
         - 4*b13 *c23 *n2*n3) + v2*( - 8*a33*b13 *c13*c23 *m2*n1*n2

                  2    4      2            2    4      2
       + 4*a33*b13 *c23 *m2*n1  - 4*a33*b13 *c23 *m2*n2

                           4   2                    5      2
       + 16*a33*b13*c13*c23 *n1 *n2 + 16*a33*b13*c23 *n1*n2

              3    2    2                 3        3   2
       - 4*b13 *c13 *c23 *n1*n2*n3 - 4*b13 *c13*c23 *n2 *n3)

        4        6    4        6    2    2      6    4      3
    + v3 *( - b13 *c13  - 2*b13 *c13 *c23  - b13 *c23 ) + v3 *(

           5    3               5    2    2           5        3
      2*b13 *c13 *c23*n2 - 2*b13 *c13 *c23 *n1 + 2*b13 *c13*c23 *n2

              5    4         2       2    2    2    2   2        2    2    4   2
       - 2*b13 *c23 *n1) + v3 *(4*a33 *b13 *c13 *c23 *m2  + 4*a33 *b13 *c23 *m2

               2        2    3               2            4
       - 16*a33 *b13*c13 *c23 *m2*n1 - 16*a33 *b13*c13*c23 *m2*n2

               2    2    4   2         2        5               2    6   2
       + 16*a33 *c13 *c23 *n1  + 32*a33 *c13*c23 *n1*n2 + 16*a33 *c23 *n2

                  3    3                      3        3
       + 4*a33*b13 *c13 *c23*m2*n3 + 4*a33*b13 *c13*c23 *m2*n3

                  2    3    2                  2    2    3
       - 8*a33*b13 *c13 *c23 *n1*n3 - 8*a33*b13 *c13 *c23 *n2*n3

                   2        4                   2    5            4    4   2
       - 16*a33*b13 *c13*c23 *n1*n3 - 16*a33*b13 *c23 *n2*n3 + b13 *c13 *n3

            4    2    2   2        4    2    2   2      4    2    2   2
       + b13 *c13 *c23 *n1  - 3*b13 *c13 *c23 *n2  + b13 *c13 *c23 *n3

              4        3            4    4   2      4    4   2
       + 4*b13 *c13*c23 *n1*n2 - b13 *c23 *n1  - b13 *c23 *n2 ) + v3*(

                  2        3                     2    4
       - 8*a33*b13 *c13*c23 *m2*n1*n3 - 8*a33*b13 *c23 *m2*n2*n3

                           4   2                    5
       + 16*a33*b13*c13*c23 *n1 *n3 + 16*a33*b13*c23 *n1*n2*n3

              3    2    2      2        3        3   2           3        3   3
       - 4*b13 *c13 *c23 *n1*n3  + 2*b13 *c13*c23 *n1 *n2 + 2*b13 *c13*c23 *n2

              3        3      2        3    4   3        3    4      2
       - 4*b13 *c13*c23 *n2*n3  - 2*b13 *c23 *n1  - 2*b13 *c23 *n1*n2 )
Stop of a subroutine.
Number of garbage collections exceeds max_gc_fac.
The factorization crashed!

which the program can not factorize further.

{HAM,FI} = 0





And again in machine readable form:



HAM=u1*v3*b13 + u1*n1 + ( - u2*v3*b13*c13)/c23 + u2*n2 + u3**2*a33 + u3*n3 + (v1
*v2*b13**2*c13)/(2*a33*c23) + 2*v1*v3*c13 + v1*(2*a33*b13*c13*c23*m2 - 4*a33*c13
*c23**2*n1 - 4*a33*c23**3*n2 + b13**2*c13**2*n3 + b13**2*c23**2*n3)/(2*a33*b13*
c23**2) + v2**2*( - b13**2*c13**2 + b13**2*c23**2)/(4*a33*c23**2) + 2*v2*v3*c23 
+ v2*m2 + v3**2*( - 2*b13**2*c13**2 - b13**2*c23**2)/(4*a33*c23**2) + v3*(b13*
c13*n2 - b13*c23*n1)/(2*a33*c23)$

FI=u1**2*v3**2*( - 4*a33**2*b13**4*c13**2*c23**2 - 4*a33**2*b13**4*c23**4) + u1
**2*v3*(8*a33**2*b13**3*c13*c23**3*n2 - 8*a33**2*b13**3*c23**4*n1) + u1**2*( - 4
*a33**2*b13**2*c23**4*n1**2 - 4*a33**2*b13**2*c23**4*n2**2) + u1*u3*v1*v3*(8*a33
**2*b13**4*c13**2*c23**2 + 8*a33**2*b13**4*c23**4) + u1*u3*v1*( - 8*a33**2*b13**
3*c13*c23**3*n2 + 8*a33**2*b13**3*c23**4*n1) + u1*u3*v2*(8*a33**2*b13**3*c13*c23
**3*n1 + 8*a33**2*b13**3*c23**4*n2) + u1*v1**2*v3*(4*a33*b13**5*c13**2*c23**2 + 
4*a33*b13**5*c23**4) + u1*v1**2*(8*a33*b13**4*c13**2*c23**2*n1 + 4*a33*b13**4*
c13*c23**3*n2 + 4*a33*b13**4*c23**4*n1) + u1*v1*v2*v3*( - 4*a33*b13**5*c13**3*
c23 - 4*a33*b13**5*c13*c23**3) + u1*v1*v2*(4*a33*b13**4*c13**2*c23**2*n2 + 4*a33
*b13**4*c23**4*n2) + u1*v2**2*(4*a33*b13**4*c13**2*c23**2*n1 + 4*a33*b13**4*c13*
c23**3*n2) + u1*v2*v3**2*( - 16*a33**2*b13**3*c13**2*c23**3 - 16*a33**2*b13**3*
c23**5) + u1*v2*v3*( - 8*a33**2*b13**3*c13**2*c23**2*m2 - 8*a33**2*b13**3*c23**4
*m2 + 16*a33**2*b13**2*c13**2*c23**3*n1 + 32*a33**2*b13**2*c13*c23**4*n2 - 16*
a33**2*b13**2*c23**5*n1 - 4*a33*b13**4*c13**3*c23*n3 - 4*a33*b13**4*c13*c23**3*
n3) + u1*v2*(8*a33**2*b13**2*c13*c23**3*m2*n2 - 8*a33**2*b13**2*c23**4*m2*n1 - 
16*a33**2*b13*c13*c23**4*n1*n2 - 16*a33**2*b13*c23**5*n2**2 + 4*a33*b13**3*c13**
2*c23**2*n2*n3 + 4*a33*b13**3*c13*c23**3*n1*n3 + 8*a33*b13**3*c23**4*n2*n3) + u1
*v3**3*(4*a33*b13**5*c13**2*c23**2 + 4*a33*b13**5*c23**4) + u1*v3**2*(8*a33*b13
**4*c13**2*c23**2*n1 + 8*a33*b13**4*c23**4*n1) + u1*v3*( - 8*a33*b13**3*c13*c23
**3*n1*n2 + 8*a33*b13**3*c23**4*n1**2) + u1*(4*a33*b13**2*c23**4*n1**3 + 4*a33*
b13**2*c23**4*n1*n2**2) + u2**2*v3**2*( - 4*a33**2*b13**4*c13**2*c23**2 - 4*a33
**2*b13**4*c23**4) + u2**2*v3*(8*a33**2*b13**3*c13*c23**3*n2 - 8*a33**2*b13**3*
c23**4*n1) + u2**2*( - 4*a33**2*b13**2*c23**4*n1**2 - 4*a33**2*b13**2*c23**4*n2
**2) + u2*u3*v1*( - 8*a33**2*b13**3*c13*c23**3*n1 - 8*a33**2*b13**3*c23**4*n2) +
 u2*u3*v2*v3*(8*a33**2*b13**4*c13**2*c23**2 + 8*a33**2*b13**4*c23**4) + u2*u3*v2
*( - 8*a33**2*b13**3*c13*c23**3*n2 + 8*a33**2*b13**3*c23**4*n1) + u2*v1**2*(4*
a33*b13**4*c13*c23**3*n1 + 4*a33*b13**4*c23**4*n2) + u2*v1*v2*v3*(4*a33*b13**5*
c13**2*c23**2 + 4*a33*b13**5*c23**4) + u2*v1*v2*(4*a33*b13**4*c13**2*c23**2*n1 +
 4*a33*b13**4*c23**4*n1) + u2*v1*v3**2*(16*a33**2*b13**3*c13**2*c23**3 + 16*a33
**2*b13**3*c23**5) + u2*v1*v3*(8*a33**2*b13**3*c13**2*c23**2*m2 + 8*a33**2*b13**
3*c23**4*m2 - 16*a33**2*b13**2*c13**2*c23**3*n1 - 32*a33**2*b13**2*c13*c23**4*n2
 + 16*a33**2*b13**2*c23**5*n1 + 4*a33*b13**4*c13**3*c23*n3 + 4*a33*b13**4*c13*
c23**3*n3) + u2*v1*( - 8*a33**2*b13**2*c13*c23**3*m2*n2 + 8*a33**2*b13**2*c23**4
*m2*n1 + 16*a33**2*b13*c13*c23**4*n1*n2 + 16*a33**2*b13*c23**5*n2**2 - 4*a33*b13
**3*c13**2*c23**2*n2*n3 - 4*a33*b13**3*c13*c23**3*n1*n3 - 8*a33*b13**3*c23**4*n2
*n3) + u2*v2**2*v3*( - 4*a33*b13**5*c13**3*c23 - 4*a33*b13**5*c13*c23**3) + u2*
v2**2*(4*a33*b13**4*c13**2*c23**2*n2 + 4*a33*b13**4*c13*c23**3*n1 + 8*a33*b13**4
*c23**4*n2) + u2*v3**3*( - 4*a33*b13**5*c13**3*c23 - 4*a33*b13**5*c13*c23**3) + 
u2*v3**2*(8*a33*b13**4*c13**2*c23**2*n2 + 8*a33*b13**4*c23**4*n2) + u2*v3*( - 8*
a33*b13**3*c13*c23**3*n2**2 + 8*a33*b13**3*c23**4*n1*n2) + u2*(4*a33*b13**2*c23
**4*n1**2*n2 + 4*a33*b13**2*c23**4*n2**3) + u3**2*v3**2*(4*a33**2*b13**4*c13**2*
c23**2 + 4*a33**2*b13**4*c23**4) + u3*v1**2*(16*a33**2*b13**2*c13*c23**4*n1 + 16
*a33**2*b13**2*c23**5*n2) + u3*v1*v3*( - 4*a33*b13**4*c13**2*c23**2*n1 - 4*a33*
b13**4*c23**4*n1) + u3*v1*(8*a33*b13**3*c13*c23**3*n1*n2 - 4*a33*b13**3*c23**4*
n1**2 + 4*a33*b13**3*c23**4*n2**2) + u3*v2**2*(16*a33**2*b13**2*c13*c23**4*n1 + 
16*a33**2*b13**2*c23**5*n2) + u3*v2*v3*( - 4*a33*b13**4*c13**2*c23**2*n2 - 4*a33
*b13**4*c23**4*n2) + u3*v2*( - 4*a33*b13**3*c13*c23**3*n1**2 + 4*a33*b13**3*c13*
c23**3*n2**2 - 8*a33*b13**3*c23**4*n1*n2) + u3*v3*( - 8*a33**2*b13**2*c13*c23**3
*m2*n1 - 8*a33**2*b13**2*c23**4*m2*n2 + 16*a33**2*b13*c13*c23**4*n1**2 + 16*a33
**2*b13*c23**5*n1*n2 - 4*a33*b13**3*c13**2*c23**2*n1*n3 - 4*a33*b13**3*c13*c23**
3*n2*n3) + u3*(4*a33*b13**2*c23**4*n1**2*n3 + 4*a33*b13**2*c23**4*n2**2*n3) + v1
**3*( - 8*a33*b13**3*c13*c23**4*n1 - 8*a33*b13**3*c23**5*n2) + v1**2*v2*(8*a33*
b13**3*c13**2*c23**3*n1 + 8*a33*b13**3*c13*c23**4*n2) + v1**2*v3**2*( - 16*a33**
2*b13**2*c13**2*c23**4 - 16*a33**2*b13**2*c23**6 - b13**6*c13**2*c23**2 - b13**6
*c23**4) + v1**2*v3*( - 16*a33**2*b13**2*c13**2*c23**3*m2 - 16*a33**2*b13**2*c23
**5*m2 + 32*a33**2*b13*c13**2*c23**4*n1 + 32*a33**2*b13*c13*c23**5*n2 - 8*a33*
b13**3*c13**3*c23**2*n3 - 8*a33*b13**3*c13*c23**4*n3 - 2*b13**5*c13**2*c23**2*n1
 - 2*b13**5*c23**4*n1) + v1*v2**2*( - 8*a33*b13**3*c13*c23**4*n1 - 8*a33*b13**3*
c23**5*n2) + v1*v2*v3**2*(2*b13**6*c13**3*c23 + 2*b13**6*c13*c23**3) + v1*v2*v3*
(2*b13**5*c13**3*c23*n1 - 2*b13**5*c13**2*c23**2*n2 + 2*b13**5*c13*c23**3*n1 - 2
*b13**5*c23**4*n2) + v1*v2*( - 4*b13**4*c13**2*c23**2*n1*n2 - 4*b13**4*c23**4*n1
*n2) + v1*v3**3*(8*a33*b13**4*c13**3*c23**2 + 8*a33*b13**4*c13*c23**4) + v1*v3**
2*(4*a33*b13**4*c13**3*c23*m2 + 4*a33*b13**4*c13*c23**3*m2 - 8*a33*b13**3*c13**3
*c23**2*n1 - 16*a33*b13**3*c13**2*c23**3*n2 - 8*a33*b13**3*c13*c23**4*n1 - 16*
a33*b13**3*c23**5*n2 + 2*b13**5*c13**4*n3 + 2*b13**5*c13**2*c23**2*n3) + v1*v3*(
 - 4*a33*b13**3*c13**2*c23**2*m2*n2 - 4*a33*b13**3*c23**4*m2*n2 + 8*a33*b13**2*
c13**2*c23**3*n1*n2 + 16*a33*b13**2*c13*c23**4*n2**2 - 8*a33*b13**2*c23**5*n1*n2
 - 2*b13**4*c13**3*c23*n2*n3 - 4*b13**4*c13**2*c23**2*n1*n3 - 6*b13**4*c13*c23**
3*n2*n3) + v1*( - 4*a33*b13**2*c13*c23**3*m2*n1**2 + 4*a33*b13**2*c13*c23**3*m2*
n2**2 - 8*a33*b13**2*c23**4*m2*n1*n2 + 8*a33*b13*c13*c23**4*n1**3 - 8*a33*b13*
c13*c23**4*n1*n2**2 + 8*a33*b13*c23**5*n1**2*n2 - 8*a33*b13*c23**5*n2**3 - 2*b13
**3*c13**2*c23**2*n1**2*n3 + 2*b13**3*c13**2*c23**2*n2**2*n3 - 4*b13**3*c13*c23
**3*n1*n2*n3 + 2*b13**3*c23**4*n1**2*n3 + 2*b13**3*c23**4*n2**2*n3) + v2**3*(8*
a33*b13**3*c13**2*c23**3*n1 + 8*a33*b13**3*c13*c23**4*n2) + v2**2*v3**2*( - 16*
a33**2*b13**2*c13**2*c23**4 - 16*a33**2*b13**2*c23**6 - b13**6*c13**4 - b13**6*
c13**2*c23**2) + v2**2*v3*( - 16*a33**2*b13**2*c13**2*c23**3*m2 - 16*a33**2*b13
**2*c23**5*m2 + 32*a33**2*b13*c13**2*c23**4*n1 + 32*a33**2*b13*c13*c23**5*n2 - 8
*a33*b13**3*c13**3*c23**2*n3 - 8*a33*b13**3*c13*c23**4*n3 + 2*b13**5*c13**3*c23*
n2 + 2*b13**5*c13*c23**3*n2) + v2**2*(2*b13**4*c13**2*c23**2*n1**2 - 2*b13**4*
c13**2*c23**2*n2**2 + 2*b13**4*c23**4*n1**2 - 2*b13**4*c23**4*n2**2) + v2*v3**3*
(8*a33*b13**4*c13**2*c23**3 + 8*a33*b13**4*c23**5) + v2*v3**2*(4*a33*b13**4*c13
**2*c23**2*m2 + 4*a33*b13**4*c23**4*m2 + 8*a33*b13**3*c13**2*c23**3*n1 + 8*a33*
b13**3*c23**5*n1 + 2*b13**5*c13**3*c23*n3 + 2*b13**5*c13*c23**3*n3) + v2*v3*(4*
a33*b13**3*c13**2*c23**2*m2*n1 + 4*a33*b13**3*c23**4*m2*n1 - 8*a33*b13**2*c13**2
*c23**3*n1**2 - 16*a33*b13**2*c13*c23**4*n1*n2 + 8*a33*b13**2*c23**5*n1**2 + 2*
b13**4*c13**3*c23*n1*n3 - 2*b13**4*c13*c23**3*n1*n3 - 4*b13**4*c23**4*n2*n3) + 
v2*( - 8*a33*b13**2*c13*c23**3*m2*n1*n2 + 4*a33*b13**2*c23**4*m2*n1**2 - 4*a33*
b13**2*c23**4*m2*n2**2 + 16*a33*b13*c13*c23**4*n1**2*n2 + 16*a33*b13*c23**5*n1*
n2**2 - 4*b13**3*c13**2*c23**2*n1*n2*n3 - 4*b13**3*c13*c23**3*n2**2*n3) + v3**4*
( - b13**6*c13**4 - 2*b13**6*c13**2*c23**2 - b13**6*c23**4) + v3**3*(2*b13**5*
c13**3*c23*n2 - 2*b13**5*c13**2*c23**2*n1 + 2*b13**5*c13*c23**3*n2 - 2*b13**5*
c23**4*n1) + v3**2*(4*a33**2*b13**2*c13**2*c23**2*m2**2 + 4*a33**2*b13**2*c23**4
*m2**2 - 16*a33**2*b13*c13**2*c23**3*m2*n1 - 16*a33**2*b13*c13*c23**4*m2*n2 + 16
*a33**2*c13**2*c23**4*n1**2 + 32*a33**2*c13*c23**5*n1*n2 + 16*a33**2*c23**6*n2**
2 + 4*a33*b13**3*c13**3*c23*m2*n3 + 4*a33*b13**3*c13*c23**3*m2*n3 - 8*a33*b13**2
*c13**3*c23**2*n1*n3 - 8*a33*b13**2*c13**2*c23**3*n2*n3 - 16*a33*b13**2*c13*c23
**4*n1*n3 - 16*a33*b13**2*c23**5*n2*n3 + b13**4*c13**4*n3**2 + b13**4*c13**2*c23
**2*n1**2 - 3*b13**4*c13**2*c23**2*n2**2 + b13**4*c13**2*c23**2*n3**2 + 4*b13**4
*c13*c23**3*n1*n2 - b13**4*c23**4*n1**2 - b13**4*c23**4*n2**2) + v3*( - 8*a33*
b13**2*c13*c23**3*m2*n1*n3 - 8*a33*b13**2*c23**4*m2*n2*n3 + 16*a33*b13*c13*c23**
4*n1**2*n3 + 16*a33*b13*c23**5*n1*n2*n3 - 4*b13**3*c13**2*c23**2*n1*n3**2 + 2*
b13**3*c13*c23**3*n1**2*n2 + 2*b13**3*c13*c23**3*n2**3 - 4*b13**3*c13*c23**3*n2*
n3**2 - 2*b13**3*c23**4*n1**3 - 2*b13**3*c23**4*n1*n2**2)$