DEtools - Maple Programming Help

Home : Support : Online Help : Mathematics : Differential Equations : Rif : DEtools/rifsimp

DEtools

 rifsimp
 simplify overdetermined polynomially nonlinear PDE or ODE systems

 Calling Sequence rifsimp(system, options) rifsimp(system, vars, options)

Parameters

 system - list or set of polynomially nonlinear PDEs or ODEs (may contain inequations) vars - (optional) list of the main dependent variables options - (optional) sequence of options to control the behavior of rifsimp

Description

 • The rifsimp function can be used to simplify or rework overdetermined systems of polynomially nonlinear PDEs or ODEs and inequations to a more useful form -- relative standard form. The rifsimp function does not solve PDE systems, but provides existence and uniqueness information, and can be used as a first step to their solution. As an example, inconsistent systems can be detected by rifsimp.
 • Basically, given an input PDE system, and a list of dependent variables or constants to solve for, rifsimp returns the simplified PDE system along with any existence conditions required for the simplified system to hold.
 • Detailed examples of the use of rifsimp for various systems (along with some explanation of the algorithm) can be found in rifsimp[overview].
 • Other options are sometimes required along with the specification of the system and its solving variables. For common options, please see rifsimp[options], and for more advanced use, please see rifsimp[adv_options].
 • For a description of all possible output configurations, see rifsimp[output].

Examples

1. Overdetermined Systems

As a first example, we have the overdetermined system of two equations in one dependent variable y(x).

 > $\mathrm{sys1}≔\left[x{\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)}^{2}{\left(\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)\right)}^{2}-2x\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)\left(\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)\right)y\left(x\right)\left(\frac{{ⅆ}^{3}}{ⅆ{x}^{3}}y\left(x\right)\right)+x{y\left(x\right)}^{2}{\left(\frac{{ⅆ}^{3}}{ⅆ{x}^{3}}y\left(x\right)\right)}^{2}-y\left(x\right)\left(\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)\right)+{\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)}^{2}=0,-\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)\left(\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)\right)+y\left(x\right)\left(\frac{{ⅆ}^{3}}{ⅆ{x}^{3}}y\left(x\right)\right)+2{y\left(x\right)}^{2}{\left(\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)\right)}^{2}-4y\left(x\right)\left(\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)\right){\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)}^{2}+2{\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)}^{4}=0\right]$
 ${\mathrm{sys1}}{≔}\left[{x}{}{\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right)}^{{2}}{}{\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right)\right)}^{{2}}{-}{2}{}{x}{}\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right)\right){}{y}{}\left({x}\right){}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}{}{y}{}\left({x}\right)\right){+}{x}{}{{y}{}\left({x}\right)}^{{2}}{}{\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}{}{y}{}\left({x}\right)\right)}^{{2}}{-}{y}{}\left({x}\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right)\right){+}{\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right)}^{{2}}{=}{0}{,}{-}\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right)\right){+}{y}{}\left({x}\right){}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}{}{y}{}\left({x}\right)\right){+}{2}{}{{y}{}\left({x}\right)}^{{2}}{}{\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right)\right)}^{{2}}{-}{4}{}{y}{}\left({x}\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right)\right){}{\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right)}^{{2}}{+}{2}{}{\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right)}^{{4}}{=}{0}\right]$ (1)

Now call rifsimp.

 > $\mathrm{with}\left(\mathrm{DEtools}\right):$
 > $\mathrm{ans1}≔\mathrm{rifsimp}\left(\mathrm{sys1}\right)$
 ${\mathrm{ans1}}{≔}{\mathrm{table}}\left(\left[{\mathrm{Solved}}{=}\left[\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right){=}\frac{{\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right)}^{{2}}}{{y}{}\left({x}\right)}\right]{,}{\mathrm{Pivots}}{=}\left[{y}{}\left({x}\right){\ne }{0}\right]{,}{\mathrm{Case}}{=}\left[\left[{y}{}\left({x}\right){\ne }{0}{,}\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}{}{y}{}\left({x}\right)\right]\right]\right]\right)$ (2)

We see that the system has been reduced (for the case y(x) not identically zero) to one much simpler equation that can now be handled with Maple's dsolve.

 > $\mathrm{dsolve}\left({\mathrm{ans1}}_{\mathrm{Solved}},\left\{y\left(x\right)\right\},\mathrm{explicit}\right)$
 $\left\{{y}{}\left({x}\right){=}{{ⅇ}}^{{\mathrm{_C1}}{}{x}}{}{\mathrm{_C2}}\right\}$ (3)

In addition to assisting dsolve in obtaining exact solutions, the simplified form can be used in combination with the initialdata function to obtain the required initial data for the solved form.

 > $\mathrm{id1}≔\mathrm{initialdata}\left({\mathrm{ans1}}_{\mathrm{Solved}}\right)$
 ${\mathrm{id1}}{≔}{\mathrm{table}}\left(\left[{\mathrm{Infinite}}{=}\left[{}\right]{,}{\mathrm{Finite}}{=}\left[{y}{}\left({{x}}_{{0}}\right){=}{\mathrm{_C1}}{,}{\mathrm{D}}{}\left({y}\right){}\left({{x}}_{{0}}\right){=}{\mathrm{_C2}}\right]\right]\right)$ (4)

In this case, it gives us the expected results for a second order ODE, but it can also calculate the required initial data for more complex PDE systems (see initialdata for more information). Numerical methods can now be successfully applied to the reduced system, with initial conditions of the type calculated by initialdata.

 > $\mathrm{dsolve}\left(\mathrm{convert}\left({\mathrm{ans1}}_{\mathrm{Solved}},\mathrm{set}\right)\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}∪\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\left\{y\left(0\right)=1,\mathrm{D}\left(y\right)\left(0\right)=2\right\},\left\{y\left(x\right)\right\},\mathrm{numeric},\mathrm{output}=\mathrm{Array}\left(\left[0,0.25,0.5,0.75,1.0\right]\right)\right)$
 $\left[\begin{array}{c}\left[\begin{array}{ccc}{x}& {y}{}\left({x}\right)& \frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\end{array}\right]\\ \left[\begin{array}{ccc}{0.}& {1.}& {2.}\\ {0.250000000000000}& {1.64872125401162}& {3.29744250802324}\\ {0.500000000000000}& {2.71828135845628}& {5.43656271691257}\\ {0.750000000000000}& {4.48168797609521}& {8.96337595219041}\\ {1.}& {7.38905394806921}& {14.7781078961384}\end{array}\right]\end{array}\right]$ (5)

Finally, we can also obtain a formal power series (Taylor series) for the local solution of this problem (see rtaylor for more information).

 > $\mathrm{tay_ser}≔\mathrm{rtaylor}\left({\mathrm{ans1}}_{\mathrm{Solved}},\mathrm{order}=4\right)$
 ${\mathrm{tay_ser}}{≔}\left[{y}{}\left({x}\right){=}{y}{}\left({{x}}_{{0}}\right){+}{\mathrm{D}}{}\left({y}\right){}\left({{x}}_{{0}}\right){}\left({x}{-}{{x}}_{{0}}\right){+}\frac{{1}}{{2}}{}\frac{{{\mathrm{D}}{}\left({y}\right){}\left({{x}}_{{0}}\right)}^{{2}}{}{\left({x}{-}{{x}}_{{0}}\right)}^{{2}}}{{y}{}\left({{x}}_{{0}}\right)}{+}\frac{{1}}{{6}}{}\frac{{{\mathrm{D}}{}\left({y}\right){}\left({{x}}_{{0}}\right)}^{{3}}{}{\left({x}{-}{{x}}_{{0}}\right)}^{{3}}}{{{y}{}\left({{x}}_{{0}}\right)}^{{2}}}{+}\frac{{1}}{{24}}{}\frac{{{\mathrm{D}}{}\left({y}\right){}\left({{x}}_{{0}}\right)}^{{4}}{}{\left({x}{-}{{x}}_{{0}}\right)}^{{4}}}{{{y}{}\left({{x}}_{{0}}\right)}^{{3}}}\right]$ (6)

The above series is fully determined when given the initial data described earlier:

 > $\mathrm{tay_ser}≔\mathrm{rtaylor}\left({\mathrm{ans1}}_{\mathrm{Solved}},\mathrm{id1},\mathrm{order}=4\right)$
 ${\mathrm{tay_ser}}{≔}\left[{y}{}\left({x}\right){=}{\mathrm{_C1}}{+}{\mathrm{_C2}}{}\left({x}{-}{{x}}_{{0}}\right){+}\frac{{1}}{{2}}{}\frac{{{\mathrm{_C2}}}^{{2}}{}{\left({x}{-}{{x}}_{{0}}\right)}^{{2}}}{{\mathrm{_C1}}}{+}\frac{{1}}{{6}}{}\frac{{{\mathrm{_C2}}}^{{3}}{}{\left({x}{-}{{x}}_{{0}}\right)}^{{3}}}{{{\mathrm{_C1}}}^{{2}}}{+}\frac{{1}}{{24}}{}\frac{{{\mathrm{_C2}}}^{{4}}{}{\left({x}{-}{{x}}_{{0}}\right)}^{{4}}}{{{\mathrm{_C1}}}^{{3}}}\right]$ (7)

2. Inconsistent Systems

rifsimp can often determine whether a system is inconsistent.

 > $\mathrm{sys2}≔\left[{\left(\frac{\partial }{\partial x}u\left(x,y\right)\right)}^{2}+\frac{{\partial }^{2}}{\partial {y}^{2}}u\left(x,y\right)=0,{u\left(x,y\right)}^{3}-\left(\frac{\partial }{\partial x}u\left(x,y\right)\right)+x=0\right]:$
 > $\mathrm{rifsimp}\left(\mathrm{sys2}\right)$
 ${\mathrm{table}}\left(\left[{\mathrm{status}}{=}{"system is inconsistent"}\right]\right)$ (8)

3. Constrained Mechanical Systems

This example shows the use of rifsimp as a preprocessor for a constrained mechanical system (that is, a Differential-Algebraic Equation or DAE system). The method of Lagrange formulates the motion of a bead of nonzero mass m on a frictionless wire of shape Phi(x,y) under gravity as follows:

 > $\mathrm{ConstrainedSys}≔\left[m\ne 0,\mathrm{Φ}\left(x\left(t\right),y\left(t\right)\right)=0,m\left(\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}x\left(t\right)\right)=-\mathrm{λ}\left(t\right)\mathrm{D}[1]\left(\mathrm{Φ}\right)\left(x\left(t\right),y\left(t\right)\right),m\left(\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}y\left(t\right)\right)=-\mathrm{λ}\left(t\right)\mathrm{D}[2]\left(\mathrm{Φ}\right)\left(x\left(t\right),y\left(t\right)\right)-mg\right]$
 ${\mathrm{ConstrainedSys}}{≔}\left[{m}{\ne }{0}{,}{\mathrm{Φ}}{}\left({x}{}\left({t}\right){,}{y}{}\left({t}\right)\right){=}{0}{,}{m}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{x}{}\left({t}\right)\right){=}{-}{\mathrm{λ}}{}\left({t}\right){}{\mathrm{D}}{[}{1}{]}{}\left({\mathrm{Φ}}\right){}\left({x}{}\left({t}\right){,}{y}{}\left({t}\right)\right){,}{m}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{y}{}\left({t}\right)\right){=}{-}{\mathrm{λ}}{}\left({t}\right){}{\mathrm{D}}{[}{2}{]}{}\left({\mathrm{Φ}}\right){}\left({x}{}\left({t}\right){,}{y}{}\left({t}\right)\right){-}{m}{}{g}\right]$ (9)

Note that a mass falling under gravity without air resistance corresponds to:

 > $\mathrm{Φ}≔\left(x,y\right)→0$
 ${\mathrm{Φ}}{≔}\left({x}{,}{y}\right){→}{0}$ (10)

We pick the example of a pendulum, so the bead moves on a circular wire.

 > $\mathrm{Φ}≔\left(x,y\right)→{x}^{2}+{y}^{2}-1$
 ${\mathrm{Φ}}{≔}\left({x}{,}{y}\right){→}{{x}}^{{2}}{+}{{y}}^{{2}}{-}{1}$ (11)
 > $\mathrm{ConstrainedSys}$
 $\left[{m}{\ne }{0}{,}{{x}{}\left({t}\right)}^{{2}}{+}{{y}{}\left({t}\right)}^{{2}}{-}{1}{=}{0}{,}{m}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{x}{}\left({t}\right)\right){=}{-}{2}{}{\mathrm{λ}}{}\left({t}\right){}{x}{}\left({t}\right){,}{m}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{y}{}\left({t}\right)\right){=}{-}{2}{}{\mathrm{λ}}{}\left({t}\right){}{y}{}\left({t}\right){-}{m}{}{g}\right]$ (12)

Such constrained systems present great difficulties for numerical solvers. Of course, we could eliminate the constraint using polar coordinates, but in general this is impossible for mechanical systems with complicated constraints. For example, just replace the constraint with another function, so that the wire has a different shape.

 > $\mathrm{Φ2}≔\left(x,y\right)→{x}^{4}+{y}^{4}-1$
 ${\mathrm{Φ2}}{≔}\left({x}{,}{y}\right){→}{{x}}^{{4}}{+}{{y}}^{{4}}{-}{1}$ (13)

(See Reid et al., 1996, given in the reference and package information page Rif). The pendulum is the classic example of such systems that have been the focus of much recent research due to their importance in applications.

 > $\mathrm{SimpSys}≔\mathrm{rifsimp}\left(\mathrm{ConstrainedSys}\right)$
 ${\mathrm{SimpSys}}{≔}{\mathrm{table}}\left(\left[{\mathrm{Constraint}}{=}\left[{g}{}{m}{}{{y}{}\left({t}\right)}^{{3}}{+}{\left(\frac{{ⅆ}}{{ⅆ}{t}}{}{y}{}\left({t}\right)\right)}^{{2}}{}{m}{+}{2}{}{\mathrm{λ}}{}\left({t}\right){}{{y}{}\left({t}\right)}^{{2}}{-}{g}{}{m}{}{y}{}\left({t}\right){-}{2}{}{\mathrm{λ}}{}\left({t}\right){=}{0}{,}{{x}{}\left({t}\right)}^{{2}}{+}{{y}{}\left({t}\right)}^{{2}}{-}{1}{=}{0}\right]{,}{\mathrm{Solved}}{=}\left[\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{y}{}\left({t}\right){=}{-}\frac{{2}{}{\mathrm{λ}}{}\left({t}\right){}{y}{}\left({t}\right){+}{m}{}{g}}{{m}}{,}\frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{λ}}{}\left({t}\right){=}{-}\frac{{3}}{{2}}{}{g}{}{m}{}\left(\frac{{ⅆ}}{{ⅆ}{t}}{}{y}{}\left({t}\right)\right){,}\frac{{ⅆ}}{{ⅆ}{t}}{}{x}{}\left({t}\right){=}{-}\frac{{y}{}\left({t}\right){}\left(\frac{{ⅆ}}{{ⅆ}{t}}{}{y}{}\left({t}\right)\right)}{{x}{}\left({t}\right)}\right]{,}{\mathrm{Pivots}}{=}\left[{m}{\ne }{0}{,}{x}{}\left({t}\right){\ne }{0}{,}{y}{}\left({t}\right){-}{1}{\ne }{0}{,}{y}{}\left({t}\right){+}{1}{\ne }{0}\right]{,}{\mathrm{Case}}{=}\left[\left[{x}{}\left({t}\right){\ne }{0}{,}\frac{{ⅆ}}{{ⅆ}{t}}{}{x}{}\left({t}\right)\right]{,}\left[\left({y}{}\left({t}\right){-}{1}\right){}\left({y}{}\left({t}\right){+}{1}\right){\ne }{0}{,}\frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{λ}}{}\left({t}\right)\right]\right]\right]\right)$ (14)
 > $\mathrm{initialdata}\left(\mathrm{SimpSys}\right)$
 ${\mathrm{table}}\left(\left[{\mathrm{Constraint}}{=}\left[{{\mathrm{_C3}}}^{{3}}{}{\mathrm{_C5}}{}{\mathrm{_C6}}{+}{2}{}{\mathrm{_C1}}{}{{\mathrm{_C3}}}^{{2}}{-}{\mathrm{_C3}}{}{\mathrm{_C5}}{}{\mathrm{_C6}}{+}{{\mathrm{_C4}}}^{{2}}{}{\mathrm{_C6}}{-}{2}{}{\mathrm{_C1}}{=}{0}{,}{{\mathrm{_C2}}}^{{2}}{+}{{\mathrm{_C3}}}^{{2}}{-}{1}{=}{0}\right]{,}{\mathrm{Pivots}}{=}\left[{\mathrm{_C6}}{\ne }{0}{,}{\mathrm{_C2}}{\ne }{0}{,}{\mathrm{_C3}}{-}{1}{\ne }{0}{,}{\mathrm{_C3}}{+}{1}{\ne }{0}\right]{,}{\mathrm{Infinite}}{=}\left[{}\right]{,}{\mathrm{Finite}}{=}\left[{\mathrm{λ}}{}\left({{t}}_{{0}}\right){=}{\mathrm{_C1}}{,}{x}{}\left({{t}}_{{0}}\right){=}{\mathrm{_C2}}{,}{y}{}\left({{t}}_{{0}}\right){=}{\mathrm{_C3}}{,}{\mathrm{D}}{}\left({y}\right){}\left({{t}}_{{0}}\right){=}{\mathrm{_C4}}{,}{g}{=}{\mathrm{_C5}}{,}{m}{=}{\mathrm{_C6}}\right]\right]\right)$ (15)

Although there is no simplification as in the previous examples, rifsimp has found additional constraints that the initial conditions must satisfy (Constraint and Pivots in the initialdata output), and obtained an equation for the time derivative of lambda(t). Maple's dsolve[numeric] can then be used to obtain a numerical solution, in the manner of the previous example.

4. Lie Symmetry Determination

This example shows the use of rifsimp to assist in determination of Lie point symmetries of an ODE. Given the ODE:

 > $\mathrm{ODE}≔2{x}^{2}y\left(x\right)\left(\frac{{ⅆ}^{2}}{ⅆ{x}^{2}}y\left(x\right)\right)-{x}^{2}\left({\left(\frac{ⅆ}{ⅆx}y\left(x\right)\right)}^{2}+1\right)+{y\left(x\right)}^{2}$
 ${\mathrm{ODE}}{≔}{2}{}{{x}}^{{2}}{}{y}{}\left({x}\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}{}{y}{}\left({x}\right)\right){-}{{x}}^{{2}}{}\left({\left(\frac{{ⅆ}}{{ⅆ}{x}}{}{y}{}\left({x}\right)\right)}^{{2}}{+}{1}\right){+}{{y}{}\left({x}\right)}^{{2}}$ (16)

We can obtain the system of PDE for it's Lie point symmetries via:

 > $\mathrm{sys}≔\left\{\mathrm{coeffs}\left(\mathrm{numer}\left(\mathrm{DEtools}[\mathrm{odepde}]\left(\mathrm{ODE},\left[\mathrm{ξ},\mathrm{η}\right]\left(x,y\right),y\left(x\right)\right)\right),\mathrm{_y1}\right)\right\}$
 ${\mathrm{sys}}{≔}\left\{{-}{2}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{y}}^{{2}}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{{y}}^{{2}}{-}\left(\frac{{\partial }}{{\partial }{y}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{y}{,}{2}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{y}}^{{2}}}{}{\mathrm{η}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{{y}}^{{2}}{-}{4}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{y}{}{\partial }{x}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{{y}}^{{2}}{-}\left(\frac{{\partial }}{{\partial }{y}}{}{\mathrm{η}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{y}{+}{\mathrm{η}}{}\left({x}{,}{y}\right){}{{x}}^{{3}}{,}{4}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{y}{}{\partial }{x}}{}{\mathrm{η}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{{y}}^{{2}}{-}{2}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{x}}^{{2}}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{{y}}^{{2}}{-}{3}{}\left(\frac{{\partial }}{{\partial }{y}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{y}{+}{3}{}\left(\frac{{\partial }}{{\partial }{y}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{x}{}{{y}}^{{3}}{-}{2}{}\left(\frac{{\partial }}{{\partial }{x}}{}{\mathrm{η}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{y}{,}{2}{}\left(\frac{{{\partial }}^{{2}}}{{\partial }{{x}}^{{2}}}{}{\mathrm{η}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{{y}}^{{2}}{-}{2}{}\left(\frac{{\partial }}{{\partial }{x}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{y}{+}{2}{}\left(\frac{{\partial }}{{\partial }{x}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right)\right){}{x}{}{{y}}^{{3}}{+}\left(\frac{{\partial }}{{\partial }{y}}{}{\mathrm{η}}{}\left({x}{,}{y}\right)\right){}{{x}}^{{3}}{}{y}{-}\left(\frac{{\partial }}{{\partial }{y}}{}{\mathrm{η}}{}\left({x}{,}{y}\right)\right){}{x}{}{{y}}^{{3}}{+}{\mathrm{η}}{}\left({x}{,}{y}\right){}{{x}}^{{3}}{+}{\mathrm{η}}{}\left({x}{,}{y}\right){}{x}{}{{y}}^{{2}}{-}{2}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right){}{{y}}^{{3}}\right\}$ (17)

We can then use rifsimp to greatly simplify this system:

 > $\mathrm{rifsys}≔\mathrm{rifsimp}\left(\mathrm{sys}\right)$
 ${\mathrm{rifsys}}{≔}{\mathrm{table}}\left(\left[{\mathrm{Solved}}{=}\left[\frac{{{\partial }}^{{2}}}{{\partial }{{x}}^{{2}}}{}{\mathrm{η}}{}\left({x}{,}{y}\right){=}\frac{{\mathrm{ξ}}{}\left({x}{,}{y}\right){}{y}{-}{\mathrm{η}}{}\left({x}{,}{y}\right){}{x}}{{{x}}^{{3}}}{,}\frac{{\partial }}{{\partial }{x}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right){=}\frac{{\mathrm{η}}{}\left({x}{,}{y}\right)}{{y}}{,}\frac{{\partial }}{{\partial }{y}}{}{\mathrm{η}}{}\left({x}{,}{y}\right){=}\frac{{\mathrm{η}}{}\left({x}{,}{y}\right)}{{y}}{,}\frac{{\partial }}{{\partial }{y}}{}{\mathrm{ξ}}{}\left({x}{,}{y}\right){=}{0}\right]\right]\right)$ (18)

From this point, pdsolve can be used to obtain the infinitesimals for the Lie point symmetries via:

 > $\mathrm{pdsolve}\left({\mathrm{rifsys}}_{'\mathrm{Solved}'}\right)$
 $\left\{{\mathrm{η}}{}\left({x}{,}{y}\right){=}\left({\mathrm{_C3}}{}{{\mathrm{ln}}{}\left({x}\right)}^{{2}}{+}{\mathrm{_C2}}{}{\mathrm{ln}}{}\left({x}\right){+}{\mathrm{_C1}}\right){}{y}{,}{\mathrm{ξ}}{}\left({x}{,}{y}\right){=}\left({\mathrm{_C3}}{}{{\mathrm{ln}}{}\left({x}\right)}^{{2}}{+}\left({\mathrm{_C2}}{-}{2}{}{\mathrm{_C3}}\right){}{\mathrm{ln}}{}\left({x}\right){+}{\mathrm{_C1}}{-}{\mathrm{_C2}}{+}{2}{}{\mathrm{_C3}}\right){}{x}\right\}$ (19)