DEtools - Maple Programming Help

Home : Support : Online Help : Mathematics : Differential Equations : Lie Symmetry Method : Commands for ODEs : DEtools/intfactor

DEtools

 intfactor
 look for integrating factors for a given ODE

 Calling Sequence intfactor(ODE, y(x), _mu = int_factor_form, try_hard=true)

Parameters

 ODE - ordinary differential equation (any order) y(x) - (optional) dependent variable; required when the ODE contains more than one function being differentiated try_hard = true - (optional) to search for integrating factors depending on n-1 variables, where n is the ODE order _mu = int_factor_form - (optional) to search for integrating factors having the indicated form

Description

 • For a given nth order ODE

$\mathrm{ODE}≔N\left(x,y,y\text{'},\mathrm{...},{y}^{\left(n-1\right)}\right){y}^{\left(n\right)}+M\left(x,y,y\text{'},\mathrm{...},{y}^{\left(n-1\right)}\right)=0$

 intfactor looks for integrating factors -- say mu -- such that mu*ODE is an exact ODE, that is, a total derivative with respect to $x$, of some function $R$ of $x$, $y$, and its derivatives up to order $n-1$. As a rule, intfactor first isolates the highest derivative, and then looks for the mus. When successful, intfactor returns mu or a sequence of them; otherwise it returns NULL. To see the PDE satisfied by mu in the general case, use odepde and gensys.
 • The knowledge of an integrating factor -- say mu -- is sufficient to reduce the order of a given ODE by one; in the case of first order ODEs this reduction of order leads directly to their solution. The reduced ODEs obtained from the knowledge of mu are also called first integrals and can be computed by using the firint command (see examples below).
 • By default intfactor looks for "as many" integrating factors as the differential order of the given ODE, and it returns an answer as soon as these integrating factors are found, or when all its schemes have already been tried.
 • In the case of first order ODEs, an integrating factor depending on only the dependent or the independent variable -- say $x$ and $y$ -- is calculated as explained in standard textbooks. Integrating factors depending on both $x$ and $y$ are built by first determining pairs of infinitesimals of underlying Lie groups (see symgen) and using this information to build the mus.
 • In the case of second order ODEs, after isolating $y\text{'}\text{'}$, intfactor is able to find integrating factors of the form $\mu \left(x,y'\right)$, $\mu \left(y,y'\right)$, and $\mathrm{\mu }\left(x,y\right)$, when they exist. In the $\mathrm{\mu }\left(x,y\right)$ case, it is sometimes necessary to solve an auxiliary second order linear ODE and the success of the command relies on the ability of dsolve in solving that ODE.
 • In the case of third and higher nth order ODEs, intfactor searches for integrating factors depending on any two of $\left\{x,y,y',y\text{'}\text{'},\mathrm{...}y\left(n-1\right)\right\}$, by formulating and tackling a linear PDE system for mu (see gensys). When integrating factors with such a dependence exist, the success of intfactor in computing them depends on the success of the dsolve and pdsolve commands in solving the differential equations of that PDE system.
 • Given an nth order ODE and calling intfactor with the optional argument try_hard = true, the routines search for mu depending on any n-1 of $\left\{x,y,y',y\text{'}\text{'},\mathrm{...}y\left(n-1\right)\right\}$ by tackling the related linear PDE system for mu as explained in the previous paragraph. Although this "hard" approach widens the chances of success significantly, the number of different mu problems being tackled is equal to the differential order and each of them involves more than two variables, so this approach could be much slower than the default one, especially for ODEs of order 4 or higher.
 • Apart from the computable integrating factor types described, in the case of second and higher order ODEs, for some particular ODE families, intfactor is able to compute mus depending on all the variables, that is, of the form $\mathrm{mu}\left(x,y,y',\mathrm{...}y\left(n-1\right)\right)$, provided that mu is polynomial in the $y\left(n-1\right)$ derivative. In the particular case of third order ODEs intfactor can in addition compute - when they exist - integrating factors of the form $\mathrm{mu}\left(y\text{'}\text{'}\right)$ without restrictions to the functional form of mu.
 • It is also possible to request to intfactor the search of integrating factors of some predefined form by using the option '_mu = int_factor_form'. In this approach, int_factor_form, could be either an algebraic expression or a list of them, containing unknowns in the form of functions of $x,y,y',\mathrm{...}$ and / or any symbols whose values are to be determined. When int_factor_form contains no unknowns, only a direct test of it using mutest is performed to validate or discard the suggestion. When the indicated form contains unknowns, intfactor sets up a related algebraic/ODE/PDE system of equations for these unknowns and attempt to solve it using Maple's solve, dsolve and pdsolve commands. This option is particularly useful when the intfactor algorithms fail in finding any integrating factor, but you have some idea of the functional form of such factors (see example at the end of this page).
 • When, for a given high order ODE, many integrating factors are returned, the sequence of mus can be used to reduce the given ODE by many orders at once. Note however that sometimes not all of the integrating factors lead to "different" reductions of order. The reason is that if mu leads to a reduced ODE -- say R1 -- then $f\left(\mathrm{R1}\right)\mathrm{\mu }$, where f is an arbitrary function, is also an integrating factor and leads to the same reduced ODE. Hence, different subroutines may find different integrating factors which are actually leading to just one reduction of order. To determine if two integrating factors lead to different reductions of order, you can first calculate the related first integrals using firint - see examples.
 • From a given ODE and an integrating factor for it, you can compute the integrating factor of another ODE, obtained from the first one by performing any arbitrary change of variables. To compute this integrating factor for the changed ODE, use muchange. To test whether an expression is an integrating factor for a given ODE, use the mutest command.
 • The intfactor command is in some sense complementary to the redode command in that the output of one of them is the input for the other one and vice versa.
 • This function is part of the DEtools package, and so it can be used in the form intfactor(..) only after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[intfactor](..).

Examples

 > $\mathrm{with}\left(\mathrm{DEtools}\right):$
 > $\mathrm{PDEtools}\left[\mathrm{declare}\right]\left(y\left(x\right),\mathrm{prime}=x\right)$
 ${y}{}\left({x}\right){}{\mathrm{will now be displayed as}}{}{y}$
 ${\mathrm{derivatives with respect to}}{}{x}{}{\mathrm{of functions of one variable will now be displayed with \text{'}}}$ (1)

Here is an ODE depending on an arbitrary function $F$.

 > $\mathrm{ODE}≔\mathrm{diff}\left(y\left(x\right),x\right)=\frac{2a}{{x}^{2}\left(-y\left(x\right)+2F\left(\frac{x{y\left(x\right)}^{2}-4a}{x}\right)a\right)}$
 ${\mathrm{ODE}}{≔}{\mathrm{y\text{'}}}{=}\frac{{2}{}{a}}{{{x}}^{{2}}{}\left({-}{y}{+}{2}{}{F}{}\left(\frac{{x}{}{{y}}^{{2}}{-}{4}{}{a}}{{x}}\right){}{a}\right)}$ (2)

A pair of symmetry infinitesimals for this ODE is given by (see symgen):

 > $\mathrm{infinitesimals}≔\mathrm{symgen}\left(\mathrm{ODE}\right)$
 ${\mathrm{infinitesimals}}{≔}\left[{\mathrm{_ξ}}{=}{y}{}{{x}}^{{2}}{,}{\mathrm{_η}}{=}{-}{2}{}{a}\right]$ (3)

From these infinitesimals, a related integrating factor for ODE is given by:

 > $\mathrm{\mu }≔\mathrm{intfactor}\left(\mathrm{ODE}\right)$
 ${\mathrm{\mu }}{≔}\frac{{-}{y}{+}{2}{}{F}{}\left(\frac{{x}{}{{y}}^{{2}}{-}{4}{}{a}}{{x}}\right){}{a}}{{F}{}\left(\frac{{x}{}{{y}}^{{2}}{-}{4}{}{a}}{{x}}\right)}$ (4)

Note that when the ODE is given as an equation $A=B$, it is converted to the form $A-B=0$ before calculating mu.

Let's see the change on the ODE type produced by mu.

 > $\mathrm{odeadvisor}\left(\mathrm{ODE}\right)$
 $\left[\left[{\mathrm{_1st_order}}{,}{\mathrm{_with_symmetry_\left[F\left(x\right),G\left(y\right)\right]}}\right]\right]$ (5)
 > $\mathrm{odeadvisor}\left(\mathrm{\mu }\mathrm{ODE}\right)$
 $\left[{\mathrm{_exact}}{,}\left[{\mathrm{_1st_order}}{,}{\mathrm{_with_symmetry_\left[F\left(x\right),G\left(y\right)\right]}}\right]\right]$ (6)

Here is a second order example with two integrating factors. Its solution is obtained using them to reduce the ODE order by two at once.

 > $\mathrm{ODE2}≔\mathrm{diff}\left(y\left(x\right),x,x\right)-\frac{\frac{1}{2}y\left(x\right)\left(4+y\left(x\right)x\right)}{{x}^{2}}-y\left(x\right)\mathrm{diff}\left(y\left(x\right),x\right)-\frac{1}{2}x{\mathrm{diff}\left(y\left(x\right),x\right)}^{2}$
 ${\mathrm{ODE2}}{≔}{\mathrm{y\text{'}\text{'}}}{-}\frac{{y}{}\left({y}{}{x}{+}{4}\right)}{{2}{}{{x}}^{{2}}}{-}{y}{}{\mathrm{y\text{'}}}{-}\frac{{{\mathrm{y\text{'}}}}^{{2}}{}{x}}{{2}}$ (7)
 > $\mathrm{\mu }≔\mathrm{intfactor}\left(\mathrm{ODE2}\right)$
 ${\mathrm{\mu }}{≔}{{x}}^{{2}}{}{{ⅇ}}^{{-}\frac{{y}{}{x}}{{2}}}{,}\frac{{{ⅇ}}^{{-}\frac{{y}{}{x}}{{2}}}}{{x}}$ (8)
 > $\mathrm{odeadvisor}\left(\mathrm{\mu }\left[1\right]\mathrm{ODE2}\right)$
 $\left[\left[{\mathrm{_2nd_order}}{,}{\mathrm{_exact}}{,}{\mathrm{_nonlinear}}\right]{,}\left[{\mathrm{_2nd_order}}{,}{\mathrm{_with_linear_symmetries}}\right]{,}\left[{\mathrm{_2nd_order}}{,}{\mathrm{_reducible}}{,}{\mathrm{_mu_xy}}\right]\right]$ (9)
 > $\mathrm{firint}\left(\mathrm{\mu }\left[1\right]\mathrm{ODE2}\right)$
 ${{ⅇ}}^{{-}\frac{{y}{}{x}}{{2}}}{}\left({y}{}{x}{+}{6}\right){+}{{x}}^{{2}}{}{{ⅇ}}^{{-}\frac{{y}{}{x}}{{2}}}{}{\mathrm{y\text{'}}}{+}{\mathrm{_C1}}{=}{0}$ (10)
 > $\mathrm{firint}\left(\mathrm{\mu }\left[2\right]\mathrm{ODE2}\right)$
 $\frac{{{ⅇ}}^{{-}\frac{{y}{}{x}}{{2}}}{}{y}}{{{x}}^{{2}}}{+}\frac{{{ⅇ}}^{{-}\frac{{y}{}{x}}{{2}}}{}{\mathrm{y\text{'}}}}{{x}}{+}{\mathrm{_C1}}{=}{0}$ (11)

Eliminating $y'$ from these two first integrals and renaming  $\mathrm{_C1}\to \mathrm{_C2}$ in one of them we obtain the ODE solution; this process is performed internally by dsolve to arrive at:

 > $\mathrm{dsolve}\left(\mathrm{ODE2}\right)$
 ${y}{=}{-}\frac{{2}{}{\mathrm{ln}}{}\left(\frac{{\mathrm{_C2}}{}{{x}}^{{3}}}{{6}}{-}\frac{{\mathrm{_C1}}}{{6}}\right)}{{x}}$ (12)

A third order ODE example and two integrating factors for it:

 > $\mathrm{ODE3}≔\mathrm{diff}\left(y\left(x\right),x,x,x\right)=\frac{\left(\mathrm{diff}\left(y\left(x\right),x\right)x+y\left(x\right)\right)\mathrm{diff}\left(y\left(x\right),x,x\right)}{y\left(x\right)x}$
 ${\mathrm{ODE3}}{≔}{\mathrm{y\text{'}\text{'}\text{'}}}{=}\frac{\left({\mathrm{y\text{'}}}{}{x}{+}{y}\right){}{\mathrm{y\text{'}\text{'}}}}{{y}{}{x}}$ (13)
 > $\mathrm{\mu }≔\mathrm{intfactor}\left(\mathrm{ODE3}\right)$
 ${\mathrm{\mu }}{≔}\frac{{1}}{{y}{}{x}}{,}\frac{{1}}{{\mathrm{y\text{'}\text{'}}}}$ (14)
 > $\mathrm{odeadvisor}\left(\frac{1}{y\left(x\right)x}\mathrm{ODE3}\right)$
 $\left[\left[{\mathrm{_3rd_order}}{,}{\mathrm{_exact}}{,}{\mathrm{_nonlinear}}\right]{,}\left[{\mathrm{_3rd_order}}{,}{\mathrm{_with_linear_symmetries}}\right]{,}\left[{\mathrm{_3rd_order}}{,}{\mathrm{_reducible}}{,}{\mathrm{_mu_y2}}\right]{,}\left[{\mathrm{_3rd_order}}{,}{\mathrm{_reducible}}{,}{\mathrm{_mu_poly_yn}}\right]\right]$ (15)

For ODE3 above, using its defaults intfactor found two integrating factors. If however we use the try_hard=true option, an integrating factor depending on an arbitrary function   is obtained.

 > $\mathrm{intfactor}\left(\mathrm{ODE3},\mathrm{try_hard}\right)$
 $\frac{{1}}{{\mathrm{y\text{'}\text{'}}}}{,}\frac{{\mathrm{_F1}}{}\left(\frac{{\mathrm{y\text{'}\text{'}}}}{{y}{}{x}}\right)}{{y}{}{x}}$ (16)

Integrating factors can be tested using mutest, which returns zero when the integrating factor is correct, meaning that it cancels the appropriate integrating factor PDE system (see odepde).

 > $\mathrm{map}\left(\mathrm{DEtools}\left[\mathrm{mutest}\right],\left[\right],\mathrm{ODE3}\right)$
 $\left[{0}{,}{0}\right]$ (17)

A fourth order ODE example and an integrating factor polynomial in the highest possible derivative $y\text{'}\text{'}\text{'}$ (for an ODE of order n, mu can depend at most on the nth-1 derivative).

 > $\mathrm{ODE4}≔\mathrm{diff}\left(y\left(x\right),x,x,x,x\right)=\frac{\frac{1}{2}\left(\left(-{x}^{2}\mathrm{diff}\left(y\left(x\right),x,x\right){\mathrm{diff}\left(y\left(x\right),x,x,x\right)}^{2}{\mathrm{ln}\left(x\right)}^{2}-\mathrm{diff}\left(y\left(x\right),x\right){\mathrm{diff}\left(y\left(x\right),x,x,x\right)}^{2}{\mathrm{ln}\left(x\right)}^{2}x\right){y\left(x\right)}^{3}-{x}^{2}{\mathrm{diff}\left(y\left(x\right),x\right)}^{2}{\mathrm{diff}\left(y\left(x\right),x,x,x\right)}^{2}{\mathrm{ln}\left(x\right)}^{2}{y\left(x\right)}^{2}+2y\left(x\right)+2\mathrm{diff}\left(y\left(x\right),x\right)\mathrm{ln}\left(x\right)x\right)}{{\mathrm{ln}\left(x\right)}^{2}{y\left(x\right)}^{3}{x}^{2}\mathrm{diff}\left(y\left(x\right),x\right)\mathrm{diff}\left(y\left(x\right),x,x,x\right)}$
 ${\mathrm{ODE4}}{≔}{\mathrm{y\text{'}\text{'}\text{'}\text{'}}}{=}\frac{\left({-}{{x}}^{{2}}{}{\mathrm{y\text{'}\text{'}}}{}{{\mathrm{y\text{'}\text{'}\text{'}}}}^{{2}}{}{{\mathrm{ln}}{}\left({x}\right)}^{{2}}{-}{\mathrm{y\text{'}}}{}{{\mathrm{y\text{'}\text{'}\text{'}}}}^{{2}}{}{{\mathrm{ln}}{}\left({x}\right)}^{{2}}{}{x}\right){}{{y}}^{{3}}{-}{{x}}^{{2}}{}{{\mathrm{y\text{'}}}}^{{2}}{}{{\mathrm{y\text{'}\text{'}\text{'}}}}^{{2}}{}{{\mathrm{ln}}{}\left({x}\right)}^{{2}}{}{{y}}^{{2}}{+}{2}{}{y}{+}{2}{}{\mathrm{y\text{'}}}{}{\mathrm{ln}}{}\left({x}\right){}{x}}{{2}{}{{\mathrm{ln}}{}\left({x}\right)}^{{2}}{}{{y}}^{{3}}{}{{x}}^{{2}}{}{\mathrm{y\text{'}}}{}{\mathrm{y\text{'}\text{'}\text{'}}}}$ (18)
 > $\mathrm{\mu }≔\mathrm{intfactor}\left(\mathrm{ODE4}\right)$
 ${\mathrm{\mu }}{≔}{\mathrm{y\text{'}}}{}{\mathrm{y\text{'}\text{'}\text{'}}}{}{y}{}{x}$ (19)

A second order ODE with three integrating factors:

 > $\mathrm{ode}≔\mathrm{diff}\left(\mathrm{diff}\left(y\left(x\right),x\right),x\right)=-\frac{\frac{1}{x}\mathrm{diff}\left(y\left(x\right),x\right)\left(2y\left(x\right)+x\mathrm{diff}\left(y\left(x\right),x\right)\right)}{y\left(x\right)}$
 ${\mathrm{ode}}{≔}{\mathrm{y\text{'}\text{'}}}{=}{-}\frac{{\mathrm{y\text{'}}}{}\left({\mathrm{y\text{'}}}{}{x}{+}{2}{}{y}\right)}{{x}{}{y}}$ (20)
 > $\mathrm{\mu }≔\mathrm{intfactor}\left(\mathrm{ode}\right)$
 ${\mathrm{\mu }}{≔}{y}{}{{x}}^{{2}}{,}{y}{}{x}{,}\frac{{1}}{{\mathrm{y\text{'}}}}$ (21)

Of the three integrating factors found, only two of them lead to different reductions of order. The reduction obtained using ${\mathrm{\mu }}_{1}$ is as follows.

 > $\mathrm{R1}≔\mathrm{firint}\left(\mathrm{\mu }\left[1\right]\mathrm{ode},y\left(x\right)\right)$
 ${\mathrm{R1}}{≔}{y}{}{\mathrm{y\text{'}}}{}{{x}}^{{2}}{+}{\mathrm{_C1}}{=}{0}$ (22)
 > $\mathrm{explicit_R1}≔\mathrm{diff}\left(y\left(x\right),x\right)=\mathrm{solve}\left(\mathrm{R1},\mathrm{diff}\left(y\left(x\right),x\right)\right)$
 ${\mathrm{explicit_R1}}{≔}{\mathrm{y\text{'}}}{=}{-}\frac{{\mathrm{_C1}}}{{{x}}^{{2}}{}{y}}$ (23)

The reduction obtained using ${\mathrm{\mu }}_{2}$:

 > $\mathrm{R2}≔\mathrm{firint}\left(\mathrm{\mu }\left[2\right]\mathrm{ode},y\left(x\right)\right)$
 ${\mathrm{R2}}{≔}\frac{{1}}{{2}}{}{{y}}^{{2}}{+}{\mathrm{y\text{'}}}{}{x}{}{y}{+}{\mathrm{_C1}}{=}{0}$ (24)
 > $\mathrm{explicit_R2}≔\mathrm{diff}\left(y\left(x\right),x\right)=\mathrm{solve}\left(\mathrm{R2},\mathrm{diff}\left(y\left(x\right),x\right)\right)$
 ${\mathrm{explicit_R2}}{≔}{\mathrm{y\text{'}}}{=}{-}\frac{{{y}}^{{2}}{+}{2}{}{\mathrm{_C1}}}{{2}{}{y}{}{x}}$ (25)

The reduction obtained using ${\mathrm{\mu }}_{3}$:

 > $\mathrm{R3}≔\mathrm{firint}\left(\mathrm{\mu }\left[3\right]\mathrm{ode},y\left(x\right)\right)$
 ${\mathrm{R3}}{≔}{\mathrm{ln}}{}\left({\mathrm{y\text{'}}}\right){+}{\mathrm{ln}}{}\left({y}\right){+}{2}{}{\mathrm{ln}}{}\left({x}\right){+}{\mathrm{_C1}}{=}{0}$ (26)
 > $\mathrm{explicit_R3}≔\mathrm{diff}\left(y\left(x\right),x\right)=\mathrm{solve}\left(\mathrm{R3},\mathrm{diff}\left(y\left(x\right),x\right)\right)$
 ${\mathrm{explicit_R3}}{≔}{\mathrm{y\text{'}}}{=}\frac{{1}}{{{ⅇ}}^{{\mathrm{_C1}}}{}{y}{}{{x}}^{{2}}}$ (27)

When the algorithms implemented in intfactor fail in computing integrating factors - or to investigate the existence of more or different type of integrating factors - we can give a hint to intfactor by using the option _mu = some_mu_form, where some_mu_form contains some unknowns - functions of the ODE variables or just symbols - whose values are to be determined. For example, for the following ODE we use intfactor to compute - if possible - integrating factors of the form $xf\left(y,y'\right)$ for some function $f\left(y,y'\right)$ to be determined.

 > $\mathrm{ode}≔\mathrm{diff}\left(y\left(x\right),x,x\right)=\frac{\left(y\left(x\right)+x\mathrm{diff}\left(y\left(x\right),x\right)+2xy\left(x\right)\mathrm{diff}\left(y\left(x\right),x\right)+{x}^{2}{\mathrm{diff}\left(y\left(x\right),x\right)}^{2}\right)\mathrm{diff}\left(y\left(x\right),x\right)}{y\left(x\right)x}$
 ${\mathrm{ode}}{≔}{\mathrm{y\text{'}\text{'}}}{=}\frac{\left({{x}}^{{2}}{}{{\mathrm{y\text{'}}}}^{{2}}{+}{2}{}{\mathrm{y\text{'}}}{}{x}{}{y}{+}{\mathrm{y\text{'}}}{}{x}{+}{y}\right){}{\mathrm{y\text{'}}}}{{x}{}{y}}$ (28)
 > $\mathrm{intfactor}\left(\mathrm{ode},\mathrm{_μ}=xf\left(y,\mathrm{_y1}\right)\right)$
 $\frac{{x}{}{y}}{{{\mathrm{y\text{'}}}}^{{2}}}{,}\frac{{x}{}{y}{}{{\mathrm{Ei}}}_{{1}}{}\left({-}{2}{}{y}\right)}{{{\mathrm{y\text{'}}}}^{{2}}}$ (29)
 > $\mathrm{map}\left(\mathrm{DEtools}\left[\mathrm{mutest}\right],\left[\right],\mathrm{ode}\right)$
 $\left[{0}{,}{0}\right]$ (30)

In this example, we indicate a form for the integrating factor which depends on 4 symbolic constants $\mathrm{k1}$ to $\mathrm{k4}$ to be determined.

 > $\mathrm{ode}≔\mathrm{diff}\left(y\left(x\right),\mathrm{}\left(x,4\right)\right)=-\frac{2{\mathrm{diff}\left(y\left(x\right),x\right)}^{2}\mathrm{diff}\left(y\left(x\right),\mathrm{}\left(x,3\right)\right)+{\mathrm{diff}\left(y\left(x\right),x\right)}^{3}\mathrm{diff}\left(y\left(x\right),\mathrm{}\left(x,3\right)\right)+\mathrm{diff}\left(y\left(x\right),\mathrm{}\left(x,3\right)\right)\mathrm{diff}\left(y\left(x\right),x,x\right){\mathrm{diff}\left(y\left(x\right),x\right)}^{2}+y\left(x\right)\mathrm{diff}\left(y\left(x\right),x\right)+{\mathrm{diff}\left(y\left(x\right),x\right)}^{2}x-\mathrm{diff}\left(y\left(x\right),x,x\right)y\left(x\right)x}{{\mathrm{diff}\left(y\left(x\right),x\right)}^{2}\left(x+y\left(x\right)+\mathrm{diff}\left(y\left(x\right),x\right)\right)}$
 ${\mathrm{ode}}{≔}{\mathrm{y\text{'}\text{'}\text{'}\text{'}}}{=}{-}\frac{{{\mathrm{y\text{'}}}}^{{3}}{}{\mathrm{y\text{'}\text{'}\text{'}}}{+}{\mathrm{y\text{'}\text{'}\text{'}}}{}{\mathrm{y\text{'}\text{'}}}{}{{\mathrm{y\text{'}}}}^{{2}}{-}{y}{}{x}{}{\mathrm{y\text{'}\text{'}}}{+}{{\mathrm{y\text{'}}}}^{{2}}{}{x}{+}{2}{}{{\mathrm{y\text{'}}}}^{{2}}{}{\mathrm{y\text{'}\text{'}\text{'}}}{+}{y}{}{\mathrm{y\text{'}}}}{{{\mathrm{y\text{'}}}}^{{2}}{}\left({x}{+}{y}{+}{\mathrm{y\text{'}}}\right)}$ (31)
 > $\mathrm{hint}≔\mathrm{_μ}=\mathrm{k1}x+\mathrm{k2}y+\mathrm{k3}\mathrm{_y1}+\mathrm{k4}\mathrm{_y2}+\mathrm{k5}\mathrm{_y3}$
 ${\mathrm{hint}}{≔}{\mathrm{_μ}}{=}{\mathrm{k3}}{}{\mathrm{_y1}}{+}{\mathrm{k4}}{}{\mathrm{_y2}}{+}{\mathrm{k5}}{}{\mathrm{_y3}}{+}{\mathrm{k1}}{}{x}{+}{\mathrm{k2}}{}{y}$ (32)
 > $\mathrm{intfactor}\left(\mathrm{ode},\mathrm{hint}\right)$
 ${x}{+}{y}{+}{\mathrm{y\text{'}}}$ (33)