restart;read("Jets.s");Assignments to partial derivativescoordinates([t,x], [u,v], 8);Example 1put('pd(U,u_x^2)' = -u_x*pd(U,u)); pd(U,u_x^2);pd(U,u_x^3);clear(pds);pd(U,u_x^3);Example 2Resolving Cauchy\342\200\223Riemann conditionsput(pd(U,u) = pd(V,v), pd(V,u) = -pd(U,v));produces the Laplace equationpd(U,u^2);i.e.pd(U,u^2)+pd(U,v^2)=0;clear(pds):Cross-derivatives and compatibility conditionsSubstitutions can contradict one another, generating chaotic and session-dependent results. For instance, the above Cauchy\342\200\223Riemann conditions can be alternatively resolved asput(
'pd(V,v)' = pd(U,u),
'pd(V,u)' = -pd(U,v)); but then the inputpd(V,u*v);may produce different results depending on whether it resulted from the \357\254\201rst or the second substitution.cc();put('pd(U,u^2)' = -pd(U,v^2));cc();Varoderingclear(pds):dependence(U(t, x, u, u_x,u_xx), V(x,u,u_x));unknowns(U,V);S := [U, pd(U,u),pd(U,u^2), pd(U,u*x), pd(U,x), pd(U,t), pd(U,u_x), pd(U,u_xx),
V, pd(V,u), pd(V,x*u), pd(V,x), pd(V,u_x)];E := V-pd(U,u_x);Varordering(degree,function,reverse);sort(S, `Vars/<<`);resolve(E);Varordering(function,degree,reverse);sort(S, `Vars/<<`);resolve(E);Trick (two or more unknowns)dependence(U(t, x,u,v,u_x), V(t, x,u,v,u_x));unknowns(U,V);S := TheExpression;Varordering(degree,function,reverse):run(S);Z := clear(pds);Varordering(function,degree,reverse):run(Z);