polysols - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

DEtools

 polysols
 find polynomial solutions of a linear ODE

 Calling Sequence polysols(ODE, v, opts) polysols(coeff_list, g, x, output=solution)

Parameters

 ODE - ordinary differential equation - linear or not - or a set or list of them v - optional - the unknown function of the ODE, or a set or list of them in the case of an ODE system coeff_list - list of coefficients of a Linear ODE, or list of lists in the case of a Linear ODE system g - right-hand side of a Linear ODE or a list of them in the case of a Linear ODE system x - independent variable of a Linear ODE or system of them output = solution - optional, for Linear ODEs, the output is returned as a solution expression, as dsolve does opts - optional, for Nonlinear ODEs and systems of them, could be any of the optional arguments accepted by casesplit

Description

 • The polysols command returns polynomial solutions of an ODE or a system of them with rational function coefficients.
 • The first calling sequence has an ODE or a system of them as the first argument, and, optionally, a variable, set or list of them indicating who are the unknown functions as the second argument - say v. If v is not given, all differentiated functions found in ODE are taken as the unknown functions of the problem. In the case of Nonlinear ODE systems, v can also be a solving ranking as described in casesplit and any of the optional arguments accepted by the casesplit command are valid for polysols.
 • The second calling sequence is valid only for Linear ODEs and is described as follows.
 Case of a single Linear ODE
 - The first argument is the list of coefficients ${c}_{n}\left(x\right)$

$\left[{c}_{0},{c}_{1},\mathrm{...},{c}_{n}\right]$

 entering the ODE, say, in y(x),

${c}_{0}y+{c}_{1}\mathrm{y\text{'}}+\mathrm{...}+{c}_{n}{y}^{\left(n\right)}=g\left(x\right)$

 - The second argument is the right-hand side of such an equation, g(x), so-called non-homogeneous term;
 - The third argument is the independent variable of the ODE, say x.
 • The output is a list of independent polynomial solutions representing a solution basis. If the ODE admits as many polynomial solutions as the differential order, the basis returned is complete. In the nonhomogeneous case, the returned value is a two-element list, with the first element a basis for the homogeneous case and the second element a particular polynomial solution (if it exists). When the optional argument output=solution is given, the output is an explicit closed form solution as the ones returned by dsolve.
 Case of a system of Linear ODEs
 - The first argument is a list of list(s) where each inner list corresponds to one equation in the system; in turn this list contains list(s) of the coefficients corresponding to each function and its derivatives in the linear ODE system, as explained in the previous paragraph for a single Linear ODE. For example, suppose a system with unknowns ${f}_{i}\left(x\right)$, this shows the correspondence between the list and function representations in typical cases (see also the examples at the end)

$\left[\left[0,1\right],\left[-1\right]\right]=\left(\frac{ⅆ}{ⅆx}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{f}_{1}\left(x\right)-{f}_{2}\left(x\right)=0\right)$

$\left[\left[\right],\left[0,1\right],\left[-1\right],\left[-1\right]\right]=\left(\frac{ⅆ}{ⅆx}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{f}_{2}\left(x\right)-{f}_{3}\left(x\right)-{f}_{4}\left(x\right)=0\right)$

$\left[\left[\right],\left[-{x}^{2}\right],\left[2\right],\left[\right],\left[\right],\left[0,1\right]\right]=\left(-{x}^{2}{f}_{2}\left(x\right)+2{f}_{3}\left(x\right)+\frac{ⅆ}{ⅆx}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{f}_{6}\left(x\right)=0\right)$

 - The second argument is a list containing the right-hand sides of each equation of the system
 - The third argument is the independent variable of the Linear ODE system
 • The output consists of a list of solution basis corresponding to each function. When no polynomial solutions exist for one function, the corresponding list is empty. In the nonhomogeneous case, the output consists of a two-element list, with the first element a list of lists of solution basis for the homogeneous case and the second element a particular polynomial solution (if it exists). When the optional argument output=solution is given, the output is an explicit closed form solution as the ones returned by dsolve, enclosed into a list.
 • The polysols command is implemented as a module and provides direct access to each of its subroutines. These subroutines are of use as programmer entry points by shortcutting the identification of the type of problem and testing of arguments for correctness, e.g., when you know a priori what the type of problem is, and are sure that test for correctness of arguments is not necessary. For each type of problem you can use the corresponding subroutine as follows.
 - Problem: an ODE system, Linear or not, represented by standard equations. Subroutine: DEtools/polysols:-ODESystem
 - Problem: a Linear ODE system represented by a list of lists of lists. Subroutine: DEtools/polysols:-ODESystemList
 - Problem: a single ODE, Linear or not, represented by an equation or algebraic expression: DEtools/polysols:-ODE
 - Problem: a single Linear ODE  represented by a list of coefficients. Subroutine: DEtools/polysols:-ODEList
 - Problem: particular solution for a single Linear ODE represented by a list of coefficients. Subroutine: DEtools/polysols:-ODEListParticularSolution
 - Problem: switch the output from a list consisting of a solution basis to the standard solution form used by dsolve. Subroutine: DEtools/polysols:-FormatOutput
 Libraries in use
 • In order to find solutions:
 - In the Linear ODE system case, or the case of a single ODE enclosed into a list, polysols uses the LinearFunctionalSystems[PolynomialSolution] command;
 - In the single Linear ODE case, polysols uses its own optimized code;
 - In the case of a Nonlinear ODE or a system of them, polysols uses the libraries for computing traveling wave solutions for autonomous PDE systems.

Examples

 > $\mathrm{with}\left(\mathrm{DEtools},\mathrm{polysols}\right):$

Linear ODEs

 Homogeneous case
 > $\mathrm{odeH}≔\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}z\left(t\right)-\frac{3\left(\frac{ⅆ}{ⅆt}z\left(t\right)\right)}{t}+\frac{3z\left(t\right)}{{t}^{2}}=0$
 ${\mathrm{odeH}}{≔}\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{z}{}\left({t}\right){-}\frac{{3}{}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{z}{}\left({t}\right)\right)}{{t}}{+}\frac{{3}{}{z}{}\left({t}\right)}{{{t}}^{{2}}}{=}{0}$ (1)
 > $\mathrm{polysols}\left(\mathrm{odeH}\right)$
 $\left[{t}{,}{{t}}^{{3}}\right]$ (2)
 The same result can be obtained using the programmer's entry point, DEtools/polysols:-ODE(ode, z(t)), which shortcuts the identification of the type of problem and avoids spending time testing the arguments for correctness.
 Same problem but using the optional argument output = solution
 > $\mathrm{polysols}\left(\mathrm{odeH},\mathrm{output}=\mathrm{solution}\right)$
 $\left[{z}{}\left({t}\right){=}{\mathrm{_C2}}{}{{t}}^{{3}}{+}{\mathrm{_C1}}{}{t}\right]$ (3)
 Adding a non-homogeneous term
 > $\mathrm{odeNH}≔\mathrm{lhs}\left(\mathrm{odeH}\right)=3{t}^{2}{\left(t+1\right)}^{2}$
 ${\mathrm{odeNH}}{≔}\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{z}{}\left({t}\right){-}\frac{{3}{}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{z}{}\left({t}\right)\right)}{{t}}{+}\frac{{3}{}{z}{}\left({t}\right)}{{{t}}^{{2}}}{=}{3}{}{{t}}^{{2}}{}{\left({t}{+}{1}\right)}^{{2}}$ (4)
 > $\mathrm{polysols}\left(\mathrm{odeNH}\right)$
 $\left[\left[{t}{,}{{t}}^{{3}}\right]{,}\frac{{1}}{{5}}{}{{t}}^{{6}}{+}\frac{{3}}{{4}}{}{{t}}^{{5}}{+}{{t}}^{{4}}\right]$ (5)
 The same Linear ODE but given as a list of coefficients (see DEtools[convertAlg]), with the right-hand-side (non-homogeneous term) passed as second argument to polysols
 > $\mathrm{odeH_list},\mathrm{NHterm}≔\mathrm{op}\left({\mathrm{DEtools}}_{\mathrm{convertAlg}}\left(\mathrm{odeNH},z\left(t\right)\right)\right)$
 ${\mathrm{odeH_list}}{,}{\mathrm{NHterm}}{≔}\left[\frac{{3}}{{{t}}^{{2}}}{,}{-}\frac{{3}}{{t}}{,}{1}\right]{,}{3}{}{{t}}^{{4}}{+}{6}{}{{t}}^{{3}}{+}{3}{}{{t}}^{{2}}$ (6)
 > $\mathrm{polysols}\left(\mathrm{odeH_list},\mathrm{NHterm},t\right)$
 $\left[\left[{t}{,}{{t}}^{{3}}\right]{,}\frac{{1}}{{5}}{}{{t}}^{{6}}{+}\frac{{3}}{{4}}{}{{t}}^{{5}}{+}{{t}}^{{4}}\right]$ (7)
 When giving a Linear ODE system, the number of equations and unknowns must be the same. The default format for the output is a list of polynomial solution basis for each function
 > $\mathrm{sys}≔\left[\frac{ⅆ}{ⅆx}\mathrm{y1}\left(x\right)-\mathrm{y2}\left(x\right),\frac{ⅆ}{ⅆx}\mathrm{y2}\left(x\right)-\mathrm{y3}\left(x\right)-\mathrm{y4}\left(x\right),\frac{ⅆ}{ⅆx}\mathrm{y3}\left(x\right)-\mathrm{y5}\left(x\right),\frac{ⅆ}{ⅆx}\mathrm{y4}\left(x\right)-2\mathrm{y1}\left(x\right)-2x\mathrm{y2}\left(x\right)-\mathrm{y5}\left(x\right),\frac{ⅆ}{ⅆx}\mathrm{y5}\left(x\right)-{x}^{2}\mathrm{y1}\left(x\right)-2x\mathrm{y3}\left(x\right)-\mathrm{y6}\left(x\right),\frac{ⅆ}{ⅆx}\mathrm{y6}\left(x\right)-{x}^{2}\mathrm{y2}\left(x\right)+2\mathrm{y3}\left(x\right)\right]$
 ${\mathrm{sys}}{≔}\left[\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{y1}}{}\left({x}\right){-}{\mathrm{y2}}{}\left({x}\right){,}\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{y2}}{}\left({x}\right){-}{\mathrm{y3}}{}\left({x}\right){-}{\mathrm{y4}}{}\left({x}\right){,}\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{y3}}{}\left({x}\right){-}{\mathrm{y5}}{}\left({x}\right){,}\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{y4}}{}\left({x}\right){-}{2}{}{\mathrm{y1}}{}\left({x}\right){-}{2}{}{x}{}{\mathrm{y2}}{}\left({x}\right){-}{\mathrm{y5}}{}\left({x}\right){,}\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{y5}}{}\left({x}\right){-}{{x}}^{{2}}{}{\mathrm{y1}}{}\left({x}\right){-}{2}{}{x}{}{\mathrm{y3}}{}\left({x}\right){-}{\mathrm{y6}}{}\left({x}\right){,}\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{y6}}{}\left({x}\right){-}{{x}}^{{2}}{}{\mathrm{y2}}{}\left({x}\right){+}{2}{}{\mathrm{y3}}{}\left({x}\right)\right]$ (8)
 > $\mathrm{vars}≔\left[\mathrm{y1}\left(x\right),\mathrm{y2}\left(x\right),\mathrm{y3}\left(x\right),\mathrm{y4}\left(x\right),\mathrm{y5}\left(x\right),\mathrm{y6}\left(x\right)\right]$
 ${\mathrm{vars}}{≔}\left[{\mathrm{y1}}{}\left({x}\right){,}{\mathrm{y2}}{}\left({x}\right){,}{\mathrm{y3}}{}\left({x}\right){,}{\mathrm{y4}}{}\left({x}\right){,}{\mathrm{y5}}{}\left({x}\right){,}{\mathrm{y6}}{}\left({x}\right)\right]$ (9)
 > $\mathrm{polysols}\left(\mathrm{sys},\mathrm{vars}\right)$
 $\left[\left[{0}{,}{1}{,}{x}\right]{,}\left[{0}{,}{0}{,}{1}\right]{,}\left[{1}{,}{-}{x}{,}{-}{{x}}^{{2}}\right]{,}\left[{-1}{,}{x}{,}{{x}}^{{2}}\right]{,}\left[{0}{,}{-1}{,}{-}{2}{}{x}\right]{,}\left[{-}{2}{}{x}{,}{{x}}^{{2}}{,}{{x}}^{{3}}{-}{2}\right]\right]$ (10)
 The same result can be obtained using the corresponding programmer's entry point, DEtools/polysols:-ODESystem(sys, vars).
 The equivalent closed form solution using the output = solution option
 > $\mathrm{polysols}\left(\mathrm{sys},\mathrm{vars},\mathrm{output}=\mathrm{solution}\right)$
 $\left[\left\{{\mathrm{y1}}{}\left({x}\right){=}{\mathrm{_C3}}{}{x}{+}{\mathrm{_C2}}{,}{\mathrm{y2}}{}\left({x}\right){=}{\mathrm{_C3}}{,}{\mathrm{y3}}{}\left({x}\right){=}{-}{\mathrm{_C3}}{}{{x}}^{{2}}{-}{\mathrm{_C2}}{}{x}{+}{\mathrm{_C1}}{,}{\mathrm{y4}}{}\left({x}\right){=}{\mathrm{_C3}}{}{{x}}^{{2}}{+}{\mathrm{_C2}}{}{x}{-}{\mathrm{_C1}}{,}{\mathrm{y5}}{}\left({x}\right){=}{-}{2}{}{\mathrm{_C3}}{}{x}{-}{\mathrm{_C2}}{,}{\mathrm{y6}}{}\left({x}\right){=}{\mathrm{_C3}}{}{{x}}^{{3}}{+}{\mathrm{_C2}}{}{{x}}^{{2}}{-}{2}{}{\mathrm{_C1}}{}{x}{-}{2}{}{\mathrm{_C3}}\right\}\right]$ (11)
 The same Linear ODE system given as a list of lists of lists
 > $\mathrm{sys}≔\left[\left[\left[0,1\right],\left[-1\right]\right],\left[\left[\right],\left[0,1\right],\left[-1\right],\left[-1\right]\right],\left[\left[\right],\left[\right],\left[0,1\right],\left[\right],\left[-1\right]\right],\left[\left[-2\right],\left[-2x\right],\left[\right],\left[0,1\right],\left[-1\right]\right],\left[\left[-{x}^{2}\right],\left[\right],\left[-2x\right],\left[\right],\left[0,1\right],\left[-1\right]\right],\left[\left[\right],\left[-{x}^{2}\right],\left[2\right],\left[\right],\left[\right],\left[0,1\right]\right]\right]$
 ${\mathrm{sys}}{≔}\left[\left[\left[{0}{,}{1}\right]{,}\left[{-1}\right]\right]{,}\left[\left[\right]{,}\left[{0}{,}{1}\right]{,}\left[{-1}\right]{,}\left[{-1}\right]\right]{,}\left[\left[\right]{,}\left[\right]{,}\left[{0}{,}{1}\right]{,}\left[\right]{,}\left[{-1}\right]\right]{,}\left[\left[{-2}\right]{,}\left[{-}{2}{}{x}\right]{,}\left[\right]{,}\left[{0}{,}{1}\right]{,}\left[{-1}\right]\right]{,}\left[\left[{-}{{x}}^{{2}}\right]{,}\left[\right]{,}\left[{-}{2}{}{x}\right]{,}\left[\right]{,}\left[{0}{,}{1}\right]{,}\left[{-1}\right]\right]{,}\left[\left[\right]{,}\left[{-}{{x}}^{{2}}\right]{,}\left[{2}\right]{,}\left[\right]{,}\left[\right]{,}\left[{0}{,}{1}\right]\right]\right]$ (12)
 > $\mathrm{polysols}\left(\mathrm{sys},\left[0,0,0,0,0,0\right],x\right)$
 $\left[\left[{0}{,}{1}{,}{x}\right]{,}\left[{0}{,}{0}{,}{1}\right]{,}\left[{1}{,}{-}{x}{,}{-}{{x}}^{{2}}\right]{,}\left[{-1}{,}{x}{,}{{x}}^{{2}}\right]{,}\left[{0}{,}{-1}{,}{-}{2}{}{x}\right]{,}\left[{-}{2}{}{x}{,}{{x}}^{{2}}{,}{{x}}^{{3}}{-}{2}\right]\right]$ (13)
 Using polysols to reduce the order of an ODE: combining the functionality of various DEtools commands
 Consider the list of four infinitesimals for symmetry generators of the form $\left[\mathrm{\xi }=0,\mathrm{\eta }=f\left(x\right)\right]$ (see symgen)
 > $S≔\left[\left[0,x\right],\left[0,{\left(x-1\right)}^{2}\right],\left[0,{x}^{3}\right],\left[0,\frac{1}{x}\right]\right]$
 ${S}{≔}\left[\left[{0}{,}{x}\right]{,}\left[{0}{,}{\left({x}{-}{1}\right)}^{{2}}\right]{,}\left[{0}{,}{{x}}^{{3}}\right]{,}\left[{0}{,}\frac{{1}}{{x}}\right]\right]$ (14)
 For linear ODEs, the function $f\left(x\right)$ entering symmetries of the form $\left[0,f\left(x\right)\right]$ is always a solution of the ODE. So we can construct the most general Linear ODE of fourth order having the four solutions (symmetries) above, three of which are polynomial in $x$, via (see DEtools[equinv])
 > $\mathrm{lin_DE}≔{\mathrm{DEtools}}_{\mathrm{equinv}}\left(S,y\left(x\right),4\right)$
 ${\mathrm{lin_DE}}{≔}\frac{{{ⅆ}}^{{4}}}{{ⅆ}{{x}}^{{4}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right){=}\frac{\left({-}{{x}}^{{2}}{+}{3}\right){}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right)}{{{x}}^{{3}}{-}{x}}{+}\frac{{3}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right)}{{{x}}^{{2}}}{-}\frac{{6}{}\left(\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right)}{{{x}}^{{3}}{-}{x}}{+}\frac{{6}{}{y}{}\left({x}\right)}{{{x}}^{{2}}{}\left({{x}}^{{2}}{-}{1}\right)}{+}{\mathrm{_F1}}{}\left({x}\right)$ (15)
 For the purpose of using polysols in this example, we remove here the arbitrary non-homogeneous and non-rational term $\mathrm{_F1}\left(x\right)$
 > $\mathrm{ode}≔\mathrm{subs}\left(\mathrm{_F1}\left(x\right)=0,\mathrm{lin_DE}\right)$
 ${\mathrm{ode}}{≔}\frac{{{ⅆ}}^{{4}}}{{ⅆ}{{x}}^{{4}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right){=}\frac{\left({-}{{x}}^{{2}}{+}{3}\right){}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right)}{{{x}}^{{3}}{-}{x}}{+}\frac{{3}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right)}{{{x}}^{{2}}}{-}\frac{{6}{}\left(\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right)}{{{x}}^{{3}}{-}{x}}{+}\frac{{6}{}{y}{}\left({x}\right)}{{{x}}^{{2}}{}\left({{x}}^{{2}}{-}{1}\right)}$ (16)
 Compute the polynomial solutions
 > $\mathrm{sol_basis}≔\mathrm{polysols}\left(\mathrm{ode},y\left(x\right)\right)$
 ${\mathrm{sol_basis}}{≔}\left[{x}{,}{{x}}^{{3}}{,}{{x}}^{{2}}{+}{1}\right]$ (17)
 Since this ODE is of fourth order and there are only three polynomial solutions, the above is not sufficient to construct the most general solution for ode. However, each ODE solution leads to a reduction of order (see DEtools[reduce_order]), so using the polynomial solutions above we can reduce the order by three
 > $R≔{\mathrm{DEtools}}_{\mathrm{reduce_order}}\left(\mathrm{ode},\mathrm{basis}=\mathrm{sol_basis}\right)$
 ${R}{≔}{y}{}\left({x}\right){=}\left({\mathrm{_a}}{}\left({\int }{2}{}{\mathrm{_a}}{}\left({\int }{-}\frac{\left({{\mathrm{_a}}}^{{2}}{-}{3}\right){}\left({\int }{\mathrm{_b}}{}\left({\mathrm{_a}}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{_a}}{+}{\mathrm{_C3}}\right)}{{2}{}{{\mathrm{_a}}}^{{4}}}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{_a}}{+}{\mathrm{_C2}}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{_a}}{+}{\mathrm{_C1}}\right)\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{&where}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\left[\left\{{-}{\mathrm{_a}}{}\left({{\mathrm{_a}}}^{{2}}{-}{1}\right){}\left({{\mathrm{_a}}}^{{2}}{-}{3}\right){}\left(\frac{{ⅆ}}{{ⅆ}{\mathrm{_a}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{_b}}{}\left({\mathrm{_a}}\right)\right){+}\left({-}{4}{}{{\mathrm{_a}}}^{{4}}{+}{6}{}{{\mathrm{_a}}}^{{2}}{-}{6}\right){}{\mathrm{_b}}{}\left({\mathrm{_a}}\right){=}{0}\right\}{,}\left\{{\mathrm{_a}}{=}{x}{,}{\mathrm{_b}}{}\left({\mathrm{_a}}\right){=}\frac{{x}{}\left({-}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right){}{{x}}^{{3}}{+}{3}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right){}{{x}}^{{2}}{-}{6}{}\left(\frac{{ⅆ}}{{ⅆ}{x}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right){}{x}{+}{3}{}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{x}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right)\right){}{x}{-}{3}{}\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{x}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{y}{}\left({x}\right){+}{6}{}{y}{}\left({x}\right)\right)}{{\left({{x}}^{{2}}{-}{3}\right)}^{{2}}}\right\}{,}\left\{{x}{=}{\mathrm{_a}}{,}{y}{}\left({x}\right){=}{\mathrm{_a}}{}\left({\int }{2}{}{\mathrm{_a}}{}\left({\int }{-}\frac{\left({{\mathrm{_a}}}^{{2}}{-}{3}\right){}\left({\int }{\mathrm{_b}}{}\left({\mathrm{_a}}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{_a}}{+}{\mathrm{_C3}}\right)}{{2}{}{{\mathrm{_a}}}^{{4}}}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{_a}}{+}{\mathrm{_C2}}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{_a}}{+}{\mathrm{_C1}}\right)\right\}\right]$ (18)
 In above we see an ODESolStructure where the solving of the original fourth order ODE problem is reduced to the solving of one first order linear ODE in $\mathrm{_b}\left(\mathrm{_a}\right)$
 > $\mathrm{op}\left(\left[2,2,1,1\right],R\right)$
 ${-}{\mathrm{_a}}{}\left({{\mathrm{_a}}}^{{2}}{-}{1}\right){}\left({{\mathrm{_a}}}^{{2}}{-}{3}\right){}\left(\frac{{ⅆ}}{{ⅆ}{\mathrm{_a}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{\mathrm{_b}}{}\left({\mathrm{_a}}\right)\right){+}\left({-}{4}{}{{\mathrm{_a}}}^{{4}}{+}{6}{}{{\mathrm{_a}}}^{{2}}{-}{6}\right){}{\mathrm{_b}}{}\left({\mathrm{_a}}\right){=}{0}$ (19)
 You can also request to reduce_order to automatically attempt integrating that remaining ODE using the option output = solution
 > $S≔{\mathrm{DEtools}}_{\mathrm{reduce_order}}\left(\mathrm{ode},\mathrm{basis}=\mathrm{sol_basis},\mathrm{output}=\mathrm{solution}\right)$
 ${S}{≔}{y}{}\left({x}\right){=}{x}{}\left({\int }{2}{}{x}{}\left({\int }{-}\frac{\left({{x}}^{{2}}{-}{3}\right){}\left({\int }\frac{\left({{x}}^{{2}}{-}{1}\right){}{\mathrm{_C4}}}{{{x}}^{{2}}{}{\left({{x}}^{{2}}{-}{3}\right)}^{{2}}}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{x}{+}{\mathrm{_C3}}\right)}{{2}{}{{x}}^{{4}}}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{x}{+}{\mathrm{_C2}}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{x}{+}{\mathrm{_C1}}\right)$ (20)
 > $\mathrm{value}\left(S\right)$
 ${y}{}\left({x}\right){=}{x}{}\left({\mathrm{_C2}}{}{{x}}^{{2}}{+}{\mathrm{_C3}}{}{x}{+}\frac{{\mathrm{_C3}}}{{x}}{+}\frac{{\mathrm{_C4}}}{{24}{}{{x}}^{{2}}}{+}{\mathrm{_C1}}\right)$ (21)
 In the solution above we see the four solutions used to construct the problem, three of which are polynomial and were used to reduce the problem from a fourth order ODE to a first order one.

Nonlinear ODEs

 A nonlinear fourth order equation, a constant and two linear solutions for it
 > $\mathrm{ODE}≔\left(\frac{1{q}^{2}{a}^{2}{\left(-1+{t}^{2}\right)}^{2}u\left(t\right)}{3s}+\frac{1pq{a}^{2}{\left(-1+{t}^{2}\right)}^{2}}{6s}\right){\left(\frac{ⅆ}{ⅆt}u\left(t\right)\right)}^{2}+\left(\frac{1{q}^{2}{a}^{2}\left(-1+{t}^{2}\right)t{u\left(t\right)}^{2}}{3s}+\frac{1pq{a}^{2}\left(-1+{t}^{2}\right)tu\left(t\right)}{3s}-\frac{1{a}^{2}\left(t-1\right)\left(3{6}^{\frac{1}{2}}arq\left({t}^{2}-\frac{1}{3}\right){\left(sq\right)}^{\frac{1}{2}}+\left(\left(12{q}^{2}{a}^{2}{t}^{2}-6{q}^{2}{a}^{2}-\frac{3{p}^{2}}{2}\right)s+{r}^{2}q\right)t\right)\left(t+1\right)}{18{s}^{2}}\right)\left(\frac{ⅆ}{ⅆt}u\left(t\right)\right)-\frac{1{a}^{3}{\left(-1+{t}^{2}\right)}^{3}\left(r{6}^{\frac{1}{2}}{\left(sq\right)}^{\frac{3}{2}}+12{q}^{2}at{s}^{2}\right)\left(\frac{{ⅆ}^{3}}{ⅆ{t}^{3}}u\left(t\right)\right)}{36{s}^{3}}+\left(\frac{1{q}^{2}{a}^{2}{\left(-1+{t}^{2}\right)}^{2}{u\left(t\right)}^{2}}{6s}+\frac{1pq{a}^{2}{\left(-1+{t}^{2}\right)}^{2}u\left(t\right)}{6s}-\frac{1{\left(-1+{t}^{2}\right)}^{2}{a}^{2}\left(12r{6}^{\frac{1}{2}}{\left(sq\right)}^{\frac{3}{2}}at+2s{r}^{2}q-3{s}^{2}{p}^{2}-12{q}^{2}{a}^{2}{s}^{2}+72{s}^{2}{q}^{2}{a}^{2}{t}^{2}\right)}{72{s}^{3}}\right)\left(\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}u\left(t\right)\right)-\frac{1{q}^{2}{a}^{4}{\left(-1+{t}^{2}\right)}^{4}\left(\frac{{ⅆ}^{4}}{ⅆ{t}^{4}}u\left(t\right)\right)}{36s}$
 ${\mathrm{ODE}}{≔}\left(\frac{{{q}}^{{2}}{}{{a}}^{{2}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}{}{u}{}\left({t}\right)}{{3}{}{s}}{+}\frac{{p}{}{q}{}{{a}}^{{2}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}}{{6}{}{s}}\right){}{\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}^{{2}}{+}\left(\frac{{{q}}^{{2}}{}{{a}}^{{2}}{}\left({{t}}^{{2}}{-}{1}\right){}{t}{}{{u}{}\left({t}\right)}^{{2}}}{{3}{}{s}}{+}\frac{{p}{}{q}{}{{a}}^{{2}}{}\left({{t}}^{{2}}{-}{1}\right){}{t}{}{u}{}\left({t}\right)}{{3}{}{s}}{-}\frac{{{a}}^{{2}}{}\left({t}{-}{1}\right){}\left({3}{}\sqrt{{6}}{}{a}{}{r}{}{q}{}\left({{t}}^{{2}}{-}\frac{{1}}{{3}}\right){}\sqrt{{s}{}{q}}{+}\left(\left({12}{}{{q}}^{{2}}{}{{a}}^{{2}}{}{{t}}^{{2}}{-}{6}{}{{q}}^{{2}}{}{{a}}^{{2}}{-}\frac{{3}}{{2}}{}{{p}}^{{2}}\right){}{s}{+}{{r}}^{{2}}{}{q}\right){}{t}\right){}\left({t}{+}{1}\right)}{{18}{}{{s}}^{{2}}}\right){}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right){-}\frac{{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{3}}{}\left({r}{}\sqrt{{6}}{}{\left({s}{}{q}\right)}^{{3}}{{2}}}{+}{12}{}{{q}}^{{2}}{}{a}{}{t}{}{{s}}^{{2}}\right){}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{t}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{36}{}{{s}}^{{3}}}{+}\left(\frac{{{q}}^{{2}}{}{{a}}^{{2}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}{}{{u}{}\left({t}\right)}^{{2}}}{{6}{}{s}}{+}\frac{{p}{}{q}{}{{a}}^{{2}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}{}{u}{}\left({t}\right)}{{6}{}{s}}{-}\frac{{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}{}{{a}}^{{2}}{}\left({12}{}{r}{}\sqrt{{6}}{}{\left({s}{}{q}\right)}^{{3}}{{2}}}{}{a}{}{t}{+}{2}{}{s}{}{{r}}^{{2}}{}{q}{-}{3}{}{{s}}^{{2}}{}{{p}}^{{2}}{-}{12}{}{{q}}^{{2}}{}{{a}}^{{2}}{}{{s}}^{{2}}{+}{72}{}{{s}}^{{2}}{}{{q}}^{{2}}{}{{a}}^{{2}}{}{{t}}^{{2}}\right)}{{72}{}{{s}}^{{3}}}\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right){-}\frac{{{q}}^{{2}}{}{{a}}^{{4}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{4}}{}\left(\frac{{{ⅆ}}^{{4}}}{{ⅆ}{{t}}^{{4}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{36}{}{s}}$ (22)
 > $\mathrm{ODE_sol}≔{\mathrm{DEtools}}_{\mathrm{polysols}}\left(\mathrm{ODE}\right)$
 ${\mathrm{ODE_sol}}{≔}\left[{u}{}\left({t}\right){=}{\mathrm{_C3}}{,}{u}{}\left({t}\right){=}{-}{a}{}{t}{-}\frac{\sqrt{{6}}{}\sqrt{{s}{}{q}}{}{r}{+}{3}{}{s}{}{p}}{{6}{}{s}{}{q}}{,}{u}{}\left({t}\right){=}{a}{}{t}{+}\frac{\sqrt{{6}}{}\sqrt{{s}{}{q}}{}{r}{-}{3}{}{s}{}{p}}{{6}{}{s}{}{q}}\right]$ (23)
 Solutions can be tested as usual with odetest
 > $\mathrm{map}\left(\mathrm{odetest},\mathrm{ODE_sol},\mathrm{ODE}\right)$
 $\left[{0}{,}{0}{,}{0}\right]$ (24)
 Using the optional argument singsol=false, the singular cases (in this example: the constant solution) are not computed
 > ${\mathrm{DEtools}}_{\mathrm{polysols}}\left(\mathrm{ODE},\mathrm{singsol}=\mathrm{false}\right)$
 $\left[{u}{}\left({t}\right){=}{-}{a}{}{t}{-}\frac{\sqrt{{6}}{}\sqrt{{s}{}{q}}{}{r}{+}{3}{}{s}{}{p}}{{6}{}{s}{}{q}}{,}{u}{}\left({t}\right){=}{a}{}{t}{+}\frac{\sqrt{{6}}{}\sqrt{{s}{}{q}}{}{r}{-}{3}{}{s}{}{p}}{{6}{}{s}{}{q}}\right]$ (25)
 A nonlinear seventh order equation depending on three parameters $\left\{a,b,c\right\}$ and polynomial solutions for it taking these parameters as solving variables:
 > $\mathrm{ODE}≔\left(\left(-a{t}^{2}+a\right)u\left(t\right)-\frac{1\left(13185{a}^{3}c{t}^{2}-4335{a}^{3}c+315{a}^{3}c{t}^{6}-525{a}^{3}c{t}^{4}+259080{a}^{5}{t}^{4}-259080{a}^{5}{t}^{2}+34544{a}^{5}+2159b\right)\left(t+1\right)\left(t-1\right)}{2159}\right)\left(\frac{ⅆ}{ⅆt}u\left(t\right)\right)-\frac{5{a}^{3}{\left(t+1\right)}^{4}{\left(t-1\right)}^{4}\left(105{t}^{2}c+17272{a}^{2}-35c\right)t\left(\frac{{ⅆ}^{4}}{ⅆ{t}^{4}}u\left(t\right)\right)}{4318}-\frac{1{a}^{3}{\left(t+1\right)}^{5}{\left(t-1\right)}^{5}\left(315{t}^{2}c+17272{a}^{2}-35c\right)\left(\frac{{ⅆ}^{5}}{ⅆ{t}^{5}}u\left(t\right)\right)}{17272}-\frac{15{a}^{3}\left(63{t}^{4}c-70{t}^{2}c+34544{a}^{2}{t}^{2}+879c-17272{a}^{2}\right)t{\left(t+1\right)}^{2}{\left(t-1\right)}^{2}\left(\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}u\left(t\right)\right)}{2159}-\frac{5{a}^{3}\left(315{t}^{4}c+103632{a}^{2}{t}^{2}-210{t}^{2}c+879c-17272{a}^{2}\right){\left(t+1\right)}^{3}{\left(t-1\right)}^{3}\left(\frac{{ⅆ}^{3}}{ⅆ{t}^{3}}u\left(t\right)\right)}{4318}-\frac{21c{a}^{3}{\left(-1+{t}^{2}\right)}^{6}t\left(\frac{{ⅆ}^{6}}{ⅆ{t}^{6}}u\left(t\right)\right)}{17272}-\frac{1c{a}^{3}{\left(-1+{t}^{2}\right)}^{7}\left(\frac{{ⅆ}^{7}}{ⅆ{t}^{7}}u\left(t\right)\right)}{34544}$
 ${\mathrm{ODE}}{≔}\left(\left({-}{a}{}{{t}}^{{2}}{+}{a}\right){}{u}{}\left({t}\right){-}\frac{\left({315}{}{{a}}^{{3}}{}{c}{}{{t}}^{{6}}{+}{259080}{}{{a}}^{{5}}{}{{t}}^{{4}}{-}{525}{}{{a}}^{{3}}{}{c}{}{{t}}^{{4}}{-}{259080}{}{{a}}^{{5}}{}{{t}}^{{2}}{+}{13185}{}{{a}}^{{3}}{}{c}{}{{t}}^{{2}}{+}{34544}{}{{a}}^{{5}}{-}{4335}{}{{a}}^{{3}}{}{c}{+}{2159}{}{b}\right){}\left({t}{+}{1}\right){}\left({t}{-}{1}\right)}{{2159}}\right){}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right){-}\frac{{5}{}{{a}}^{{3}}{}{\left({t}{+}{1}\right)}^{{4}}{}{\left({t}{-}{1}\right)}^{{4}}{}\left({105}{}{{t}}^{{2}}{}{c}{+}{17272}{}{{a}}^{{2}}{-}{35}{}{c}\right){}{t}{}\left(\frac{{{ⅆ}}^{{4}}}{{ⅆ}{{t}}^{{4}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{4318}}{-}\frac{{{a}}^{{3}}{}{\left({t}{+}{1}\right)}^{{5}}{}{\left({t}{-}{1}\right)}^{{5}}{}\left({315}{}{{t}}^{{2}}{}{c}{+}{17272}{}{{a}}^{{2}}{-}{35}{}{c}\right){}\left(\frac{{{ⅆ}}^{{5}}}{{ⅆ}{{t}}^{{5}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{17272}}{-}\frac{{15}{}{{a}}^{{3}}{}\left({63}{}{{t}}^{{4}}{}{c}{+}{34544}{}{{a}}^{{2}}{}{{t}}^{{2}}{-}{70}{}{{t}}^{{2}}{}{c}{-}{17272}{}{{a}}^{{2}}{+}{879}{}{c}\right){}{t}{}{\left({t}{+}{1}\right)}^{{2}}{}{\left({t}{-}{1}\right)}^{{2}}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{2159}}{-}\frac{{5}{}{{a}}^{{3}}{}\left({315}{}{{t}}^{{4}}{}{c}{+}{103632}{}{{a}}^{{2}}{}{{t}}^{{2}}{-}{210}{}{{t}}^{{2}}{}{c}{-}{17272}{}{{a}}^{{2}}{+}{879}{}{c}\right){}{\left({t}{+}{1}\right)}^{{3}}{}{\left({t}{-}{1}\right)}^{{3}}{}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{t}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{4318}}{-}\frac{{21}{}{c}{}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{6}}{}{t}{}\left(\frac{{{ⅆ}}^{{6}}}{{ⅆ}{{t}}^{{6}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{17272}}{-}\frac{{c}{}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{7}}{}\left(\frac{{{ⅆ}}^{{7}}}{{ⅆ}{{t}}^{{7}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right)}{{34544}}$ (26)
 > ${\mathrm{DEtools}}_{\mathrm{polysols}}\left(\mathrm{ODE},\mathrm{parameters}=\left\{a,b,c\right\}\right)$
 $\left[\left\{{a}{=}{a}{,}{b}{=}{b}{,}{c}{=}{c}{,}{u}{}\left({t}\right){=}{\mathrm{_C3}}\right\}{,}\left\{{a}{=}{a}{,}{b}{=}{-}\frac{{10952}}{{5}}{}{{a}}^{{5}}{-}{\mathrm{_C3}}{}{a}{,}{c}{=}{-}\frac{{2159}{}{{a}}^{{2}}}{{25}}{,}{u}{}\left({t}\right){=}\frac{{8316}}{{5}}{}{{a}}^{{4}}{}{{t}}^{{6}}{-}\frac{{33264}}{{5}}{}{{a}}^{{4}}{}{{t}}^{{4}}{+}{8316}{}{{a}}^{{4}}{}{{t}}^{{2}}{+}{\mathrm{_C3}}\right\}{,}\left\{{a}{=}{0}{,}{b}{=}{0}{,}{c}{=}{c}{,}{u}{}\left({t}\right){=}{\mathrm{_C9}}{}{{t}}^{{6}}{+}{\mathrm{_C8}}{}{{t}}^{{5}}{+}{\mathrm{_C7}}{}{{t}}^{{4}}{+}{\mathrm{_C6}}{}{{t}}^{{3}}{+}{\mathrm{_C5}}{}{{t}}^{{2}}{+}{\mathrm{_C4}}{}{t}{+}{\mathrm{_C3}}\right\}\right]$ (27)
 A nonlinear ODE system
 > ${\mathrm{ode}}_{1}≔\left(\left(-12+12{t}^{2}\right)au\left(t\right)+\left(2+6{t}^{4}-8{t}^{2}\right){a}^{3}-4b{t}^{2}+4b\right)\left(\frac{ⅆ}{ⅆt}u\left(t\right)\right)+\left(-12+12{t}^{2}\right)a\left(\frac{ⅆ}{ⅆt}w\left(t\right)\right)+\left(24-24{t}^{2}\right)av\left(t\right)\left(\frac{ⅆ}{ⅆt}v\left(t\right)\right)+6{a}^{3}{\left(-1+{t}^{2}\right)}^{2}t\left(\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}u\left(t\right)\right)+{a}^{3}{\left(-1+{t}^{2}\right)}^{3}\left(\frac{{ⅆ}^{3}}{ⅆ{t}^{3}}u\left(t\right)\right)=0$
 ${{\mathrm{ode}}}_{{1}}{≔}\left(\left({12}{}{{t}}^{{2}}{-}{12}\right){}{a}{}{u}{}\left({t}\right){+}\left({6}{}{{t}}^{{4}}{-}{8}{}{{t}}^{{2}}{+}{2}\right){}{{a}}^{{3}}{-}{4}{}{b}{}{{t}}^{{2}}{+}{4}{}{b}\right){}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right){+}\left({12}{}{{t}}^{{2}}{-}{12}\right){}{a}{}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{w}{}\left({t}\right)\right){+}\left({-}{24}{}{{t}}^{{2}}{+}{24}\right){}{a}{}{v}{}\left({t}\right){}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{v}{}\left({t}\right)\right){+}{6}{}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}{}{t}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right){+}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{3}}{}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{t}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{u}{}\left({t}\right)\right){=}{0}$ (28)
 > ${\mathrm{ode}}_{2}≔\left(\left(6-6{t}^{2}\right)au\left(t\right)+\left(-2-6{t}^{4}+8{t}^{2}\right){a}^{3}-2b{t}^{2}+2b\right)\left(\frac{ⅆ}{ⅆt}w\left(t\right)\right)-6{a}^{3}{\left(-1+{t}^{2}\right)}^{2}t\left(\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}w\left(t\right)\right)-{a}^{3}{\left(-1+{t}^{2}\right)}^{3}\left(\frac{{ⅆ}^{3}}{ⅆ{t}^{3}}w\left(t\right)\right)=0$
 ${{\mathrm{ode}}}_{{2}}{≔}\left(\left({-}{6}{}{{t}}^{{2}}{+}{6}\right){}{a}{}{u}{}\left({t}\right){+}\left({-}{6}{}{{t}}^{{4}}{+}{8}{}{{t}}^{{2}}{-}{2}\right){}{{a}}^{{3}}{-}{2}{}{b}{}{{t}}^{{2}}{+}{2}{}{b}\right){}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{w}{}\left({t}\right)\right){-}{6}{}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}{}{t}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{w}{}\left({t}\right)\right){-}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{3}}{}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{t}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{w}{}\left({t}\right)\right){=}{0}$ (29)
 > ${\mathrm{ode}}_{3}≔\left(\left(6-6{t}^{2}\right)au\left(t\right)+\left(-2-6{t}^{4}+8{t}^{2}\right){a}^{3}-2b{t}^{2}+2b\right)\left(\frac{ⅆ}{ⅆt}v\left(t\right)\right)-{a}^{3}{\left(-1+{t}^{2}\right)}^{3}\left(\frac{{ⅆ}^{3}}{ⅆ{t}^{3}}v\left(t\right)\right)-6{a}^{3}{\left(-1+{t}^{2}\right)}^{2}t\left(\frac{{ⅆ}^{2}}{ⅆ{t}^{2}}v\left(t\right)\right)=0$
 ${{\mathrm{ode}}}_{{3}}{≔}\left(\left({-}{6}{}{{t}}^{{2}}{+}{6}\right){}{a}{}{u}{}\left({t}\right){+}\left({-}{6}{}{{t}}^{{4}}{+}{8}{}{{t}}^{{2}}{-}{2}\right){}{{a}}^{{3}}{-}{2}{}{b}{}{{t}}^{{2}}{+}{2}{}{b}\right){}\left(\frac{{ⅆ}}{{ⅆ}{t}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{v}{}\left({t}\right)\right){-}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{3}}{}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{t}}^{{3}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{v}{}\left({t}\right)\right){-}{6}{}{{a}}^{{3}}{}{\left({{t}}^{{2}}{-}{1}\right)}^{{2}}{}{t}{}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}\phantom{\rule[-0.0ex]{0.4em}{0.0ex}}{v}{}\left({t}\right)\right){=}{0}$ (30)
 > $\mathrm{sys}≔\left[{\mathrm{ode}}_{1},{\mathrm{ode}}_{2},{\mathrm{ode}}_{3}\right]:$
 Set the infolevel to 1 to see basic information on the computation of a non-singular polynomial solutions for it
 > ${\mathrm{infolevel}}_{\mathrm{polysols}}≔1$
 ${{\mathrm{infolevel}}}_{{\mathrm{polysols}}}{≔}{1}$ (31)
 > ${\mathrm{DEtools}}_{\mathrm{polysols}}\left(\mathrm{sys},\mathrm{parameters}=\left\{a,b\right\},\mathrm{singsol}=\mathrm{false}\right)$
 * Power series solution [1]: {u(t) = t^2*A[1,2]+t*A[1,1]+A[1,0], v(t) = t*A[2,1]+A[2,0], w(t) = t^4*A[3,4]+t^3*A[3,3]+t^2*A[3,2]+t*A[3,1]+A[3,0]}, splitting into cases with respect to the parameters {a, b} * Power series solution [2]: {u(t) = t^2*A[1,2]+t*A[1,1]+A[1,0], v(t) = t^2*A[2,2]+t*A[2,1]+A[2,0], w(t) = t*A[3,1]+A[3,0]}, splitting into cases with respect to the parameters {a, b}
 $\left[\left\{{a}{=}{a}{,}{b}{=}{-}\frac{{{a}}^{{4}}{-}{3}{}{{\mathrm{_C3}}}^{{2}}}{{2}{}{a}}{,}{u}{}\left({t}\right){=}{-}{{a}}^{{2}}{}{{t}}^{{2}}{+}\frac{{{a}}^{{4}}{-}{{\mathrm{_C3}}}^{{2}}}{{2}{}{{a}}^{{2}}}{,}{v}{}\left({t}\right){=}{\mathrm{_C3}}{}{t}{+}\frac{{\mathrm{_C5}}}{{2}{}{\mathrm{_C3}}}{,}{w}{}\left({t}\right){=}{\mathrm{_C5}}{}{t}{+}{\mathrm{_C4}}\right\}\right]$ (32)
 To see the whole set of solutions including the singular ones try the input above without singsol=false.