restart; read("Jets.s"); SUtKRVRTfk5ld35DQ35mb3J+TWFwbGV+MTN+YXN+b2Z+MDV+T2N0fjIwMTBHNiI= `Blimit = 25000 ressize = 1000 putsize = 200 maxsize = 100`
<Text-field style="Heading 1" layout="Heading 1">Assignments to partial derivatives</Text-field> coordinates([t,x], [u,v], 8); NmRwSSV1Xzh0RzYiSSd1Xzd0X3hHRiRJKHVfNnRfeHhHRiRJKXVfNXRfeHh4R0YkSSh1XzR0XzR4R0YkSSl1X3R0dF81eEdGJEkodV90dF82eEdGJEkndV90Xzd4R0YkSSV1Xzh4R0YkSSV1Xzd0R0YkSSd1XzZ0X3hHRiRJKHVfNXRfeHhHRiRJKXVfNHRfeHh4R0YkSSl1X3R0dF80eEdGJEkodV90dF81eEdGJEkndV90XzZ4R0YkSSV1Xzd4R0YkSSV1XzZ0R0YkSSd1XzV0X3hHRiRJKHVfNHRfeHhHRiRJKXVfdHR0eHh4R0YkSSh1X3R0XzR4R0YkSSd1X3RfNXhHRiRJJXVfNnhHRiRJJXVfNXRHRiRJJ3VfNHRfeEdGJEkodV90dHR4eEdGJEkodV90dHh4eEdGJEkndV90XzR4R0YkSSV1XzV4R0YkSSV1XzR0R0YkSSd1X3R0dHhHRiRJJ3VfdHR4eEdGJEkndV90eHh4R0YkSSV1XzR4R0YkSSZ1X3R0dEdGJEkmdV90dHhHRiRJJnVfdHh4R0YkSSZ1X3h4eEdGJEkldV90dEdGJEkldV90eEdGJEkldV94eEdGJEkkdV90R0YkSSR1X3hHRiRJJXZfOHRHRiRJJ3ZfN3RfeEdGJEkodl82dF94eEdGJEkpdl81dF94eHhHRiRJKHZfNHRfNHhHRiRJKXZfdHR0XzV4R0YkSSh2X3R0XzZ4R0YkSSd2X3RfN3hHRiRJJXZfOHhHRiRJJXZfN3RHRiRJJ3ZfNnRfeEdGJEkodl81dF94eEdGJEkpdl80dF94eHhHRiRJKXZfdHR0XzR4R0YkSSh2X3R0XzV4R0YkSSd2X3RfNnhHRiRJJXZfN3hHRiRJJXZfNnRHRiRJJ3ZfNXRfeEdGJEkodl80dF94eEdGJEkpdl90dHR4eHhHRiRJKHZfdHRfNHhHRiRJJ3ZfdF81eEdGJEkldl82eEdGJEkldl81dEdGJEkndl80dF94R0YkSSh2X3R0dHh4R0YkSSh2X3R0eHh4R0YkSSd2X3RfNHhHRiRJJXZfNXhHRiRJJXZfNHRHRiRJJ3ZfdHR0eEdGJEkndl90dHh4R0YkSSd2X3R4eHhHRiRJJXZfNHhHRiRJJnZfdHR0R0YkSSZ2X3R0eEdGJEkmdl90eHhHRiRJJnZfeHh4R0YkSSV2X3R0R0YkSSV2X3R4R0YkSSV2X3h4R0YkSSR2X3RHRiRJJHZfeEdGJA==
<Text-field style="Heading 2" layout="Heading 2">Example 1</Text-field> put('pd(U,u_x^2)' = -u_x*pd(U,u)); pd(U,u_x^2); LCQqJkkkdV94RzYiIiIiLUkjcGRHRiU2JEkiVUdGJUkidUdGJUYmISIi pd(U,u_x^3); LCYtSSNwZEc2IjYkSSJVR0YlSSJ1R0YlISIiKiZJJHVfeEdGJSIiIi1GJDYkRicqJkYrRixGKEYsRixGKQ== clear(pds); NyMsJi1JI3BkRzYiNiRJIlVHNiIqJEkkdV94RzYiIiIjIiIiKiZJJHVfeEc2IiIiIi1JI3BkRzYiNiRJIlVHNiJJInVHNiIiIiIiIiI= pd(U,u_x^3); LUkjcGRHNiI2JEkiVUdGJCokSSR1X3hHRiQiIiQ=
<Text-field style="Heading 2" layout="Heading 2">Example 2</Text-field> Resolving Cauchy\342\200\223Riemann conditions put(pd(U,u) = pd(V,v), pd(V,u) = -pd(U,v)); produces the Laplace equation pd(U,u^2); i.e. pd(U,u^2)+pd(U,v^2)=0; clear(pds):
<Text-field style="Heading 2" layout="Heading 2">Cross-derivatives and compatibility conditions</Text-field> Substitutions can contradict one another, generating chaotic and session-dependent results. For instance, the above Cauchy\342\200\223Riemann conditions can be alternatively resolved as put( 'pd(V,v)' = pd(U,u), 'pd(V,u)' = -pd(U,v)); but then the input pd(V,u*v); LUkjcGRHNiI2JEkiVUdGJCokSSJ1R0YkIiIj may produce different results depending on whether it resulted from the \357\254\201rst or the second substitution. cc(); PCMsJi1JI3BkRzYiNiRJIlVHNiIqJEkidkc2IiIiIyEiIi1JI3BkRzYiNiRJIlVHNiIqJEkidUc2IiIiIyEiIg== put('pd(U,u^2)' = -pd(U,v^2)); cc(); PCI=
<Text-field style="Heading 2" layout="Heading 2">Varodering</Text-field> clear(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);
<Text-field style="Heading 2" layout="Heading 2">Trick (two or more unknowns)</Text-field> 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);