PDEtools - Maple Programming Help

Home : Support : Online Help : Mathematics : Differential Equations : PDEtools : Symmetries : PDEtools/ConsistencyTest

PDEtools

 ConsistencyTest
 tests whether a given system of equations is or not consistent.

 Calling Sequence ConsistencyTest(S, PDESYS, DepVars)

Parameters

 PDESYS - a PDE or a set or list of them; it can include ODEs and not differential equations or inequations DepVars - optional - may be required; a function or a list of them indicating the dependent variables of the problem outputthesystem - optional - to return a system equivalent to PDESYS simplified with respect to its integrability conditions no_Fn - optional, related to outputthesystem, to specify that all auxiliary functions $\mathrm{_Fn}$ (n integer) introduced in the process are to be eliminated before outputting

Description

 • ConsistencyTest tests whether a system of equations, possibly containing inequations and involving or not differential equations, is or not consistent. A system is inconsistent when one of its equations contradicts one other one. In order to be solvable, a system must be consistent in the first place.
 • If DepVars is not specified, ConsistencyTest will consider all the differentiated unknown functions in PDESYS as unknown of the problems.
 • By default, ConsistencyTest returns true or false. If the system is consistent when the option argument outputthesystem is specified, ConsistencyTest returns three lists with equations, inequations and function replacements to represent PDESYS in differential polynomial form, where the equations and inequations are equivalent to PDESYS simplified with respect to its integrability conditions. These three lists  returned can also be computed using dpolyform, but the number of operations performed to compute it is smaller when using ConsistencyTest. If PDESYS happens to be inconsistent and outputthesystem is specified ConsistencyTest returns NULL.
 • To obtain output without auxiliary functions when outputthesystem is specified use in addition the optional argument no_Fn.

Examples

 > $\mathrm{with}\left(\mathrm{PDEtools}\right):$
 > $\mathrm{declare}\left(f\left(x,y\right),g\left(x,y\right)\right)$
 ${f}{}\left({x}{,}{y}\right){}{\mathrm{will now be displayed as}}{}{f}$
 ${g}{}\left({x}{,}{y}\right){}{\mathrm{will now be displayed as}}{}{g}$ (1)
 > ${\mathrm{sys}}_{1}≔\left[\frac{{\partial }^{2}}{\partial {y}^{2}}f\left(x,y\right)=0,-6\left(\frac{\partial }{\partial y}f\left(x,y\right)\right)y+\frac{{\partial }^{2}}{\partial {y}^{2}}g\left(x,y\right)-2\left(\frac{{\partial }^{2}}{\partial y\partial x}f\left(x,y\right)\right)=0,-9\left(\frac{\partial }{\partial y}f\left(x,y\right)\right)a{y}^{2}-12\left(\frac{\partial }{\partial y}f\left(x,y\right)\right){a}^{2}y-3\left(\frac{\partial }{\partial y}f\left(x,y\right)\right)b-3\left(\frac{\partial }{\partial x}f\left(x,y\right)\right)y+2\left(\frac{{\partial }^{2}}{\partial y\partial x}g\left(x,y\right)\right)-\left(\frac{{\partial }^{2}}{\partial {x}^{2}}f\left(x,y\right)\right)-3g\left(x,y\right)=0,-6g\left(x,y\right)ay-6\left(\frac{\partial }{\partial x}f\left(x,y\right)\right)a{y}^{2}-8\left(\frac{\partial }{\partial x}f\left(x,y\right)\right){a}^{2}y+3\left(\frac{\partial }{\partial y}g\left(x,y\right)\right)a{y}^{2}+4\left(\frac{\partial }{\partial y}g\left(x,y\right)\right){a}^{2}y-4g\left(x,y\right){a}^{2}-2\left(\frac{\partial }{\partial x}f\left(x,y\right)\right)b+\left(\frac{\partial }{\partial y}g\left(x,y\right)\right)b-3y\left(\frac{\partial }{\partial x}g\left(x,y\right)\right)+\frac{{\partial }^{2}}{\partial {x}^{2}}g\left(x,y\right)=0\right]$
 ${{\mathrm{sys}}}_{{1}}{≔}\left[{{f}}_{{y}{,}{y}}{=}{0}{,}{-}{6}{}{{f}}_{{y}}{}{y}{-}{2}{}{{f}}_{{x}{,}{y}}{+}{{g}}_{{y}{,}{y}}{=}{0}{,}{-}{12}{}{{f}}_{{y}}{}{{a}}^{{2}}{}{y}{-}{9}{}{{f}}_{{y}}{}{a}{}{{y}}^{{2}}{-}{3}{}{{f}}_{{y}}{}{b}{-}{3}{}{{f}}_{{x}}{}{y}{-}{3}{}{g}{-}{{f}}_{{x}{,}{x}}{+}{2}{}{{g}}_{{x}{,}{y}}{=}{0}{,}{-}{8}{}{{f}}_{{x}}{}{{a}}^{{2}}{}{y}{+}{4}{}{{g}}_{{y}}{}{{a}}^{{2}}{}{y}{-}{6}{}{{f}}_{{x}}{}{a}{}{{y}}^{{2}}{+}{3}{}{{g}}_{{y}}{}{a}{}{{y}}^{{2}}{-}{4}{}{g}{}{{a}}^{{2}}{-}{6}{}{g}{}{a}{}{y}{-}{2}{}{{f}}_{{x}}{}{b}{+}{{g}}_{{y}}{}{b}{-}{3}{}{y}{}{{g}}_{{x}}{+}{{g}}_{{x}{,}{x}}{=}{0}\right]$ (2)
 > $\mathrm{ConsistencyTest}\left({\mathrm{sys}}_{1}\right)$
 ${\mathrm{true}}$ (3)

Let's see the actual form of the consistent system sys after its integrability conditions are taken into account

 > $\mathrm{ConsistencyTest}\left({\mathrm{sys}}_{1},\mathrm{outputthesystem}\right)$
 $\left[{g}{=}{0}{,}{{f}}_{{x}}{=}{0}{,}{{f}}_{{y}}{=}{0}\right]{,}\left[\right]{,}\left[\right]$ (4)

Construct an inconsistent system departing from sys, by multiplying any two of its equations (both are equal to zero) and equating the result to 1

 > ${\mathrm{sys}}_{2}≔\left[\mathrm{op}\left({\mathrm{sys}}_{1}\right),\mathrm{expand}\left(\mathrm{lhs}\left({{\mathrm{sys}}_{1}}_{1}\right)\mathrm{lhs}\left({{\mathrm{sys}}_{1}}_{2}\right)\right)=1\right]$
 ${{\mathrm{sys}}}_{{2}}{≔}\left[{{f}}_{{y}{,}{y}}{=}{0}{,}{-}{6}{}{{f}}_{{y}}{}{y}{-}{2}{}{{f}}_{{x}{,}{y}}{+}{{g}}_{{y}{,}{y}}{=}{0}{,}{-}{12}{}{{f}}_{{y}}{}{{a}}^{{2}}{}{y}{-}{9}{}{{f}}_{{y}}{}{a}{}{{y}}^{{2}}{-}{3}{}{{f}}_{{y}}{}{b}{-}{3}{}{{f}}_{{x}}{}{y}{-}{3}{}{g}{-}{{f}}_{{x}{,}{x}}{+}{2}{}{{g}}_{{x}{,}{y}}{=}{0}{,}{-}{8}{}{{f}}_{{x}}{}{{a}}^{{2}}{}{y}{+}{4}{}{{g}}_{{y}}{}{{a}}^{{2}}{}{y}{-}{6}{}{{f}}_{{x}}{}{a}{}{{y}}^{{2}}{+}{3}{}{{g}}_{{y}}{}{a}{}{{y}}^{{2}}{-}{4}{}{g}{}{{a}}^{{2}}{-}{6}{}{g}{}{a}{}{y}{-}{2}{}{{f}}_{{x}}{}{b}{+}{{g}}_{{y}}{}{b}{-}{3}{}{y}{}{{g}}_{{x}}{+}{{g}}_{{x}{,}{x}}{=}{0}{,}{-}{6}{}{{f}}_{{y}{,}{y}}{}{{f}}_{{y}}{}{y}{-}{2}{}{{f}}_{{y}{,}{y}}{}{{f}}_{{x}{,}{y}}{+}{{f}}_{{y}{,}{y}}{}{{g}}_{{y}{,}{y}}{=}{1}\right]$ (5)
 > $\mathrm{ConsistencyTest}\left({\mathrm{sys}}_{2}\right)$
 ${\mathrm{false}}$ (6)

When the system is inconsistent and the option outputthesystem` is received, ConsistencyTest returns NULL and a related warning message is displayed

 > $\mathrm{ConsistencyTest}\left({\mathrm{sys}}_{2},\mathrm{outputthesystem}\right)$
 Warning: System is inconsistent

An example involving mathematical functions (exp)

 > ${\mathrm{sys}}_{3}≔\left[\frac{{\partial }^{2}}{\partial {y}^{2}}f\left(x,y\right)=0,-6\left(\frac{\partial }{\partial y}f\left(x,y\right)\right)y+\frac{{\partial }^{2}}{\partial {y}^{2}}g\left(x,y\right)-2\left(\frac{{\partial }^{2}}{\partial y\partial x}f\left(x,y\right)\right)=0,-9\left(\frac{\partial }{\partial y}f\left(x,y\right)\right)a{y}^{2}-12\left(\frac{\partial }{\partial y}f\left(x,y\right)\right){a}^{2}y-3\left(\frac{\partial }{\partial y}f\left(x,y\right)\right)b-3\left(\frac{\partial }{\partial x}f\left(x,y\right)\right)y+2\left(\frac{{\partial }^{2}}{\partial y\partial x}g\left(x,y\right)\right)-\left(\frac{{\partial }^{2}}{\partial {x}^{2}}f\left(x,y\right)\right)-3{ⅇ}^{x}-3g\left(x,y\right)=0,-6\left({ⅇ}^{x}+g\left(x,y\right)\right)ay-6\left(\frac{\partial }{\partial x}f\left(x,y\right)\right)a{y}^{2}-8\left(\frac{\partial }{\partial x}f\left(x,y\right)\right){a}^{2}y+3\left(\frac{\partial }{\partial y}g\left(x,y\right)\right)a{y}^{2}+4\left(\frac{\partial }{\partial y}g\left(x,y\right)\right){a}^{2}y-4\left({ⅇ}^{x}+g\left(x,y\right)\right){a}^{2}-2\left(\frac{\partial }{\partial x}f\left(x,y\right)\right)b+\left(\frac{\partial }{\partial y}g\left(x,y\right)\right)b-3y\left({ⅇ}^{x}+\frac{\partial }{\partial x}g\left(x,y\right)\right)+{ⅇ}^{x}+\frac{{\partial }^{2}}{\partial {x}^{2}}g\left(x,y\right)=0\right]$
 ${{\mathrm{sys}}}_{{3}}{≔}\left[{{f}}_{{y}{,}{y}}{=}{0}{,}{-}{6}{}{{f}}_{{y}}{}{y}{-}{2}{}{{f}}_{{x}{,}{y}}{+}{{g}}_{{y}{,}{y}}{=}{0}{,}{-}{9}{}{{f}}_{{y}}{}{a}{}{{y}}^{{2}}{-}{12}{}{{f}}_{{y}}{}{{a}}^{{2}}{}{y}{-}{3}{}{{f}}_{{y}}{}{b}{-}{3}{}{{f}}_{{x}}{}{y}{+}{2}{}{{g}}_{{x}{,}{y}}{-}{{f}}_{{x}{,}{x}}{-}{3}{}{{ⅇ}}^{{x}}{-}{3}{}{g}{=}{0}{,}{-}{6}{}\left({{ⅇ}}^{{x}}{+}{g}\right){}{a}{}{y}{-}{6}{}{{f}}_{{x}}{}{a}{}{{y}}^{{2}}{-}{8}{}{{f}}_{{x}}{}{{a}}^{{2}}{}{y}{+}{3}{}{{g}}_{{y}}{}{a}{}{{y}}^{{2}}{+}{4}{}{{g}}_{{y}}{}{{a}}^{{2}}{}{y}{-}{4}{}\left({{ⅇ}}^{{x}}{+}{g}\right){}{{a}}^{{2}}{-}{2}{}{{f}}_{{x}}{}{b}{+}{{g}}_{{y}}{}{b}{-}{3}{}{y}{}\left({{ⅇ}}^{{x}}{+}{{g}}_{{x}}\right){+}{{ⅇ}}^{{x}}{+}{{g}}_{{x}{,}{x}}{=}{0}\right]$ (7)
 > $\mathrm{ConsistencyTest}\left({\mathrm{sys}}_{3}\right)$
 ${\mathrm{true}}$ (8)

The output with three lists, typical of dpolyform, showing the consistent form of sys, involves an auxiliary function _F1 to represent in differential polynomial form the system originally containing exp

 > $\mathrm{ConsistencyTest}\left({\mathrm{sys}}_{3},\mathrm{outputthesystem}\right)$
 $\left[{g}{=}{-}{\mathrm{_F1}}{}\left({x}{,}{y}\right){,}{{f}}_{{x}}{=}{0}{,}{{f}}_{{y}}{=}{0}\right]{,}\left[\right]{,}\left[{\mathrm{_F1}}{}\left({x}{,}{y}\right){=}{{ⅇ}}^{{x}}\right]$ (9)

The output typical of casesplit for the same system does not contain auxiliary functions; to obtain it use the optional argument no_Fn

 > $\mathrm{ConsistencyTest}\left({\mathrm{sys}}_{3},\mathrm{outputthesystem},\mathrm{no_Fn}\right)$
 $\left[{g}{=}{-}{{ⅇ}}^{{x}}{,}{{f}}_{{x}}{=}{0}{,}{{f}}_{{y}}{=}{0}\right]\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{&where}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\left[\right]$ (10)

To compute a form of the system which does not involve auxiliary functions _Fn and also is entirely differential polynomial, so it does not involve mathematical functions, use dpolyform

 > $\mathrm{dpolyform}\left({\mathrm{sys}}_{3},\mathrm{no_Fn}\right)$
 $\left[{{f}}_{{x}}{=}{0}{,}{{f}}_{{y}}{=}{0}{,}{{g}}_{{x}}{=}{g}{,}{{g}}_{{y}}{=}{0}\right]\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{&where}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\left[{g}{\ne }{0}\right]$ (11)

Note however that, depending on the example, the elimination of the auxiliary functions as in the output above may be an expensive computational process.