restart;
read("Jets.s");
SUtKRVRTfk5ld35DQ35mb3J+TWFwbGV+MTN+YXN+b2Z+MDV+T2N0fjIwMTBHNiI=
`Blimit = 25000 ressize = 1000 putsize = 200 maxsize = 100`
Linear equations in total derivatives
Example: Symmetries of KdV
coordinates([t,x], [u], 5);
NjZJJXVfNXRHNiJJJ3VfNHRfeEdGJEkodV90dHR4eEdGJEkodV90dHh4eEdGJEkndV90XzR4R0YkSSV1XzV4R0YkSSV1XzR0R0YkSSd1X3R0dHhHRiRJJ3VfdHR4eEdGJEkndV90eHh4R0YkSSV1XzR4R0YkSSZ1X3R0dEdGJEkmdV90dHhHRiRJJnVfdHh4R0YkSSZ1X3h4eEdGJEkldV90dEdGJEkldV90eEdGJEkldV94eEdGJEkkdV90R0YkSSR1X3hHRiQ=
equation (u_t = u_xxx + 6*u*u_x);
L0kkdV90RzYiLCZJJnVfeHh4R0YkIiIiKiZJInVHRiRGJ0kkdV94R0YkRiciIic=
Symmetries (more precisely, their generating functions U ) are solutions of the determining equation S = 0 where S is the universal linearization operator:
S := symmetries(u = U);
LCotSSNUREc2IjYkSSJVR0YlSSJ0R0YlIiIiKiZGJ0YpSSR1X3hHRiVGKSEiJyomLUYkNiRGJ0kieEdGJUYpSSJ1R0YlRilGLC1GJDYkRicqJEYwIiIkISIi
The determining equation for generating functions of conservation laws is L = 0, where
L := laws(1 = V);
LCgtSSNUREc2IjYkSSJWR0YlSSJ0R0YlISIiKiZJInVHRiUiIiItRiQ2JEYnSSJ4R0YlRiwiIictRiQ2JEYnKiRGLyIiJEYs
As a smooth function, U can depend only on a \357\254\201nite number of variables, i.e., there must be an upper bound on the jet order of U.
Assume that the bound is \357\254\201ve. A \357\254\201fth-order generating function is one that depends on x, t, u, ux , . . . , u5x .
Accordingly, our assumption is made known to the computer as follows:
dependence(U(x,t,u,u_x,u_xx,u_xxx,u_4x,u_5x));
L0kiVUc2IjwqSSJ0R0YkSSJ1R0YkSSJ4R0YkSSR1X3hHRiRJJXVfeHhHRiRJJnVfeHh4R0YkSSV1XzR4R0YkSSV1XzV4R0Yk
To run in automatic regime, Jets also needs to know what are the unknowns:
unknowns(U);
SSJVRzYi
Now launch the computation:
run(S);
<0> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkldV81eEdGJSIiIyIiIQ==
<0> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV80eEdGJSIiIkkldV81eEdGJUYqIiIh
<0> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkmdV94eHhHRiUiIiJJJXVfNXhHRiVGKiIiIQ==
<0> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV94eEdGJSIiIkkldV81eEdGJUYqIiIh
<1> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkkdV94R0YlIiIiSSV1XzV4R0YlRioiIiE=
<1> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidUdGJSIiIkkldV81eEdGJUYqIiIh
<1> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkieEdGJSIiIkkldV81eEdGJUYqIiIh
<2> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkldV80eEdGJSIiIyIiIQ==
<2> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV80eEdGJSIiIkkmdV94eHhHRiVGKiIiIQ==
<2> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV94eEdGJSIiIkkldV80eEdGJUYqIiIh
<2> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkkdV94R0YlIiIiSSV1XzR4R0YlRioiIiE=
<3> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidUdGJSIiIkkldV80eEdGJUYqIiIh
<3> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV80eEdGJSIiIkkieEdGJUYqIiIh
<3> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkmdV94eHhHRiUiIiMiIiE=
<4> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV94eEdGJSIiIkkmdV94eHhHRiVGKiIiIQ==
<4> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkkdV94R0YlIiIiSSZ1X3h4eEdGJUYqIiIh
<4> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkmdV94eHhHRiUiIiJJInVHRiUiIiMiIiE=
<4> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidUdGJSIiIkkmdV94eHhHRiVGKiwkLUYkNiRGJ0kldV81eEdGJSIjNQ==
<5> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidEdGJSIiIkkldV81eEdGJUYqLCQtRiQ2JEYnKiZJJnVfeHh4R0YlRipJInhHRiVGKiIiJA==
<5> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkldV94eEdGJSIiIyIiIQ==
Ly1JI3BkRzYiNiRJIlVHRiUqJkkmdV94eHhHRiUiIiJJInhHRiUiIiMiIiE=
<6> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV94eEdGJSIiIkkkdV94R0YlIiIjIiIh
<6> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqKEkldV94eEdGJSIiIkkidUdGJUYqSSR1X3hHRiVGKiIiIQ==
<7> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkkdV94R0YlIiIiSSV1X3h4R0YlRiosJC1GJDYkRidJJXVfNXhHRiUiIz8=
<7> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV94eEdGJSIiIkkidUdGJSIiIyIiIQ==
<7> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidUdGJSIiIkkldV94eEdGJUYqLCQtRiQ2JEYnSSV1XzR4R0YlIiIp
<8> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidEdGJSIiIkkldV80eEdGJUYqLCQtRiQ2JEYnKiZJJXVfeHhHRiVGKkkieEdGJUYqIiIk
<8> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkkdV94R0YlIiIkIiIh
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV94eEdGJSIiIkkieEdGJSIiIyIiIQ==
<9> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkkdV94R0YlIiIjSSJ1R0YlIiIiIiIh
<9> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkkdV94R0YlIiIjLCQtRiQ2JEYnSSV1XzR4R0YlIiM3
<10> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidUdGJSIiIkkkdV94R0YlRiosJC1GJDYkRidJJnVfeHh4R0YlIiIn
<10> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkmdV94eHhHRiUiIiJJInhHRiVGKiIiIQ==
<11> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidEdGJSIiIkkmdV94eHhHRiVGKiwkLUYkNiRGJyomSSR1X3hHRiVGKkkieEdGJUYqIiIk
<11> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkidUdGJSIiJSIiIQ==
Ly1JI3BkRzYiNiRJIlVHRiUqJkkieEdGJSIiI0kkdV94R0YlIiIiIiIh
<12> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkidUdGJSIiJCwkLUYkNiRGJ0kldV80eEdGJSIjSw==
<12> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJEkidUdGJSIiIywmKiYtRiQ2JEYnSSV1XzV4R0YlIiIiSSR1X3hHRiVGMCEjPy1GJDYkRidJJXVfeHhHRiUiIiU=
<13> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkldV94eEdGJSIiIkkieEdGJUYqIiIh
<13> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidEdGJSIiIkkldV94eEdGJUYqLCQtRiQ2JEYnKiZJInVHRiVGKkkieEdGJUYqIiIk
<14> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkieEdGJSIiI0kidUdGJSIiIiIiIQ==
Ly1JI3BkRzYiNiRJIlVHRiVJJXVfNHhHRiUiIiE=
<14> Put:
Ly1JI3BkRzYiNiRJIlVHRiVJJXVfeHhHRiUsJComLUYkNiRGJ0kldV81eEdGJSIiIkkkdV94R0YlRi4iIz8=
<15> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkkdV94R0YlIiIiSSJ4R0YlRiosKComLUYkNiRGJ0kmdV94eHhHRiVGKkYpRiohIiQqJi1GJDYkRidJJXVfNXhHRiVGKkYwRiohIiYtRiQ2JEYnSSJ1R0YlI0YqIiIj
<15> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidEdGJSIiIkkkdV94R0YlRiosMEYnIiInKiYtRiQ2JEYnRitGKkYrRiohIicqKEkidUdGJUYqLUYkNiRGJ0kldV81eEdGJUYqSSZ1X3h4eEdGJUYqISNJKiZGNkYqRjRGKkYxKiYtRiQ2JEYnRjdGKkY3RipGMSooRjtGKkYzRipGK0YqISM9KiZGM0YqLUYkNiRGJ0YzRioiIiQ=
<16> Put:
Ly1JI3BkRzYiNiRJIlVHRiUqJkkidUdGJSIiIkkidEdGJUYqLCotRiQ2JEYnSSJ4R0YlIiM3KiZJJHVfeEdGJSIiIy1GJDYkRidJJnVfeHh4R0YlRiohIz0qKC1GJDYkRidJJXVfNXhHRiVGKkY2RipGMkYqISNJKiZGMkYqLUYkNiRGJ0YpRioiIiQ=
Ly1JI3BkRzYiNiRJIlVHRiUqJkkieEdGJSIiIkkidUdGJUYqIiIh
<16> Put:
Ly1JI3BkRzYiNiRJIlVHRiVJInhHRiUsKComSSR1X3hHRiUiIiItRiQ2JEYnSSJ1R0YlRiwjRiwiIiMqJkYrRjEtRiQ2JEYnSSZ1X3h4eEdGJUYsISIkKigtRiQ2JEYnSSV1XzV4R0YlRixGNUYsRitGLCEiJg==
Ly1JI3BkRzYiNiRJIlVHRiUqJEkieEdGJSIiIyIiIQ==
<16> Put:
Ly1JI3BkRzYiNiRJIlVHRiVJInRHRiVRXHY2KlUqdV94LTYqdV94XjIqcGQoVSx1X3gpLTYqcGQoVSx1XzV4KSp1XzV4KnVfeC0zMCpwZChVLHVfNXgpKnVfeCp1KnVfeHh4LTE1KnBkKFUsdV81eCkqdV94eHheMi0xNSp1X3gqdV94eHgqcGQoVSx1X3h4eCktMTgqcGQoVSx1X3h4eCkqdSp1X3heMiszKnBkKFUsdSkqdSp1X3grMy8yKnVfeHh4KnBkKFUsdSlGJQ==
<17> Success!
To retrieve the results, on should proceed as follows.
dependence();
L0kiVUc2IjwpSSJ0R0YkSSJ1R0YkSSJ4R0YkSSR1X3hHRiRJJXVfeHhHRiRJJnVfeHh4R0YkSSV1XzV4R0Yk
Z := clear(pds);
Ny8tSSNwZEc2IjYkSSJVR0YlKiRJJHVfeEdGJSIiIy1GJDYkRicqJEkldV81eEdGJUYqLUYkNiRGJyokSSZ1X3h4eEdGJUYqLUYkNiRGJyomSSJ1R0YlIiIiRi5GNy1GJDYkRicqJkYpRjdGLkY3LUYkNiRGJyomRilGN0YyRjctRiQ2JEYnKiZGMkY3Ri5GNywmLUYkNiRGJyomRjZGN0YpRjdGNy1GJDYkRidGMiEiJywmLUYkNiRGJyomRjZGN0YyRjdGNy1GJDYkRidGLiEjNSwmLUYkNiRGJyokRjZGKkY3KiZGTEY3RilGNyEjZywmLUYkNiRGJ0kldV94eEdGJUY3RlMhIz8sKi1GJDYkRidJInhHRiVGNyomRilGNy1GJDYkRidGNkY3IyEiIkYqKiZGKUYqRkVGNyIiJCooRkxGN0YyRjdGKUY3IiImLC4tRiQ2JEYnSSJ0R0YlRjcqJkYnRjdGKUY3RkcqJkYpRiotRiQ2JEYnRilGNyIiJyomLChGOkZHKihGMkY3RjZGN0YpRjchI0lGMSEjOkY3RkxGN0ZcbyomLCZGPUZecComRjZGN0YpRiohIz1GN0ZFRjdGXG8qJiwmRkRGXm9GMiNGXm9GKkY3RmluRjdGXG8=
The above command also removes all assignments to partial derivatives (otherwise Z would be a set of zeroes). This is why clear(pds) should not be called repeatedly.
First step
U := U1*u_5x + U0;
dependence(U1(t,x,u,u_x,u_xx,u_xxx), U0(t,x,u,u_x,u_xx,u_xxx));
LCYqJkkjVTFHNiIiIiJJJXVfNXhHRiVGJkYmSSNVMEdGJUYm
NiQvSSNVMEc2IjwoSSJ0R0YlSSJ1R0YlSSJ4R0YlSSR1X3hHRiVJJXVfeHhHRiVJJnVfeHh4R0YlL0kjVTFHRiVGJg==
map(expand, Z);
Ny8sJiomLUkjcGRHNiI2JEkjVTFHRicqJEkkdV94R0YnIiIjIiIiSSV1XzV4R0YnRi1GLS1GJjYkSSNVMEdGJ0YqRi0iIiEsJiomLUYmNiRGKSokSSZ1X3h4eEdGJ0YsRi1GLkYtRi0tRiY2JEYxRjdGLS1GJjYkRilJInVHRictRiY2JEYpRissJiomLUYmNiRGKSomRitGLUY4Ri1GLUYuRi1GLS1GJjYkRjFGREYtLUYmNiRGKUY4LCoqJi1GJjYkRikqJkY9Ri1GK0YtRi1GLkYtRi0tRiY2JEYxRk1GLSomRkdGLUYuRi0hIictRiY2JEYxRjhGUSwoKiYtRiY2JEYpKiZGPUYtRjhGLUYtRi5GLUYtLUYmNiRGMUZYRi1GKSEjNSwoKiYtRiY2JEYpKiRGPUYsRi1GLkYtRi0tRiY2JEYxRmpuRi0qJkYpRi1GK0YtISNnLCgqJi1GJjYkRilJJXVfeHhHRidGLUYuRi1GLS1GJjYkRjFGY29GLUZdbyEjPywwKiYtRiY2JEYpSSJ4R0YnRi1GLkYtRi0tRiY2JEYxRltwRi0qKEYrRi1GO0YtRi5GLSMhIiJGLComRitGLS1GJjYkRjFGPUYtRl9wKihGK0YsRkdGLUYuRi0iIiQqJkYrRixGUkYtRmVwKihGKUYtRjhGLUYrRi0iIiYsQComLUYmNiRGKUkidEdGJ0YtRi5GLUYtLUYmNiRGMUZdcUYtKiZGK0YtRjFGLUZRKihGK0YsRj5GLUYuRi0iIicqJkYrRiwtRiY2JEYxRitGLUZicSoqRilGLUY4Ri1GPUYtRitGLSIjSSomRilGLUY4RiwiIzoqKkYrRi1GOEYtRkdGLUYuRi1GaXEqKEYrRi1GOEYtRlJGLUZpcSoqRj1GLUYrRixGR0YtRi5GLSIjPSooRj1GLUYrRixGUkYtRl1yKipGPUYtRitGLUY7Ri1GLkYtISIkKihGPUYtRitGLUZicEYtRmByKihGOEYtRjtGLUYuRi0jRmByRiwqJkY4Ri1GYnBGLUZjcg==
map(expand@pd, Z, u_5x);
Ny8tSSNwZEc2IjYkSSNVMUdGJSokSSR1X3hHRiUiIiMiIiEtRiQ2JEYnKiRJJnVfeHh4R0YlRipGK0YrLUYkNiRGJyomRikiIiJGL0YzRissJi1GJDYkRicqJkkidUdGJUYzRilGM0YzLUYkNiRGJ0YvISInLUYkNiRGJyomRjhGM0YvRjMtRiQ2JEYnKiRGOEYqLUYkNiRGJ0kldV94eEdGJSwoLUYkNiRGJ0kieEdGJUYzKiZGKUYzLUYkNiRGJ0Y4RjMjISIiRioqJkYpRipGOUYzIiIkLC4tRiQ2JEYnSSJ0R0YlRjMqJkYpRiotRiQ2JEYnRilGMyIiJyooRjlGM0YpRjNGL0YzIiM6KihGOUYzRjhGM0YpRioiIz0qKEZKRjNGOEYzRilGMyEiJComRkpGM0YvRjMjRmduRio=
unknowns(U1,U0):
run(Z);
<0> Put:
Ly1JI3BkRzYiNiRJI1UxR0YlSSJ0R0YlIiIh
Ly1JI3BkRzYiNiRJI1UxR0YlSSJ1R0YlIiIh
Ly1JI3BkRzYiNiRJI1UxR0YlSSJ4R0YlIiIh
Ly1JI3BkRzYiNiRJI1UxR0YlKiRJInVHRiUiIiMiIiE=
Ly1JI3BkRzYiNiRJI1UxR0YlKiRJJHVfeEdGJSIiIyIiIQ==
Ly1JI3BkRzYiNiRJI1UxR0YlKiRJJnVfeHh4R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJI1UxR0YlKiZJInVHRiUiIiJJJHVfeEdGJUYqIiIh
Ly1JI3BkRzYiNiRJI1UxR0YlKiZJInVHRiUiIiJJJnVfeHh4R0YlRioiIiE=
Ly1JI3BkRzYiNiRJI1UxR0YlKiZJJHVfeEdGJSIiIkkmdV94eHhHRiVGKiIiIQ==
Ly1JI3BkRzYiNiRJI1UxR0YlSSR1X3hHRiUiIiE=
Ly1JI3BkRzYiNiRJI1UxR0YlSSV1X3h4R0YlIiIh
Ly1JI3BkRzYiNiRJI1UxR0YlSSZ1X3h4eEdGJSIiIQ==
<0> Put:
Ly1JI3BkRzYiNiRJI1UwR0YlSSJ0R0YlUVx0Nip1X3gqVTArMyp1KnVfeCpwZChVMCx1KSszLzIqdV94eHgqcGQoVTAsdSktMTUqVTEqdV94eHheMi0zMCpVMSp1X3h4eCp1KnVfeC0xNSp1X3gqdV94eHgqcGQoVTAsdV94eHgpLTE4KnUqdV94XjIqcGQoVTAsdV94eHgpLTYqdV94XjIqcGQoVTAsdV94KUYl
Ly1JI3BkRzYiNiRJI1UwR0YlSSJ4R0YlLCgqJkkkdV94R0YlIiIiLUYkNiRGJ0kidUdGJUYsI0YsIiIjKiZGK0YxLUYkNiRGJ0kmdV94eHhHRiVGLCEiJCooSSNVMUdGJUYsRjVGLEYrRiwhIiY=
Ly1JI3BkRzYiNiRJI1UwR0YlKiRJJHVfeEdGJSIiIyIiIQ==
Ly1JI3BkRzYiNiRJI1UwR0YlKiRJJnVfeHh4R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJInVHRiUiIiJJJHVfeEdGJUYqLCQtRiQ2JEYnSSZ1X3h4eEdGJSIiJw==
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJJHVfeEdGJSIiIkkmdV94eHhHRiVGKiIiIQ==
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJJXVfeHhHRiUiIiJJJnVfeHh4R0YlRioiIiE=
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJJnVfeHh4R0YlIiIiSSJ1R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJJnVfeHh4R0YlIiIjSSJ1R0YlIiIiIiIh
<0> Put:
Ly1JI3BkRzYiNiRJI1UwR0YlKiRJInVHRiUiIiMsJComSSNVMUdGJSIiIkkkdV94R0YlRi4iI2c=
Ly1JI3BkRzYiNiRJI1UwR0YlKiRJInVHRiUiIiQiIiE=
Ly1JI3BkRzYiNiRJI1UwR0YlKiRJJXVfeHhHRiUiIiMiIiE=
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJInVHRiUiIiJJJXVfeHhHRiVGKiIiIQ==
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJInVHRiUiIiJJJnVfeHh4R0YlRiosJEkjVTFHRiUiIzU=
Ly1JI3BkRzYiNiRJI1UwR0YlKiZJJXVfeHhHRiUiIiJJInVHRiUiIiMiIiE=
<0> Put:
Ly1JI3BkRzYiNiRJI1UwR0YlSSV1X3h4R0YlLCQqJkkjVTFHRiUiIiJJJHVfeEdGJUYsIiM/
<1> Success!
dependence();
NiQvSSNVMEc2IjwoSSJ0R0YlSSJ1R0YlSSJ4R0YlSSR1X3hHRiVJJXVfeHhHRiVJJnVfeHh4R0YlL0kjVTFHRiU8Ig==
Z := clear(pds);
NystSSNwZEc2IjYkSSNVMEdGJSokSSR1X3hHRiUiIiMtRiQ2JEYnKiRJJnVfeHh4R0YlRiotRiQ2JEYnKiZGKSIiIkYuRjIsJi1GJDYkRicqJkkidUdGJUYyRi5GMkYySSNVMUdGJSEjNSwmLUYkNiRGJyokRjdGKkYyKiZGOEYyRilGMiEjZywmLUYkNiRGJ0kldV94eEdGJUYyRj4hIz8sJi1GJDYkRicqJkY3RjJGKUYyRjItRiQ2JEYnRi4hIicsKi1GJDYkRidJInhHRiVGMiomRilGMi1GJDYkRidGN0YyIyEiIkYqKiZGKUYqRklGMiIiJCooRjhGMkYuRjJGKUYyIiImLC4tRiQ2JEYnSSJ0R0YlRjIqJkYpRjJGJ0YyRksqJiwmRkhGVkYuI0ZWRipGMkZRRjJGVComLCZGLSEjOiooRi5GMkY3RjJGKUYyISNJRjJGOEYyRlQqJiwmRjFGXW8qJkY3RjJGKUYqISM9RjJGSUYyRlQqJkYpRiotRiQ2JEYnRilGMiIiJw==
Second step
unksselect(Z, {U0});
NyYtSSNwZEc2IjYkSSNVMEdGJSokSSR1X3hHRiUiIiMtRiQ2JEYnKiRJJnVfeHh4R0YlRiotRiQ2JEYnKiZGKSIiIkYuRjIsJi1GJDYkRicqJkkidUdGJUYyRilGMkYyLUYkNiRGJ0YuISIn
U0 := U03*u_x+U04*u_xxx+U05;
dependence(U03(t, u, x, u_xx), U04(t, u, x, u_xx), U05(t, u, x, u_xx));
LCgqJkkkVTAzRzYiIiIiSSR1X3hHRiVGJkYmKiZJJFUwNEdGJUYmSSZ1X3h4eEdGJUYmRiZJJFUwNUdGJUYm
NiYvSSRVMDNHNiI8JkkidEdGJUkidUdGJUkieEdGJUkldV94eEdGJS9JJFUwNEdGJUYmL0kkVTA1R0YlRiYvSSNVMUdGJTwi
map(simpl, {Diff(U0, u_x, u_x), Diff(U0, u_xxx, u_x), Diff(U0, u_xxx, u_xxx)});
PCMiIiE=
map(simpl,Z);
NysiIiFGI0YjLCYtSSNwZEc2IjYkSSRVMDRHRidJInVHRiciIiJJI1UxR0YnISM1LCoqJi1GJjYkSSRVMDNHRicqJEYqIiIjRitJJHVfeEdGJ0YrRisqJi1GJjYkRilGM0YrSSZ1X3h4eEdGJ0YrRistRiY2JEkkVTA1R0YnRjNGKyomRixGK0Y1RishI2csKiomLUYmNiRGMkkldV94eEdGJ0YrRjVGK0YrKiYtRiY2JEYpRkNGK0Y5RitGKy1GJjYkRjxGQ0YrRj0hIz8sJi1GJjYkRjJGKkYrRikhIicsMiomLUYmNiRGMkkieEdGJ0YrRjVGK0YrKiYtRiY2JEYpRlJGK0Y5RitGKy1GJjYkRjxGUkYrKiZGS0YrRjVGNCMhIiJGNCooRjVGK0YlRitGOUYrRlkqJkY1RistRiY2JEY8RipGK0ZZKiZGNUY0RilGKyIiJCooRixGK0Y5RitGNUYrIiImLD4qJi1GJjYkRjJJInRHRidGK0Y1RitGKyomLUYmNiRGKUZhb0YrRjlGK0YrLUYmNiRGPEZhb0YrKihGNUYrRilGK0Y5RisiIioqJkY1RitGPEYrRk0qKEYqRitGNUY0RktGKyEiJCoqRipGK0Y1RitGJUYrRjlGK0ZbcCooRipGK0Y1RitGZ25GK0ZbcCooRjlGK0ZLRitGNUYrI0ZbcEY0KiZGJUYrRjlGNEZfcComRjlGK0ZnbkYrRl9wKiZGLEYrRjlGNCIjOioqRixGK0Y5RitGKkYrRjVGKyIjSSooRilGK0YqRitGNUY0IiM9
unknowns(U1,U03, U04, U05):
run(Z);
<0> Put:
L0kkVTA0RzYiLCQtSSNwZEdGJDYkSSRVMDNHRiRJInVHRiQjIiIiIiIn
Ly1JI3BkRzYiNiRJJFUwNEdGJSokSSJ1R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJJFUwNEdGJSomSSJ1R0YlIiIiSSV1X3h4R0YlRioiIiE=
Ly1JI3BkRzYiNiRJJFUwNEdGJSomSSV1X3h4R0YlIiIiSSJ4R0YlRioiIiE=
Ly1JI3BkRzYiNiRJJFUwNEdGJUkldV94eEdGJSIiIQ==
<0> Put:
Ly1JI3BkRzYiNiRJJFUwM0dGJSokSSJ1R0YlIiIkIiIh
Ly1JI3BkRzYiNiRJJFUwM0dGJSomSSJ1R0YlIiIiSSV1X3h4R0YlRioiIiE=
Ly1JI3BkRzYiNiRJJFUwM0dGJSomSSV1X3h4R0YlIiIiSSJ1R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJJFUwM0dGJSooSSJ4R0YlIiIiSSJ1R0YlRipJJXVfeHhHRiVGKiIiIQ==
<0> Put:
Ly1JI3BkRzYiNiRJJFUwM0dGJSokSSJ1R0YlIiIjLCRJI1UxR0YlIiNn
Ly1JI3BkRzYiNiRJJFUwNUdGJSokSSV1X3h4R0YlIiIjLCQqJi1GJDYkSSRVMDNHRiVGKCIiIkkkdV94R0YlRjAhIiI=
Ly1JI3BkRzYiNiRJJFUwNUdGJSomSSV1X3h4R0YlIiIiSSJ1R0YlIiIjIiIh
<0> Put:
Ly1JI3BkRzYiNiRJJFUwM0dGJSokSSV1X3h4R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJJFUwM0dGJSomSSJ4R0YlIiIiSSJ1R0YlRioiIiE=
Ly1JI3BkRzYiNiRJJFUwM0dGJUkldV94eEdGJSwkSSNVMUdGJSIjPw==
Ly1JI3BkRzYiNiRJJFUwNUdGJSokSSJ1R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJJFUwNUdGJSomSSJ1R0YlIiIiSSV1X3h4R0YlRioiIiE=
<0> Put:
Ly1JI3BkRzYiNiRJJFUwNUdGJUkidUdGJSwkLUYkNiRJJFUwM0dGJSomRigiIiJJInRHRiVGLiNGLiIiKg==
Ly1JI3BkRzYiNiRJJFUwNUdGJSomSSV1X3h4R0YlIiIiSSJ4R0YlRioiIiE=
Ly1JI3BkRzYiNiRJJFUwNUdGJUkldV94eEdGJSIiIQ==
<0> Put:
L0kkVTA1RzYiLCYtSSNwZEdGJDYkSSRVMDNHRiRJInRHRiQjIiIiIiInKiZJInVHRiRGLC1GJzYkRilJInhHRiRGLCEiIg==
Ly1JI3BkRzYiNiRJJFUwM0dGJSomSSJ1R0YlIiIiSSJ0R0YlRiosJC1GJDYkRidJInhHRiUiIz0=
<0> Put:
Ly1JI3BkRzYiNiRJJFUwM0dGJSokSSJ4R0YlIiIjIiIh
Ly1JI3BkRzYiNiRJJFUwM0dGJSomSSJ0R0YlIiIiSSJ4R0YlRioiIiE=
<0> Put:
Ly1JI3BkRzYiNiRJJFUwM0dGJSokSSJ0R0YlIiIjIiIh
<0> Success!
dependence();
NiQvSSRVMDNHNiI8JkkidEdGJUkidUdGJUkieEdGJUkldV94eEdGJS9JI1UxR0YlPCI=
Third step
Z := clear(pds);
NyktSSNwZEc2IjYkSSRVMDNHRiUqJEkieEdGJSIiIy1GJDYkRicqJEkidEdGJUYqLUYkNiRGJyomRi4iIiJGKUYyLUYkNiRGJyomRilGMkkidUdGJUYyLCZJI1UxR0YlISNnLUYkNiRGJyokRjZGKkYyLCYtRiQ2JEYnSSV1X3h4R0YlRjJGOCEjPywmLUYkNiRGJyomRjZGMkYuRjJGMi1GJDYkRidGKSEjPQ==
unksselect(Z, {U03});
NyctSSNwZEc2IjYkSSRVMDNHRiUqJEkieEdGJSIiIy1GJDYkRicqJEkidEdGJUYqLUYkNiRGJyomRi4iIiJGKUYyLUYkNiRGJyomRilGMkkidUdGJUYyLCYtRiQ2JEYnKiZGNkYyRi5GMkYyLUYkNiRGJ0YpISM9
lprint(map(convert, {Diff(U03, t, t), Diff(U03, x, t), Diff(U03, x, u), Diff(U03, x, x)}, diff));
{diff(U03(t, u, x, u_xx), t, t), diff(U03(t, u, x, u_xx), t, x), diff(U03(t, u, x, u_xx), x, u), diff(U03(t, u, x, u_xx), x, x)}
U03 := U031*t+U032*x+U033;
dependence(U031(u, u_xx),U032(u_xx),U033(u, u_xx));
LCgqJkklVTAzMUc2IiIiIkkidEdGJUYmRiYqJkklVTAzMkdGJUYmSSJ4R0YlRiZGJkklVTAzM0dGJUYm
NiYvSSNVMUc2IjwiL0klVTAzMUdGJTwkSSJ1R0YlSSV1X3h4R0YlL0klVTAzMkdGJTwjRisvSSVVMDMzR0YlRik=
map(simpl, {Diff(U03, t, t), Diff(U03, x, t), Diff(U03, x, u), Diff(U03, x, x)});
PCMiIiE=
map(simpl, Z);
NykiIiFGI0YjRiMsKEkjVTFHNiIhI2cqJi1JI3BkR0YmNiRJJVUwMzFHRiYqJEkidUdGJiIiIyIiIkkidEdGJkYwRjAtRio2JEklVTAzM0dGJkYtRjAsKiomLUYqNiRGLEkldV94eEdGJkYwRjFGMEYwKiYtRio2JEklVTAzMkdGJkY5RjBJInhHRiZGMEYwLUYqNiRGNEY5RjBGJSEjPywmLUYqNiRGLEYuRjBGPSEjPQ==
unknowns(U1, U032, U031, U033);
NiZJI1UxRzYiSSVVMDMyR0YkSSVVMDMxR0YkSSVVMDMzR0Yk
run(Z);
<0> Put:
Ly1JI3BkRzYiNiRJJVUwMzFHRiVJInVHRiUsJEklVTAzMkdGJSIjPQ==
Ly1JI3BkRzYiNiRJJVUwMzFHRiVJJXVfeHhHRiUsJComSSJ1R0YlIiIiLUYkNiRJJVUwMzJHRiVGKEYsIiIn
Ly1JI3BkRzYiNiRJJVUwMzNHRiUqJEkldV94eEdGJSIiIywmKiYtRiQ2JEklVTAzMUdGJUYoIiIiSSJ0R0YlRjAhIiIqJi1GJDYkSSVVMDMyR0YlRihGMEkieEdGJUYwRjI=
Ly1JI3BkRzYiNiRJJVUwMzNHRiUqJkkidUdGJSIiIkkldV94eEdGJUYqLCQqJi1GJDYkSSVVMDMxR0YlRihGKkkidEdGJUYqISIi
Ly1JI3BkRzYiNiRJJVUwMzNHRiUqJkkldV94eEdGJSIiIkkidUdGJSIiIywkKiYtRiQ2JEklVTAzMUdGJUYoRipJInRHRiVGKiEiIg==
<0> Put:
Ly1JI3BkRzYiNiRJJVUwMzJHRiUqJEkldV94eEdGJSIiIyIiIQ==
Ly1JI3BkRzYiNiRJJVUwMzJHRiVJJXVfeHhHRiUiIiE=
Ly1JI3BkRzYiNiRJJVUwMzNHRiUqJEkidUdGJSIiJCIiIQ==
<0> Put:
Ly1JI3BkRzYiNiRJJVUwMzNHRiUqJEkidUdGJSIiIywkSSNVMUdGJSIjZw==
Ly1JI3BkRzYiNiRJJVUwMzNHRiVJJXVfeHhHRiUsJEkjVTFHRiUiIz8=
<0> Success!
dependence();
NiYvSSNVMUc2IjwiL0klVTAzMUdGJTwjSSJ1R0YlL0klVTAzMkdGJUYmL0klVTAzM0dGJTwkRipJJXVfeHhHRiU=
Final step
Z := clear(pds);
NyUsJi1JI3BkRzYiNiRJJVUwMzFHRiZJInVHRiYiIiJJJVUwMzJHRiYhIz0sJkkjVTFHRiYhI2ctRiU2JEklVTAzM0dGJiokRikiIiNGKiwmLUYlNiRGMkkldV94eEdGJkYqRi4hIz8=
lprint(map(convert, Z, diff));
[diff(U031(u), u)-18*U032, -60*U1+diff(U033(u, u_xx), u, u), diff(U033(u, u_xx), u_xx)-20*U1]
U031 := 18*U032*u+C3;
U033 := (20*u_xx+30*u^2)*U1+C1*u+C2;
dependence(C1(), C2(),C3());
LCYqJkkidUc2IiIiIkklVTAzMkdGJUYmIiM9SSNDM0dGJUYm
LCgqJiwmSSV1X3h4RzYiIiM/KiRJInVHRiYiIiMiI0kiIiJJI1UxR0YmRixGLComSSNDMUdGJkYsRilGLEYsSSNDMkdGJkYs
NicvSSNDMUc2IjwiL0kjQzJHRiVGJi9JI0MzR0YlRiYvSSNVMUdGJUYmL0klVTAzMkdGJUYm
map(simpl,Z);
NyUiIiFGI0Yj