File r38/packages/crack/liepde.rlg artifact 317cd24e71 part of check-in b5833487d7


Tue Feb 10 12:26:55 2004 run on Linux
%*******************************************************************%
%                                                                   %
%                      L I E P D E . T S T                          %
%                      -------------------                          %
%  liepde.tst contains test examples for the program liepde.red.    %
%                                                                   %
%  Author of this file: Thomas Wolf                                 %
%  Date:  21. April 1998, 6. May 2003                               %
%                                                                   %
%  Details about the syntax of liepde.red are given in liepde.tex.  %
%                                                                   %
%  To run this demo you need to load liepde and crack through       %
%     load crack,liepde$                                            %
%  and to read in this file as                                      %
%     in "liepde.tst";                                              %
%  If you got the source code of a newer version of liepde then     %
%  either read it in through                                        %
%     in "crack.red","liepde.red"$                                  %
%  (with the appropriate directory name in front of liepde.red)     %
%  or, to speed up the calculation, you compile before with         %
%     faslout "crack"$                                              %
%     in "crack.red"$                                               %
%     faslend$                                                      %
%     faslout "liepde"$                                             %
%     in "liepde.red"$                                              %
%     faslend$                                                      %
%  and then load both it with                                       %
%     load crack,liepde$                                            %
%                                                                   %
%*******************************************************************%

load crack;



lisp(depl!*:=nil)$

     % clearing of all dependences
setcrackflags()$


lisp(print_:=nil)$


on dfprint$



comment
-------------------------------------------------------
The following runs demonstrate the program LIEPDE for 
the computation of infinitesimal symmetries. Times given
below refer to a 8 MB session under LINUX on a 133 MHz
Pentium PC with the CRACK version of April 1998 running
PSL Reduce.
-------------------------------------------------------;


lisp(prelim_:=nil)$

      % not necessary as this is the default value
lisp(individual_:=nil)$

  % not necessary as this is the default value

comment
-------------------------------------------------------
The first example is a single ODE with a parametric
function f=f(x) for which point symmetries are to be
determined.
(Time ~ 6 sec.);

write"-------------------------------------------------------";


-------------------------------------------------------


lisp(freeint_:=nil)$

 % This enables the solution of differential equ.s in
                     % which unevaluated integrals remain. This becomes
                     % necessary in this example through the parametric
                     % function f=f(x)
depend y,x$


depend f,x$


liepde({df(y,x,2)=-(y+3*f)*df(y,x)+y**3-f*y**2-(2*f**2+df(f,x))*y, 
        {y}, {x}},
       {"point"},{},{})$



The ODE under investigation is :

                                 2        2    3
y  = - f *y - 3*y *f - y *y - 2*f *y - f*y  + y
 2x     x        x      x

for the function(s) : 

y(x)



The symmetries are:

-------- 1. Symmetry:

         int(f,x)
xi_x= - e

       int(f,x)
eta_y=e        *f*y

-------- 2. Symmetry:

         int(f,x)          1
xi_x= - e        *int(-----------,x)
                        int(f,x)
                       e

       int(f,x)          1
eta_y=e        *int(-----------,x)*f*y + y
                      int(f,x)
                     e

-------- 


{{},

           int(f,x)          1                int(f,x)
 {xi_x= - e        *int(-----------,x)*c_5 - e        *c_6,
                          int(f,x)
                         e

         int(f,x)          1                    int(f,x)
  eta_y=e        *int(-----------,x)*c_5*f*y + e        *c_6*f*y + c_5*y},
                        int(f,x)
                       e

 {c_6,c_5}}

nodepnd {y,f}$


lisp(freeint_:=t)$

   % Because the simplification of differential
                     % expressions which involve unevaluated integrals
                     % may provide difficulties such solutions involving
                     % unevaluated integrals are disabled.

comment
-------------------------------------------------------
The following example demonstrates a number of things.
The Burgers equation is investigated concerning third
order symmetries. The equation is used to substitute
df(u,t) and all derivatives of df(u,t). This computation
also shows that any equations that remain unsolved are
returned, like in this case the heat quation.
(Time ~ 15 sec.);

write"-------------------------------------------------------";


-------------------------------------------------------


nodepnd {u}$


depend u,t,x$


liepde({df(u,t)=df(u,x,2)+df(u,x)**2,{u},{t,x}},{"general",3},{},{})$



The PDE under investigation is :

           2
u =u   + u
 t  2x    x

for the function(s) : 

u(t,x)



The symmetries are:

-------- 1. Symmetry:

xi_t=0

xi_x=0

                    2
eta_u=2*u  *t + 2*u  *t + u *x
         2x        x       x

-------- 2. Symmetry:

xi_t=0

xi_x=0

             2       2  2                     2
eta_u=4*u  *t  + 4*u  *t  + 4*u *t*x + 2*t + x
         2x         x          x

-------- 3. Symmetry:

xi_t=0

xi_x=0

              2
eta_u=u   + u
       2x    x

-------- 4. Symmetry:

xi_t=0

xi_x=0

                                         3       2
eta_u=2*u  *t + 6*u  *u *t + u  *x + 2*u  *t + u  *x
         3x        2x  x      2x        x       x

-------- 5. Symmetry:

xi_t=0

xi_x=0

             3              3           2         3  3        2  2            2
eta_u=8*u  *t  + 24*u  *u *t  + 12*u  *t *x + 8*u  *t  + 12*u  *t *x + 12*u *t
         3x          2x  x          2x           x           x             x

           2            3
 + 6*u *t*x  + 6*t*x + x
      x

-------- 6. Symmetry:

xi_t=0

xi_x=0

             2              2                   3  2       2           2
eta_u=4*u  *t  + 12*u  *u *t  + 4*u  *t*x + 4*u  *t  + 4*u  *t*x + u *x  - x
         3x          2x  x         2x          x          x         x

-------- 7. Symmetry:

xi_t=0

xi_x=0

                         3
eta_u=u   + 3*u  *u  + u
       3x      2x  x    x

-------- 8. Symmetry:

xi_t=0

xi_x=0

eta_u=2*u *t + x
         x

-------- 9. Symmetry:

xi_t=0

xi_x=0

eta_u=u
       x

-------- 10. Symmetry:

xi_t=0

xi_x=0

eta_u=1

-------- 

Further symmetries:

xi_t=0

xi_x=0

       c_25
eta_u=------
         u
        e

with c_25(t,x)
which still have to satisfy: 

0=c_25  - c_25
      t       2x


{{c_25  - c_25  },
      t       2x

 {xi_t=0,

  xi_x=0,

          u               u           2      u           3      u
  eta_u=(e *u  *c_31 + 4*e *u  *c_32*t  + 8*e *u  *c_33*t  + 2*e *u  *c_34*t
             3x              3x                 3x                 3x

               u                   u              2       u              3
          + 3*e *u  *u *c_31 + 12*e *u  *u *c_32*t  + 24*e *u  *u *c_33*t
                  2x  x               2x  x                  2x  x

               u                  u               u           2
          + 6*e *u  *u *c_34*t + e *u  *c_28 + 4*e *u  *c_29*t
                  2x  x              2x              2x

               u                 u                    u           2
          + 2*e *u  *c_30*t + 4*e *u  *c_32*t*x + 12*e *u  *c_33*t *x
                  2x                2x                   2x

             u               u   3           u   3       2      u   3       3
          + e *u  *c_34*x + e *u  *c_31 + 4*e *u  *c_32*t  + 8*e *u  *c_33*t
                2x              x               x                  x

               u   3           u   2           u   2       2      u   2
          + 2*e *u  *c_34*t + e *u  *c_28 + 4*e *u  *c_29*t  + 2*e *u  *c_30*t
                  x               x               x                  x

               u   2                u   2       2      u   2
          + 4*e *u  *c_32*t*x + 12*e *u  *c_33*t *x + e *u  *c_34*x
                  x                    x                  x

               u                u              u          2       u          2
          + 4*e *u *c_29*t*x + e *u *c_30*x + e *u *c_32*x  + 12*e *u *c_33*t
                  x                x              x                  x

               u            2    u              u                u
          + 6*e *u *c_33*t*x  + e *u *c_35 + 2*e *u *c_36*t + 2*e *c_29*t
                  x                 x              x

             u       2    u             u             u       3    u
          + e *c_29*x  - e *c_32*x + 6*e *c_33*t*x + e *c_33*x  + e *c_36*x

             u               u
          + e *c_37 + c_25)/e },

 {c_25,

  c_30,

  c_29,

  c_28,

  c_34,

  c_33,

  c_32,

  c_31,

  c_36,

  c_35,

  c_37}}


comment
-------------------------------------------------------
Now the same equation is investigated, this time only
df(u,x,2) and its derivatives are substituted. As a
consequence less jet-variables (u-derivatives of lower
order) are generated in the process of formulating the
symmetry conditions. Less jet-variables in which the
conditions have to be fulfilled identically means less
overdetermined conditions and more solutions which to
compute takes longer than before.
(Time ~ 85 sec.);

write"-------------------------------------------------------";


-------------------------------------------------------


liepde({df(u,x,2)=df(u,t)-df(u,x)**2,{u},{t,x}},{"general",3},{},{})$



The PDE under investigation is :

           2
u  =u  - u
 2x  t    x

for the function(s) : 

u(t,x)



The symmetries are:

-------- 1. Symmetry:

xi_t=0

xi_x=0

                                                                     2
eta_u= - 8*u   *u *t + u    - 4*u    *t - 4*u  *u *t - 2*u  *x - 4*u  *u *t
            t,x  t      t,x      2t,x        2t  x        2t        t   x

       2
 - 2*u  *x + u *u
      t       t  x

-------- 2. Symmetry:

xi_t=0

xi_x=0

                     3            2              2             3              3
eta_u= - 32*u   *u *t  + 12*u   *t  - 12*u   *t*x  - 16*u    *t  - 16*u  *u *t
             t,x  t          t,x          t,x            2t,x          2t  x

           2          2     3        2  2               2               2
 - 24*u  *t *x - 16*u  *u *t  - 24*u  *t *x + 12*u *u *t  - 12*u *u *t*x
       2t            t   x          t             t  x          t  x

                     3         2
 + 24*u *t*x - 2*u *x  + 9*u *x  - 9*x
       t          t         x

-------- 3. Symmetry:

xi_t=0

xi_x=0

                     5             4            3  2             5
eta_u= - 64*u   *u *t  - 160*u   *t  - 80*u   *t *x  - 32*u    *t
             t,x  t           t,x          t,x             2t,x

              5           4          2     5        2  4                4
 - 32*u  *u *t  - 80*u  *t *x - 32*u  *u *t  - 80*u  *t *x - 160*u *u *t
       2t  x          2t            t   x          t              t  x

             3  2           3            2  3           3           2  2
 - 80*u *u *t *x  - 240*u *t *x - 40*u *t *x  - 120*u *t  - 120*u *t *x
       t  x              t            t              x           x

            4       2           3    5
 - 10*u *t*x  - 60*t *x - 20*t*x  - x
       x

-------- 4. Symmetry:

xi_t=0

xi_x=0

                     4            3            2  2             4              4
eta_u= - 32*u   *u *t  + 16*u   *t  - 24*u   *t *x  - 16*u    *t  - 16*u  *u *t
             t,x  t          t,x          t,x             2t,x          2t  x

           3          2     4        2  3               3             2  2
 - 32*u  *t *x - 16*u  *u *t  - 32*u  *t *x + 16*u *u *t  - 24*u *u *t *x
       2t            t   x          t             t  x          t  x

          2             3          2            2       4               3
 + 48*u *t *x - 8*u *t*x  + 84*u *t  + 36*u *t*x  - u *x  + 48*t*x + 8*x
       t           t            x          x         x

-------- 5. Symmetry:

xi_t=0

xi_x=0

                    2                    2            2             2
eta_u= - 8*u   *u *t  + 2*u   *t - u   *x  - 4*u    *t  - 4*u  *u *t
            t,x  t         t,x      t,x         2t,x         2t  x

                   2     2       2                          2
 - 4*u  *t*x - 4*u  *u *t  - 4*u  *t*x + 2*u *u *t - u *u *x  + 2*u *x
      2t          t   x         t           t  x      t  x         t

-------- 6. Symmetry:

xi_t=0

xi_x=0

                                        2
eta_u= - 2*u   *u  - u     - u  *u  - u  *u
            t,x  t    2t,x    2t  x    t   x

-------- 7. Symmetry:

xi_t=0

xi_x=0

eta_u

                                                                  3       2
= - 2*u   *u *x - 2*u  *t - u    *x - 6*u  *u *t - u  *u *x - 2*u  *t - u  *u *x
       t,x  t        3t      2t,x        2t  t      2t  x        t       t   x

-------- 8. Symmetry:

xi_t=0

xi_x=0

                     2                           3           3             2
eta_u= - 48*u   *u *t *x + 12*u   *t*x - 2*u   *x  - 16*u  *t  - 24*u    *t *x
             t,x  t            t,x          t,x          3t          2t,x

              3              2               2        3  3        2     2
 - 48*u  *u *t  - 24*u  *u *t *x - 12*u  *t*x  - 16*u  *t  - 24*u  *u *t *x
       2t  t          2t  x            2t            t           t   x

        2    2                           3         2
 - 12*u  *t*x  + 12*u *u *t*x - 2*u *u *x  + 6*u *x  - 6*u *x + 3
       t             t  x          t  x         t         x

-------- 9. Symmetry:

xi_t=0

xi_x=0

                      4              3              2  3           5
eta_u= - 160*u   *u *t *x + 80*u   *t *x - 40*u   *t *x  - 32*u  *t
              t,x  t            t,x            t,x             3t

             4                5              4             3  2        3  5
 - 80*u    *t *x - 96*u  *u *t  - 80*u  *u *t *x - 80*u  *t *x  - 32*u  *t
       2t,x            2t  t          2t  x            2t             t

        2     4          2  3  2             3               2  3           3
 - 80*u  *u *t *x - 80*u  *t *x  + 80*u *u *t *x - 40*u *u *t *x  + 360*u *t
       t   x            t              t  x            t  x              t

           2  2            4           2              3       5        2
 + 120*u *t *x  - 10*u *t*x  + 420*u *t *x + 60*u *t*x  - u *x  + 120*t
        t             t             x            x         x

          2       4
 + 120*t*x  + 10*x

-------- 10. Symmetry:

xi_t=0

xi_x=0

                      5               4               3  3           6
eta_u= - 384*u   *u *t *x - 960*u   *t *x - 160*u   *t *x  - 64*u  *t
              t,x  t             t,x             t,x             3t

              5                 6               5              5            4  2
 - 192*u    *t *x - 192*u  *u *t  - 192*u  *u *t *x - 480*u  *t  - 240*u  *t *x
        2t,x             2t  t           2t  x             2t           2t

        3  6         2     5           2  5         2  4  2              4
 - 64*u  *t  - 192*u  *u *t *x - 480*u  *t  - 240*u  *t *x  - 960*u *u *t *x
       t            t   x             t            t               t  x

              3  3           4           3  2          2  4           3
 - 160*u *u *t *x  - 720*u *t  - 720*u *t *x  - 60*u *t *x  - 720*u *t *x
        t  x              t           t             t              x

           2  3            5        3        2  2         4    6
 - 240*u *t *x  - 12*u *t*x  - 120*t  - 180*t *x  - 30*t*x  - x
        x             x

-------- 11. Symmetry:

xi_t=0

xi_x=0

                     3              2               3           4
eta_u= - 64*u   *u *t *x + 24*u   *t *x - 8*u   *t*x  - 16*u  *t
             t,x  t            t,x           t,x            3t

             3                4              3             2  2        3  4
 - 32*u    *t *x - 48*u  *u *t  - 32*u  *u *t *x - 24*u  *t *x  - 16*u  *t
       2t,x            2t  t          2t  x            2t             t

        2     3          2  2  2             2                3            2
 - 32*u  *u *t *x - 24*u  *t *x  + 24*u *u *t *x - 8*u *u *t*x  + 24*u *t*x
       t   x            t              t  x           t  x            t

       4                     3              2
 - u *x  - 24*u *t*x + 6*u *x  - 30*t - 15*x
    t          x          x

-------- 12. Symmetry:

xi_t=0

xi_x=0

                                         2                            2
eta_u= - 8*u   *u *t*x + u   *x - 4*u  *t  - 4*u    *t*x - 12*u  *u *t
            t,x  t        t,x        3t         2t,x           2t  t

                       2       3  2       2            2  2
 - 4*u  *u *t*x - u  *x  - 4*u  *t  - 4*u  *u *t*x - u  *x  + u *u *x
      2t  x        2t         t          t   x        t        t  x

-------- 13. Symmetry:

xi_t=0

xi_x=0

                            3
eta_u= - u   - 3*u  *u  - u
          3t      2t  t    t

-------- 14. Symmetry:

xi_t=0

xi_x=0

                                2
eta_u= - u   *x - 2*u  *t - 2*u  *t - u *u *x
          t,x        2t        t       t  x

-------- 15. Symmetry:

xi_t=0

xi_x=0

                  2            3       2  3             2             2
eta_u= - 12*u   *t *x - 8*u  *t  - 8*u  *t  - 12*u *u *t *x - 6*u *t*x
             t,x           2t         t           t  x           t

                  3            2
 + 6*u *t*x - u *x  + 6*t + 3*x
      x        x

-------- 16. Symmetry:

xi_t=0

xi_x=0

                  3             4        2  4             3            3
eta_u= - 32*u   *t *x - 16*u  *t  - 16*u  *t  - 32*u *u *t *x - 48*u *t
             t,x            2t          t           t  x            t

          2  2          2             3       2         2    4
 - 24*u *t *x  - 48*u *t *x - 8*u *t*x  - 12*t  - 12*t*x  - x
       t             x           x

-------- 17. Symmetry:

xi_t=0

xi_x=0

                             2       2  2                       2
eta_u= - 8*u   *t*x - 8*u  *t  - 8*u  *t  - 8*u *u *t*x - 2*u *x  + 2*u *x - 1
            t,x          2t         t          t  x          t         x

-------- 18. Symmetry:

xi_t=0

xi_x=0

              2
eta_u=u   + u
       2t    t

-------- 19. Symmetry:

xi_t=0

xi_x=0

eta_u= - 2*u   *t - 2*u *u *t - u *x
            t,x        t  x      t

-------- 20. Symmetry:

xi_t=0

xi_x=0

                 3            3          2            2           2            3
eta_u= - 8*u   *t  - 8*u *u *t  - 12*u *t *x - 12*u *t  - 6*u *t*x  - 6*t*x - x
            t,x         t  x          t            x         x

-------- 21. Symmetry:

xi_t=0

xi_x=0

                 2            2                  2
eta_u= - 4*u   *t  - 4*u *u *t  - 4*u *t*x - u *x  + x
            t,x         t  x         t        x

-------- 22. Symmetry:

xi_t=0

xi_x=0

eta_u=u    + u *u
       t,x    t  x

-------- 23. Symmetry:

xi_t=0

xi_x=0

eta_u= - 2*u *t - x
            x

-------- 24. Symmetry:

xi_t=0

xi_x=0

eta_u=u
       x

-------- 25. Symmetry:

xi_t=0

xi_x=0

eta_u= - 4*u *t - 2*u *x + 1
            t        x

-------- 26. Symmetry:

xi_t=0

xi_x=0

               2                     2
eta_u= - 4*u *t  - 4*u *t*x - 2*t - x
            t         x

-------- 27. Symmetry:

xi_t=0

xi_x=0

eta_u=u
       t

-------- 28. Symmetry:

xi_t=0

xi_x=0

eta_u=1

-------- 

Further symmetries:

xi_t=0

xi_x=0

       c_94
eta_u=------
         u
        e

with c_94(t,x)
which still have to satisfy: 

0= - c_94  + c_94
         t       2x


{{ - c_94  + c_94  },
         t       2x

 {xi_t=0,

  xi_x=0,

               u                    u                2       u                4
  eta_u=( - 2*e *u   *u *c_100 - 8*e *u   *u *c_101*t  - 32*e *u   *u *c_102*t
                  t,x  t               t,x  t                   t,x  t

                u                5       u                3
          - 64*e *u   *u *c_103*t  - 32*e *u   *u *c_104*t
                   t,x  t                   t,x  t

               u                      u
          - 8*e *u   *u *c_105*t - 8*e *u   *u *c_107*t*x
                  t,x  t                 t,x  t

                u                3          u                5
          - 64*e *u   *u *c_108*t *x - 384*e *u   *u *c_109*t *x
                   t,x  t                      t,x  t

                 u                4         u                2
          - 160*e *u   *u *c_110*t *x - 48*e *u   *u *c_111*t *x
                    t,x  t                     t,x  t

               u                      u                 u             2
          - 2*e *u   *u *c_112*x + 2*e *u   *c_101*t - e *u   *c_101*x
                  t,x  t                 t,x               t,x

                u             3       u             2  2        u             4
          + 16*e *u   *c_102*t  - 24*e *u   *c_102*t *x  - 160*e *u   *c_103*t
                   t,x                   t,x                       t,x

                u             3  2       u             2       u               2
          - 80*e *u   *c_103*t *x  + 12*e *u   *c_104*t  - 12*e *u   *c_104*t*x
                   t,x                      t,x                   t,x

             u               u                    u             2
          + e *u   *c_105 + e *u   *c_107*x + 24*e *u   *c_108*t *x
                t,x             t,x                  t,x

               u               3        u             4
          - 8*e *u   *c_108*t*x  - 960*e *u   *c_109*t *x
                  t,x                      t,x

                 u             3  3       u             3
          - 160*e *u   *c_109*t *x  + 80*e *u   *c_110*t *x
                    t,x                      t,x

                u             2  3       u                     u             3
          - 40*e *u   *c_110*t *x  + 12*e *u   *c_111*t*x - 2*e *u   *c_111*x
                   t,x                      t,x                   t,x

               u                      u             3         u             2
          - 8*e *u   *c_115*t*x - 32*e *u   *c_116*t *x - 12*e *u   *c_117*t *x
                  t,x                    t,x                     t,x

             u                 u                 u             2
          - e *u   *c_118*x + e *u   *c_119 - 4*e *u   *c_120*t
                t,x               t,x               t,x

               u             3      u                 u
          - 8*e *u   *c_121*t  - 2*e *u   *c_122*t - e *u  *c_106
                  t,x                  t,x               3t

               u            2       u            4       u            6
          - 4*e *u  *c_107*t  - 16*e *u  *c_108*t  - 64*e *u  *c_109*t
                  3t                   3t                   3t

                u            5       u            3      u
          - 32*e *u  *c_110*t  - 16*e *u  *c_111*t  - 2*e *u  *c_112*t
                   3t                   3t                  3t

             u                  u              2       u              4
          - e *u    *c_100 - 4*e *u    *c_101*t  - 16*e *u    *c_102*t
                2t,x               2t,x                   2t,x

                u              5       u              3      u
          - 32*e *u    *c_103*t  - 16*e *u    *c_104*t  - 4*e *u    *c_105*t
                   2t,x                   2t,x                  2t,x

               u                       u              3
          - 4*e *u    *c_107*t*x - 32*e *u    *c_108*t *x
                  2t,x                    2t,x

                 u              5         u              4
          - 192*e *u    *c_109*t *x - 80*e *u    *c_110*t *x
                    2t,x                     2t,x

                u              2      u                    u
          - 24*e *u    *c_111*t *x - e *u    *c_112*x - 3*e *u  *u *c_106
                   2t,x                  2t,x                 2t  t

                u               2       u               4
          - 12*e *u  *u *c_107*t  - 48*e *u  *u *c_108*t
                   2t  t                   2t  t

                 u               6       u               5
          - 192*e *u  *u *c_109*t  - 96*e *u  *u *c_110*t
                    2t  t                   2t  t

                u               3      u                   u
          - 48*e *u  *u *c_111*t  - 6*e *u  *u *c_112*t - e *u  *u *c_100
                   2t  t                  2t  t               2t  x

               u               2       u               4       u               5
          - 4*e *u  *u *c_101*t  - 16*e *u  *u *c_102*t  - 32*e *u  *u *c_103*t
                  2t  x                   2t  x                   2t  x

                u               3      u                     u
          - 16*e *u  *u *c_104*t  - 4*e *u  *u *c_105*t - 4*e *u  *u *c_107*t*x
                   2t  x                  2t  x                 2t  x

                u               3          u               5
          - 32*e *u  *u *c_108*t *x - 192*e *u  *u *c_109*t *x
                   2t  x                      2t  x

                u               4         u               2
          - 80*e *u  *u *c_110*t *x - 24*e *u  *u *c_111*t *x
                   2t  x                     2t  x

             u                     u                     u            3
          - e *u  *u *c_112*x - 4*e *u  *c_101*t*x - 32*e *u  *c_102*t *x
                2t  x                 2t                    2t

                u            4         u            2        u
          - 80*e *u  *c_103*t *x - 24*e *u  *c_104*t *x - 2*e *u  *c_105*x
                   2t                     2t                    2t

             u            2       u            2  2        u            5
          - e *u  *c_107*x  - 24*e *u  *c_108*t *x  - 480*e *u  *c_109*t
                2t                   2t                       2t

                 u            4  2       u            3  2
          - 240*e *u  *c_109*t *x  - 80*e *u  *c_110*t *x
                    2t                      2t

                u              2    u                u            2
          - 12*e *u  *c_111*t*x  + e *u  *c_114 - 8*e *u  *c_115*t
                   2t                  2t               2t

                u            4      u            3      u
          - 16*e *u  *c_116*t  - 8*e *u  *c_117*t  - 2*e *u  *c_118*t
                   2t                  2t                  2t

             u   3            u   3        2       u   3        4
          - e *u  *c_106 - 4*e *u  *c_107*t  - 16*e *u  *c_108*t
                t                t                    t

                u   3        6       u   3        5       u   3        3
          - 64*e *u  *c_109*t  - 32*e *u  *c_110*t  - 16*e *u  *c_111*t
                   t                    t                    t

               u   3            u   2               u   2           2
          - 2*e *u  *c_112*t - e *u  *u *c_100 - 4*e *u  *u *c_101*t
                  t                t   x               t   x

                u   2           4       u   2           5
          - 16*e *u  *u *c_102*t  - 32*e *u  *u *c_103*t
                   t   x                   t   x

                u   2           3      u   2                 u   2
          - 16*e *u  *u *c_104*t  - 4*e *u  *u *c_105*t - 4*e *u  *u *c_107*t*x
                   t   x                  t   x                 t   x

                u   2           3          u   2           5
          - 32*e *u  *u *c_108*t *x - 192*e *u  *u *c_109*t *x
                   t   x                      t   x

                u   2           4         u   2           2
          - 80*e *u  *u *c_110*t *x - 24*e *u  *u *c_111*t *x
                   t   x                     t   x

             u   2                 u   2                 u   2        3
          - e *u  *u *c_112*x - 4*e *u  *c_101*t*x - 32*e *u  *c_102*t *x
                t   x                 t                     t

                u   2        4         u   2        2        u   2
          - 80*e *u  *c_103*t *x - 24*e *u  *c_104*t *x - 2*e *u  *c_105*x
                   t                      t                     t

             u   2        2       u   2        2  2        u   2        5
          - e *u  *c_107*x  - 24*e *u  *c_108*t *x  - 480*e *u  *c_109*t
                t                    t                        t

                 u   2        4  2       u   2        3  2
          - 240*e *u  *c_109*t *x  - 80*e *u  *c_110*t *x
                    t                       t

                u   2          2    u   2            u   2        2
          - 12*e *u  *c_111*t*x  + e *u  *c_114 - 8*e *u  *c_115*t
                   t                   t                t

                u   2        4      u   2        3      u   2
          - 16*e *u  *c_116*t  - 8*e *u  *c_117*t  - 2*e *u  *c_118*t
                   t                   t                   t

               u                  u              2       u              3
          + 2*e *u *u *c_101*t - e *u *u *c_101*x  + 16*e *u *u *c_102*t
                  t  x               t  x                   t  x

                u              2  2        u              4
          - 24*e *u *u *c_102*t *x  - 160*e *u *u *c_103*t
                   t  x                       t  x

                u              3  2       u              2
          - 80*e *u *u *c_103*t *x  + 12*e *u *u *c_104*t
                   t  x                      t  x

                u                2    u                u
          - 12*e *u *u *c_104*t*x  + e *u *u *c_105 + e *u *u *c_107*x
                   t  x                  t  x             t  x

                u              2        u                3
          + 24*e *u *u *c_108*t *x - 8*e *u *u *c_108*t*x
                   t  x                    t  x

                 u              4          u              3  3
          - 960*e *u *u *c_109*t *x - 160*e *u *u *c_109*t *x
                    t  x                      t  x

                u              3         u              2  3
          + 80*e *u *u *c_110*t *x - 40*e *u *u *c_110*t *x
                   t  x                     t  x

                u                      u              3      u
          + 12*e *u *u *c_111*t*x - 2*e *u *u *c_111*x  - 8*e *u *u *c_115*t*x
                   t  x                   t  x                  t  x

                u              3         u              2      u
          - 32*e *u *u *c_116*t *x - 12*e *u *u *c_117*t *x - e *u *u *c_118*x
                   t  x                     t  x                  t  x

             u                  u              2      u              3
          + e *u *u *c_119 - 4*e *u *u *c_120*t  - 8*e *u *u *c_121*t
                t  x               t  x                  t  x

               u                    u                  u           2
          - 2*e *u *u *c_122*t + 2*e *u *c_101*x + 48*e *u *c_102*t *x
                  t  x                 t                  t

               u             3        u           3         u           2  3
          - 8*e *u *c_102*t*x  - 240*e *u *c_103*t *x - 40*e *u *c_103*t *x
                  t                      t                     t

                u                   u           3       u             2
          + 24*e *u *c_104*t*x - 2*e *u *c_104*x  + 24*e *u *c_108*t*x
                   t                   t                   t

             u           4        u           4        u           3  2
          - e *u *c_108*x  - 720*e *u *c_109*t  - 720*e *u *c_109*t *x
                t                    t                    t

                u           2  4        u           3        u           2  2
          - 60*e *u *c_109*t *x  + 360*e *u *c_110*t  + 120*e *u *c_110*t *x
                   t                       t                    t

                u             4      u           2      u           2
          - 10*e *u *c_110*t*x  + 6*e *u *c_111*x  - 2*e *u *c_115*x
                   t                    t                  t

                u           3       u           2  2      u             2
          - 48*e *u *c_116*t  - 24*e *u *c_116*t *x  - 6*e *u *c_117*t*x
                   t                   t                     t

               u                    u           2      u
          - 4*e *u *c_120*t*x - 12*e *u *c_121*t *x - e *u *c_122*x
                  t                    t                  t

             u               u           2      u                  u           2
          + e *u *c_125 - 4*e *u *c_126*t  - 4*e *u *c_127*t + 84*e *u *c_102*t
                t               t                  t                  x

                u             2    u           4        u           3
          + 36*e *u *c_102*t*x  - e *u *c_102*x  - 120*e *u *c_103*t
                   x                  x                    x

                 u           2  2       u             4      u           2
          - 120*e *u *c_103*t *x  - 10*e *u *c_103*t*x  + 9*e *u *c_104*x
                    x                      x                    x

                u                   u           3        u           3
          - 24*e *u *c_108*t*x + 6*e *u *c_108*x  - 720*e *u *c_109*t *x
                   x                   x                    x

                 u           2  3       u             5        u           2
          - 240*e *u *c_109*t *x  - 12*e *u *c_109*t*x  + 420*e *u *c_110*t *x
                    x                      x                      x

                u             3    u           5      u
          + 60*e *u *c_110*t*x  - e *u *c_110*x  - 6*e *u *c_111*x
                   x                  x                  x

               u                  u           2        u             3
          + 2*e *u *c_115*x - 48*e *u *c_116*t *x - 8*e *u *c_116*t*x
                  x                  x                    x

               u                 u           3    u           2
          + 6*e *u *c_117*t*x - e *u *c_117*x  - e *u *c_120*x
                  x                 x                x

                u           2      u             2    u
          - 12*e *u *c_121*t  - 6*e *u *c_121*t*x  + e *u *c_123
                   x                  x                  x

               u                 u                   u
          - 2*e *u *c_124*t - 4*e *u *c_126*t*x - 2*e *u *c_127*x
                  x                 x                   x

                u                u        3       u        2
          + 48*e *c_102*t*x + 8*e *c_102*x  - 60*e *c_103*t *x

                u          3    u        5      u               u
          - 20*e *c_103*t*x  - e *c_103*x  - 9*e *c_104*x - 30*e *c_108*t

                u        2        u        3        u        2  2
          - 15*e *c_108*x  - 120*e *c_109*t  - 180*e *c_109*t *x

                u          4    u        6        u        2        u          2
          - 30*e *c_109*t*x  - e *c_109*x  + 120*e *c_110*t  + 120*e *c_110*t*x

                u        4      u          u             u        2
          + 10*e *c_110*x  + 3*e *c_111 - e *c_115 - 12*e *c_116*t

                u          2    u        4      u              u        2
          - 12*e *c_116*t*x  - e *c_116*x  + 6*e *c_117*t + 3*e *c_117*x

             u              u              u        3    u
          + e *c_120*x - 6*e *c_121*t*x - e *c_121*x  - e *c_124*x

               u            u        2    u          u                u
          - 2*e *c_126*t - e *c_126*x  + e *c_127 + e *c_128 + c_94)/e },

 {c_94,

  c_105,

  c_104,

  c_103,

  c_102,

  c_101,

  c_100,

  c_112,

  c_111,

  c_110,

  c_109,

  c_108,

  c_107,

  c_106,

  c_118,

  c_117,

  c_116,

  c_115,

  c_114,

  c_122,

  c_121,

  c_120,

  c_119,

  c_124,

  c_123,

  c_127,

  c_126,

  c_125,

  c_128}}

nodepnd {u}$



comment
-------------------------------------------------------
The following example includes the Karpman equations
for three unknown functions in 4 variables. 

If point symmetries are to be computed for a single
equation or a system of equations of higher than first
order then there is the option to formulate at first
preliminary conditions for each equation, have CRACK
solving these conditions before the full set of conditions
is formulated and solved. This strategy is adopted if a
lisp flag prelim_ has the value t. The default value
is nil. 

Similarly, if a system of equations is to be investigated
and a flag individual_ has the value t then symmetry
conditions are formulated and investigated for each
individual equation successively. The default value is nil.

It is advantageous to split a large set of conditions
into smaller sets to be investigated successively if
each set is sufficiently overdetermined to be solvable
quickly. Then any substitutions are done in the smaller
set and the next set of conditions is shorter. For
example, for the Karpman equations below the speedup for
prelim_:=t and individual_:=t is a factor of 10.
(Time ~ 1 min.);

write"-------------------------------------------------------";


-------------------------------------------------------


lisp(prelim_:=t)$


lisp(individual_:=t)$



depend r,x,y,z,t;


depend f,x,y,z,t;


depend v,x,y,z,t;



on time$


Time: 8610 ms  plus GC time: 150 ms

liepde({

first 
solve(
        {df(r,t) + w1*df(r,z)
         + s1*(df(r,x)*df(f,x)+df(r,y)*df(f,y)+r*df(f,x,2)/2+r*df(f,y,2)/2) 
         + s2*(df(r,z)*df(f,z)+r*df(f,z,2)/2),
 
         df(f,t) + w1*df(f,z) 
         - (s1*(df(r,x,2)/r+df(r,y,2)/r-df(f,x)**2-df(f,y)**2) +
            s2*(df(r,z,2)/r-df(f,z)**2))/2 + a1*v,
       
         df(v,t,2) - w2**2*(df(v,x,2)+df(v,y,2)+df(v,z,2))
         - 2*a2*r*(df(r,x,2)+df(r,y,2)+df(r,z,2))
         - 2*a2*(df(r,x)**2+df(r,y)**2+df(r,z)**2)},
         
        {df(v,x,2), df(r,x,2), df(f,x,2)}         
 
     )         
         , {r,f,v}, {x,y,z,t}},

       {"point"}, 
       
       {},{})$



The PDE-system under investigation is :

                 2       2     2          2     2          2     2
v  =( - 4*f *a2*r  - 2*f  *a2*r *s1 - 2*f  *a2*r *s1 - 2*f  *a2*r *s2
 2x        t            x                y                z

                 2          2             2
      - 4*f *a2*r *w1 - 2*r  *a2*s1 - 2*r  *a2*s1 - 2*r  *a2*r*s1
           z               x             y             2z

                            2                           2            2
      + 2*r  *a2*r*s2 - 2*r  *a2*s1 + v  *s1 - v  *s1*w2  - v  *s1*w2
           2z              z           2t       2y           2z

                 2          2
      - 4*a1*a2*r *v)/(s1*w2 )

                2          2          2
r  =(2*f *r + f  *r*s1 + f  *r*s1 + f  *r*s2 + 2*f *r*w1 - r  *s1 - r  *s2
 2x     t      x          y          z            z         2y       2z

      + 2*a1*r*v)/s1

f  =
 2x

  - 2*f *r *s1 - f  *r*s1 - 2*f *r *s1 - f  *r*s2 - 2*f *r *s2 - 2*r  - 2*r *w1
       x  x       2y           y  y       2z           z  z         t      z
--------------------------------------------------------------------------------
                                      r*s1

for the function(s) : 

r(z,y,t,x), f(z,y,t,x), v(z,y,t,x)




=============== Initializations
time for initializations: 30 ms    GC time : 0 ms


=============== Preconditions for the 1. equation
time to formulate conditions: 260 ms    GC time : 20 ms
CRACK needed :  721 ms    GC time : 10 ms


=============== Preconditions for the 2. equation


=============== Preconditions for the 3. equation
time to formulate conditions: 130 ms    GC time : 0 ms
CRACK needed :  230 ms    GC time : 0 ms


=============== Full conditions for the 1. equation
time to formulate conditions: 50 ms    GC time : 0 ms
CRACK needed :  1150 ms    GC time : 20 ms


=============== Full conditions for the 2. equation
time to formulate conditions: 20 ms    GC time : 0 ms
CRACK needed :  70 ms    GC time : 0 ms


=============== Full conditions for the 3. equation
time to formulate conditions: 40 ms    GC time : 0 ms
CRACK needed :  50 ms    GC time : 0 ms

The symmetries are:

-------- 1. Symmetry:

xi_x=0

xi_y=0

xi_z=0

xi_t=0

eta_r=0

eta_f=1

eta_v=0

-------- 2. Symmetry:

xi_x=0

xi_y=0

xi_z=0

xi_t=0

eta_r=0

          2
eta_f= - t

       2*t
eta_v=-----
       a1

-------- 3. Symmetry:

xi_x=0

xi_y=0

xi_z=0

xi_t=0

eta_r=0

eta_f= - t

       1
eta_v=----
       a1

-------- 4. Symmetry:

xi_x=0

      1
xi_y=----
      s2

xi_z=0

xi_t=0

eta_r=0

eta_f=0

eta_v=0

-------- 5. Symmetry:

      1
xi_x=----
      s2

xi_y=0

xi_z=0

xi_t=0

eta_r=0

eta_f=0

eta_v=0

-------- 6. Symmetry:

      y
xi_x=----
      s2

       - x
xi_y=------
       s2

xi_z=0

xi_t=0

eta_r=0

eta_f=0

eta_v=0

-------- 7. Symmetry:

xi_x=0

xi_y=0

      1
xi_z=----
      s1

xi_t=0

eta_r=0

eta_f=0

eta_v=0

-------- 8. Symmetry:

xi_x=0

xi_y=0

xi_z=0

xi_t=1

eta_r=0

eta_f=0

eta_v=0

-------- 


{{},

        c_211 + c_212*y
 {xi_x=-----------------,
              s2

         - c_212*x + c_219
  xi_y=--------------------,
                s2

        c_208
  xi_z=-------,
         s1

  xi_t=c_221,

  eta_r=0,

                       2
  eta_f=c_222 - c_223*t  - c_224*t,

         2*c_223*t + c_224
  eta_v=-------------------},
                a1

 {c_222,

  c_223,

  c_224,

  c_219,

  c_211,

  c_212,

  c_208,

  c_221}}

Time: 2931 ms  plus GC time: 50 ms


off time$


nodepnd {r,f,v}$



comment
-------------------------------------------------------
In the following example a system of two equations (by
V.Sokolov) is investigated concerning a special ansatz for
4th order symmetries. The ansatz for the symmetries includes
two unknown functions f,g. Because x is the second variable
in the list of variables {t,x}, the name u!`2 stands for
df(u,x).
Because higher order symmetries are investigated we have
to set prelim_:=nil. The symmetries to be calculated are
lengthy and therefore conditions are not very overdetermined.
In that case CRACK can take long to solve a single 
subset of conditions. The complete set of conditions would
have been more overdetermined and easier to solve. Therefore
the advantage of first formulating all conditions and then
solving them together with one CRACK call is that having
more equations, the chance of finding short integrable
equations among then is higher, i.e. CRACK has more freedom
in optimizing the computation. Therefore individual_:=nil
is more appropriate in this example.

Because 4th order conditions are to be computed the
`binding stack size' is increased.
(Time ~ 5 min.);

write"-------------------------------------------------------";


-------------------------------------------------------


lisp(prelim_:=nil)$


lisp(individual_:=nil)$


lisp(if getd 'set_bndstk_size then set_bndstk_size(7000))$



nodepnd {u,v}$


depend  u,x,t;


depend  v,x,t;



des:={df(u,t)=+df(u,x,2) + (u + v)*df(u,x) + u*df(v,x),
      df(v,t)=-df(v,x,2) + (u + v)*df(v,x) + v*df(u,x)
     }$



nodepnd {f,g}$


depend f,t,x,u,v,u!`2,v!`2,u!`2!`2,v!`2!`2,u!`2!`2!`2,v!`2!`2!`2$


depend g,t,x,u,v,u!`2,v!`2,u!`2!`2,v!`2!`2,u!`2!`2!`2,v!`2!`2!`2$


liepde({des,{u,v},{t,x}},
       {xi_t=0,
        xi_x=0,
        eta_u=+df(u,x,4)+f,
        eta_v=-df(v,x,4)+g
       },
       {f,g},{}
      )$



The PDE-system under investigation is :

u =u   + u *u + u *v + v *u
 t  2x    x      x      x

v =u *v - v   + v *u + v *v
 t  x      2x    x      x

for the function(s) : 

u(x,t), v(x,t)



The symmetries are:

-------- 1. Symmetry:

xi_t=0

xi_x=0

eta_u=u   + 4*u  *t + 2*u  *u + 2*u  *v + 5*u  *u  + 3*u  *v  + 6*u  *t*u
       4x      3x        3x        3x        2x  x      2x  x      2x

                3       2    9             3       2                 2
 + 6*u  *t*v + ---*u  *u  + ---*u  *u*v + ---*u  *v  + 2*u  *x + 6*u  *t
      2x        2   2x       2   2x        2   2x         2x        x

       2      9    2                             9
 + 3*u  *u + ---*u  *v + 2*u *v   + 6*u *v *t + ---*u *v *u + 3*u *v *v
      x       2   x         x  2x      x  x      2   x  x        x  x

           2                         2    1      3    9      2      9        2
 + 3*u *t*u  + 12*u *t*u*v + 3*u *t*v  + ---*u *u  + ---*u *u *v + ---*u *u*v
      x            x            x         2   x       2   x         2   x

               1      3                                     2
 + 2*u *u*x + ---*u *v  + 2*u *v*x + 8*u  + v  *u + 6*v *t*u  + 6*v *t*u*v
      x        2   x         x          x    3x        x           x

    3      3    9      2      3        2               2
 + ---*v *u  + ---*v *u *v + ---*v *u*v  + 2*v *u*x + u  + 3*u*v
    2   x       2   x         2   x           x

                                                             9
eta_v=u  *v + 2*u  *v  + 3*u *v   - 6*u *v *t - 3*u *v *u - ---*u *v *v
       3x        2x  x      x  2x      x  x        x  x      2   x  x

                        2    3      2      9        2    3      3
 + 6*u *t*u*v + 6*u *t*v  + ---*u *u *v + ---*u *u*v  + ---*u *v  + 2*u *v*x
      x            x         2   x         2   x         2   x         x

 - v   + 4*v  *t + 2*v  *u + 2*v  *v + 5*v  *v  - 6*v  *t*u - 6*v  *t*v
    4x      3x        3x        3x        2x  x      2x          2x

    3       2    9             3       2                 2      9    2
 - ---*v  *u  - ---*v  *u*v - ---*v  *v  - 2*v  *x - 6*v  *t - ---*v  *u
    2   2x       2   2x        2   2x         2x        x       2   x

       2             2                         2    1      3    9      2
 - 3*v  *v + 3*v *t*u  + 12*v *t*u*v + 3*v *t*v  + ---*v *u  + ---*v *u *v
      x         x            x            x         2   x       2   x

    9        2               1      3                       2
 + ---*v *u*v  + 2*v *u*x + ---*v *v  + 2*v *v*x + 3*u*v + v
    2   x           x        2   x         x

-------- 2. Symmetry:

xi_t=0

xi_x=0

                                                               3       2
eta_u=u   + 2*u  *u + 2*u  *v + 4*u   + 5*u  *u  + 3*u  *v  + ---*u  *u
       4x      3x        3x        3x      2x  x      2x  x    2   2x

    9                       3       2                 2      9    2         2
 + ---*u  *u*v + 6*u  *u + ---*u  *v  + 6*u  *v + 3*u  *u + ---*u  *v + 6*u
    2   2x          2x      2   2x         2x        x       2   x         x

               9                                   1      3    9      2
 + 2*u *v   + ---*u *v *u + 3*u *v *v + 6*u *v  + ---*u *u  + ---*u *u *v
      x  2x    2   x  x        x  x        x  x    2   x       2   x

         2    9        2                1      3         2            3      3
 + 3*u *u  + ---*u *u*v  + 12*u *u*v + ---*u *v  + 3*u *v  + v  *u + ---*v *u
      x       2   x            x        2   x         x       3x      2   x

    9      2           2    3        2
 + ---*v *u *v + 6*v *u  + ---*v *u*v  + 6*v *u*v
    2   x           x       2   x           x

                                                 9
eta_v=u  *v + 2*u  *v  + 3*u *v   - 3*u *v *u - ---*u *v *v - 6*u *v
       3x        2x  x      x  2x      x  x      2   x  x        x  x

    3      2      9        2               3      3         2
 + ---*u *u *v + ---*u *u*v  + 6*u *u*v + ---*u *v  + 6*u *v  - v   + 2*v  *u
    2   x         2   x           x        2   x         x       4x      3x

                                 3       2    9                       3       2
 + 2*v  *v + 4*v   + 5*v  *v  - ---*v  *u  - ---*v  *u*v - 6*v  *u - ---*v  *v
      3x        3x      2x  x    2   2x       2   2x          2x      2   2x

              9    2         2         2    1      3    9      2           2
 - 6*v  *v - ---*v  *u - 3*v  *v - 6*v   + ---*v *u  + ---*v *u *v + 3*v *u
      2x      2   x         x         x     2   x       2   x           x

    9        2                1      3         2
 + ---*v *u*v  + 12*v *u*v + ---*v *v  + 3*v *v
    2   x            x        2   x         x

-------- 3. Symmetry:

xi_t=0

xi_x=0

                                                                 3       2
eta_u=u   + 2*u  *u + 2*u  *v + 5*u  *u  + 3*u  *v  + 2*u  *t + ---*u  *u
       4x      3x        3x        2x  x      2x  x      2x      2   2x

    9             3       2       2      9    2                 9
 + ---*u  *u*v + ---*u  *v  + 3*u  *u + ---*u  *v + 2*u *v   + ---*u *v *u
    2   2x        2   2x         x       2   x         x  2x    2   x  x

                                      1      3    9      2      9        2
 + 3*u *v *v + 2*u *t*u + 2*u *t*v + ---*u *u  + ---*u *u *v + ---*u *u*v
      x  x        x          x        2   x       2   x         2   x

    1      3                              3      3    9      2      3        2
 + ---*u *v  + u *x + v  *u + 2*v *t*u + ---*v *u  + ---*v *u *v + ---*v *u*v
    2   x       x      3x        x        2   x       2   x         2   x

 + u

                                                 9
eta_v=u  *v + 2*u  *v  + 3*u *v   - 3*u *v *u - ---*u *v *v + 2*u *t*v
       3x        2x  x      x  2x      x  x      2   x  x        x

    3      2      9        2    3      3
 + ---*u *u *v + ---*u *u*v  + ---*u *v  - v   + 2*v  *u + 2*v  *v + 5*v  *v
    2   x         2   x         2   x       4x      3x        3x        2x  x

              3       2    9             3       2    9    2         2
 - 2*v  *t - ---*v  *u  - ---*v  *u*v - ---*v  *v  - ---*v  *u - 3*v  *v
      2x      2   2x       2   2x        2   2x       2   x         x

                          1      3    9      2      9        2    1      3
 + 2*v *t*u + 2*v *t*v + ---*v *u  + ---*v *u *v + ---*v *u*v  + ---*v *v
      x          x        2   x       2   x         2   x         2   x

 + v *x + v
    x

-------- 4. Symmetry:

xi_t=0

xi_x=0

                                                       3       2    9
eta_u=u   + 2*u  *u + 2*u  *v + 5*u  *u  + 3*u  *v  + ---*u  *u  + ---*u  *u*v
       4x      3x        3x        2x  x      2x  x    2   2x       2   2x

    3       2             2      9    2                 9
 + ---*u  *v  + u   + 3*u  *u + ---*u  *v + 2*u *v   + ---*u *v *u + 3*u *v *v
    2   2x       2x      x       2   x         x  2x    2   x  x        x  x

    1      3    9      2      9        2           1      3
 + ---*u *u  + ---*u *u *v + ---*u *u*v  + u *u + ---*u *v  + u *v + v  *u
    2   x       2   x         2   x         x      2   x       x      3x

    3      3    9      2      3        2
 + ---*v *u  + ---*v *u *v + ---*v *u*v  + v *u
    2   x       2   x         2   x         x

                                                 9             3      2
eta_v=u  *v + 2*u  *v  + 3*u *v   - 3*u *v *u - ---*u *v *v + ---*u *u *v
       3x        2x  x      x  2x      x  x      2   x  x      2   x

    9        2    3      3
 + ---*u *u*v  + ---*u *v  + u *v - v   + 2*v  *u + 2*v  *v + 5*v  *v
    2   x         2   x       x      4x      3x        3x        2x  x

    3       2    9             3       2          9    2         2      1      3
 - ---*v  *u  - ---*v  *u*v - ---*v  *v  - v   - ---*v  *u - 3*v  *v + ---*v *u
    2   2x       2   2x        2   2x       2x    2   x         x       2   x

    9      2      9        2           1      3
 + ---*v *u *v + ---*v *u*v  + v *u + ---*v *v  + v *v
    2   x         2   x         x      2   x       x

-------- 5. Symmetry:

xi_t=0

xi_x=0

                                                       3       2    9
eta_u=u   + 2*u  *u + 2*u  *v + 5*u  *u  + 3*u  *v  + ---*u  *u  + ---*u  *u*v
       4x      3x        3x        2x  x      2x  x    2   2x       2   2x

    3       2       2      9    2                 9
 + ---*u  *v  + 3*u  *u + ---*u  *v + 2*u *v   + ---*u *v *u + 3*u *v *v
    2   2x         x       2   x         x  2x    2   x  x        x  x

    1      3    9      2      9        2    1      3                 3      3
 + ---*u *u  + ---*u *u *v + ---*u *u*v  + ---*u *v  + u  + v  *u + ---*v *u
    2   x       2   x         2   x         2   x       x    3x      2   x

    9      2      3        2
 + ---*v *u *v + ---*v *u*v
    2   x         2   x

                                                 9             3      2
eta_v=u  *v + 2*u  *v  + 3*u *v   - 3*u *v *u - ---*u *v *v + ---*u *u *v
       3x        2x  x      x  2x      x  x      2   x  x      2   x

    9        2    3      3                                         3       2
 + ---*u *u*v  + ---*u *v  - v   + 2*v  *u + 2*v  *v + 5*v  *v  - ---*v  *u
    2   x         2   x       4x      3x        3x        2x  x    2   2x

    9             3       2    9    2         2      1      3    9      2
 - ---*v  *u*v - ---*v  *v  - ---*v  *u - 3*v  *v + ---*v *u  + ---*v *u *v
    2   2x        2   2x       2   x         x       2   x       2   x

    9        2    1      3
 + ---*v *u*v  + ---*v *v  + v
    2   x         2   x       x

-------- 


{{},

 {xi_t=0,

  xi_x=0,

  eta_u=(2*u   + 8*u  *c_326 + 8*u  *c_328*t + 4*u  *u + 4*u  *v + 10*u  *u
            4x      3x            3x              3x        3x         2x  x

          + 6*u  *v  + 12*u  *c_326*u + 12*u  *c_326*v + 12*u  *c_328*t*u
               2x  x       2x               2x               2x

          + 12*u  *c_328*t*v + 4*u  *c_328*x + 2*u  *c_332 + 4*u  *c_333*t
                2x                2x              2x            2x

                   2                      2        2              2
          + 3*u  *u  + 9*u  *u*v + 3*u  *v  + 12*u  *c_326 + 12*u  *c_328*t
               2x         2x          2x          x              x

                2         2
          + 6*u  *u + 9*u  *v + 4*u *v   + 12*u *v *c_326 + 12*u *v *c_328*t
               x         x         x  2x       x  x             x  x

                                                2
          + 9*u *v *u + 6*u *v *v + 6*u *c_326*u  + 24*u *c_326*u*v
               x  x        x  x        x                x

                        2                 2
          + 6*u *c_326*v  + 6*u *c_328*t*u  + 24*u *c_328*t*u*v
               x               x                  x

                          2
          + 6*u *c_328*t*v  + 4*u *c_328*u*x + 4*u *c_328*v*x + 16*u *c_328
               x                 x                x                 x

          + 2*u *c_332*u + 2*u *c_332*v + 4*u *c_333*t*u + 4*u *c_333*t*v
               x              x              x                x

                                            3         2             2       3
          + 2*u *c_333*x + 2*u *c_334 + u *u  + 9*u *u *v + 9*u *u*v  + u *v
               x              x          x         x           x         x

                                   2                                    2
          + 2*v  *u + 12*v *c_326*u  + 12*v *c_326*u*v + 12*v *c_328*t*u
               3x         x                x                 x

          + 12*v *c_328*t*u*v + 4*v *c_328*u*x + 2*v *c_332*u + 4*v *c_333*t*u
                x                  x                x              x

                  3         2             2            2
          + 3*v *u  + 9*v *u *v + 3*v *u*v  + 2*c_328*u  + 6*c_328*u*v
               x         x           x

          + 2*c_333*u)/2,

  eta_v=(2*u  *v + 4*u  *v  + 6*u *v   - 12*u *v *c_326 - 12*u *v *c_328*t
            3x        2x  x      x  2x       x  x             x  x

                                                                   2
          - 6*u *v *u - 9*u *v *v + 12*u *c_326*u*v + 12*u *c_326*v
               x  x        x  x         x                 x

                                               2
          + 12*u *c_328*t*u*v + 12*u *c_328*t*v  + 4*u *c_328*v*x + 2*u *c_332*v
                x                   x                 x                x

                                   2             2         3
          + 4*u *c_333*t*v + 3*u *u *v + 9*u *u*v  + 3*u *v  - 2*v
               x                x           x           x         4x

          + 8*v  *c_326 + 8*v  *c_328*t + 4*v  *u + 4*v  *v + 10*v  *v
               3x            3x              3x        3x         2x  x

          - 12*v  *c_326*u - 12*v  *c_326*v - 12*v  *c_328*t*u
                2x               2x               2x

          - 12*v  *c_328*t*v - 4*v  *c_328*x - 2*v  *c_332 - 4*v  *c_333*t
                2x                2x              2x            2x

                   2                      2        2              2
          - 3*v  *u  - 9*v  *u*v - 3*v  *v  - 12*v  *c_326 - 12*v  *c_328*t
               2x         2x          2x          x              x

                2         2                 2                                 2
          - 9*v  *u - 6*v  *v + 6*v *c_326*u  + 24*v *c_326*u*v + 6*v *c_326*v
               x         x         x                x                x

                          2                                     2
          + 6*v *c_328*t*u  + 24*v *c_328*t*u*v + 6*v *c_328*t*v
               x                  x                  x

          + 4*v *c_328*u*x + 4*v *c_328*v*x + 2*v *c_332*u + 2*v *c_332*v
               x                x                x              x

                                                                              3
          + 4*v *c_333*t*u + 4*v *c_333*t*v + 2*v *c_333*x + 2*v *c_334 + v *u
               x                x                x              x          x

                  2             2       3                          2
          + 9*v *u *v + 9*v *u*v  + v *v  + 6*c_328*u*v + 2*c_328*v  + 2*c_333*v
               x           x         x

         )/2},

 {c_328,

  c_326,

  c_333,

  c_332,

  c_334}}

nodepnd {f,g}$



end$


Time for test: 20401 ms, plus GC time: 360 ms


REDUCE Historical
REDUCE Sourceforge Project | Historical SVN Repository | GitHub Mirror | SourceHut Mirror | NotABug Mirror | Chisel Mirror | Chisel RSS ]