The Linear Homogeneous PDE commands in the LieAlgebrasOfVectorFields package provide a general purpose toolkit for dealing with linear homogeneous partial differential equations (LHPDE). To be more specific, with these commands you can:

reduce a LHPDEs system using the rifsimp command and with ranking recorded with the system.

find properties of a LHPDEs system such as solution dimension, parametric derivatives, and whether a total degree ranking was used.

explore relationships between solutions of two LHPDEs systems including whether one solution space is a subspace of the other, finding intersection of solution spaces of two LHPDEs, and finding vector space sum of solution spaces of two LHPDEs.

access to utility commands for manipulating LHPDEs.

manipulate as a partial differential operator form.

manipulate a LHPDEs system that has partial dependencies (i.e. the dependent variables may depend on only some of the independent variables).

manipulate a LHPDEs system that is in partially or fully integrated form (i.e. a system with constants of integration involved).

The linear homogeneous PDE functionality have been written focusing on the key LHPDE object, and almost all methods are within the LHPDE object. There only five exports relating to Linear Homogeneous PDE in the LieAlgebrasOfVectorFields package, four of them are for constructing a Maple object, and one is a subpackage containing utility commands.

The existing Maple builtins normal, expand, simplify, indets, has, hastype, type and convert are extended to work for LHPDE and LHPDO objects.

Some of the internal routines for Linear Homogeneous PDE are also available for use as LieAlgebrasOfVectorFields:-LHLibrary:-Routine. For more details, see overview of the LHLibrary subpackage.

To use this functionality, first a LHPDE object must be constructed (e.g. via LieAlgebrasOfVectorFields[LHPDE] command). Then once a LHPDE object is constructed successfully, a variety of methods will be available. For more detail, see Overview of the LHPDE object.

There are five related exports, namely:

The following methods are available for a LHPDE object.

The following methods are available for a LHPDO object.

A LHPDO object can also act as a differential operator. See the LHPDO Object as Operator for more detail.

The following Maple builtins are extended for allowing both LHPDE and LHPDO objects: type, expand, has, hastype, indets, normal, simplify, convert. See LHPDE Object Overloaded Builtins and LHPDO Object Overloaded Builtins for more detail.

The following methods are available for a IDBasis object.

$\mathrm{with}\left(\mathrm{LieAlgebrasOfVectorFields}\right)$

$\left[\mathrm{Differential}\,\mathrm{DisplayStructure}\,\mathrm{Distribution}\,\mathrm{EliminationLAVF}\,\mathrm{EliminationSystem}\,\mathrm{IDBasis}\,\mathrm{IsLinearizable}\,\mathrm{LAVF}\,\mathrm{LHLibrary}\,\mathrm{LHPDE}\,\mathrm{LHPDO}\,\mathrm{MapDE}\,\mathrm{OneForm}\,\mathrm{SymmetryLAVF}\,\mathrm{VFPDO}\,\mathrm{VectorField}\right]$

$\mathrm{detSys}≔\left[\mathrm{diff}\left(\mathrm{\xi }\left(x\,y\right)\,y\,y\right)=0\,\mathrm{diff}\left(\mathrm{\eta }\left(x\,y\right)\,x\right)+\mathrm{diff}\left(\mathrm{\xi }\left(x\,y\right)\,y\right)=0\,\mathrm{diff}\left(\mathrm{\eta }\left(x\,y\right)\,y\right)=0\,\mathrm{diff}\left(\mathrm{\xi }\left(x\,y\right)\,x\right)=0\right]$

$\mathrm{detSys}≔\left[\frac{{\partial }^2}{\partial y^2}\mathrm{\xi }\left(x\,y\right)=0\,\frac{\partial }{\partial x}\mathrm{\eta }\left(x\,y\right)+\frac{\partial }{\partial y}\mathrm{\xi }\left(x\,y\right)=0\,\frac{\partial }{\partial y}\mathrm{\eta }\left(x\,y\right)=0\,\frac{\partial }{\partial x}\mathrm{\xi }\left(x\,y\right)=0\right]$

$\mathrm{E2}≔\mathrm{LHPDE}\left(\mathrm{detSys}\,\mathrm{indep}=\left[x\,y\right]\,\mathrm{dep}=\left[\mathrm{\xi }\,\mathrm{\eta }\right]\right)$

$\mathrm{E2}≔\left[\frac{{\partial }^2}{\partial y^2}\mathrm{\xi }\left(x\,y\right)=0\,\frac{\partial }{\partial x}\mathrm{\eta }\left(x\,y\right)+\frac{\partial }{\partial y}\mathrm{\xi }\left(x\,y\right)=0\,\frac{\partial }{\partial y}\mathrm{\eta }\left(x\,y\right)=0\,\frac{\partial }{\partial x}\mathrm{\xi }\left(x\,y\right)=0\right],\mathrm{indep}=\left[x\,y\right],\mathrm{dep}=\left[\mathrm{\xi }\left(x\,y\right)\,\mathrm{\eta }\left(x\,y\right)\right]$

$\mathrm{exports}\left(\mathrm{E2}\,'\mathrm{static}'\right)$

$\mathrm{type},\mathrm{convert},\mathrm{normal},\mathrm{expand},\mathrm{simplify},\mathrm{indets},\mathrm{has},\mathrm{hastype},\mathrm{GetIndependents},\mathrm{GetDependents},\mathrm{GetSystem},\mathrm{GetRanking},\mathrm{SetIDBasis},\mathrm{GetIDBasis},\mathrm{Copy},\mathrm{Augment},\mathrm{SolutionDimension},\mathrm{IsFiniteType},\mathrm{IsTrivial},\mathrm{ParametricDerivatives},\mathrm{OrderOfInvolution},\mathrm{IsRifReduced},\mathrm{IsTotalDegreeRanking},\mathrm{AreSameSpace},\mathrm{AreSame},\mathrm{RifReduce},\mathrm{ReducedForm},\mathrm{AdjustDependencies},\mathrm{Intersection},\mathrm{VectorSpaceSum},\mathrm{EliminationSystem},\mathrm{IsSubspace},\mathrm{DChange},\mathrm{dchange},\mathrm{LHSolve},\mathrm{initialdata},\mathrm{InitialData},\mathrm{ModulePrint},\mathrm{ModuleCopy},\mathrm{ModuleApply}$

$\mathrm{type}\left(\mathrm{E2}\,'\mathrm{LHPDE}'\right)$

$\mathrm{true}$

$\mathrm{\Delta }≔\mathrm{convert}\left(\mathrm{E2}\,'\mathrm{LHPDO}'\right)$

$\mathrm{\Delta }≔\left(\mathrm{\ξ}\,\mathrm{\η}\right)\→\left[\frac{\partial }{\partial y}\left(\frac{\partial }{\partial y}\mathrm{\ξ}\right)\,\frac{\partial }{\partial x}\mathrm{\η}\+\frac{\partial }{\partial y}\mathrm{\ξ}\,\frac{\partial }{\partial y}\mathrm{\η}\,\frac{\partial }{\partial x}\mathrm{\ξ}\right]$

$\mathrm{\Delta }\left(\left[-\left(y-\mathrm{y0}\right)\,x-\mathrm{x0}\right]\right)$

$\left[0\,0\,0\,0\right]$

$\mathrm{E2p}≔\mathrm{RifReduce}\left(\mathrm{E2}\right)$

$\mathrm{E2p}≔\left[\frac{{\partial }^2}{\partial y^2}\mathrm{\xi }\left(x\,y\right)=0\,\frac{\partial }{\partial x}\mathrm{\xi }\left(x\,y\right)=0\,\frac{\partial }{\partial x}\mathrm{\eta }\left(x\,y\right)=-\frac{\partial }{\partial y}\mathrm{\xi }\left(x\,y\right)\,\frac{\partial }{\partial y}\mathrm{\eta }\left(x\,y\right)=0\right],\mathrm{indep}=\left[x\,y\right],\mathrm{dep}=\left[\mathrm{\xi }\left(x\,y\right)\,\mathrm{\eta }\left(x\,y\right)\right]$

$\mathrm{SolutionDimension}\left(\mathrm{E2p}\right)$

$3$

$\mathrm{LHLibrary}:-\mathrm{IsHomogeneous}\left(\mathrm{GetSystem}\left(\mathrm{E2p}\right)\,\left[\mathrm{\xi }\left(x\,y\right)\,\mathrm{\eta }\left(x\,y\right)\right]\right)$

$\mathrm{true}$

$\mathrm{AdjustDependencies}\left(\mathrm{E2p}\,\mathrm{dep}=''least''\right)$

$\left[\frac{{\ⅆ}^2}{\ⅆy^2}\mathrm{\xi }\left(y\right)=0\,\frac{\ⅆ}{\ⅆx}\mathrm{\eta }\left(x\right)=-\frac{\ⅆ}{\ⅆy}\mathrm{\xi }\left(y\right)\right],\mathrm{indep}=\left[x\,y\right],\mathrm{dep}=\left[\mathrm{\xi }\left(y\right)\,\mathrm{\eta }\left(x\right)\right]$

$\mathrm{LHSolve}\left(\mathrm{E2p}\right)$

$\left[\mathrm{\xi }\left(x\,y\right)=-\mathrm{c\_\_1}y+\mathrm{c\_\_3}\,\mathrm{\eta }\left(x\,y\right)=\mathrm{c\_\_1}x+\mathrm{c\_\_2}\right]$