Factorization method for solving differential equations

Description of the method

The factorization method is a method that can be used to solve certain kinds of differential equations. The idea behind the method is to start with a differential equation:

$F(x,y,y',\dots ,y^{(k)})=0$

and try to factor the expression $F$ as a product of two expressions, say $G$ and $H$. The differential equation now becomes:

$G(x,y,y',\dots ,y^{(k)})H(x,y,y',\dots ,y^{(k)})=0$

We now consider the two equations:

$G(x,y,y',\dots ,y^{(k)})=0$

$H(x,y,y',\dots ,y^{(k)})=0$

It is clear that any functional (respectively, relational) solution for either of the equations is a functional (respectively, relational) solution for $F$. Conversely, any functional (respectively, relational) solution for $F$ must, at every point in the domain, satisfy one of the differential equations. However, we need to worry about the possibility of a functional solution for $F$ that is a solution for $G$ on part of its domain and a solution for $H$ on part of its domain. It is quite rare for these "mixed" solutions to occur, but it does make sense to check for them.

The general solution for $F$ can thus be expressed as:

(General solution for $G$) $\bigcup$ (General solution for $H$) $\bigcup$ (Any solution for $F$ that is a solution for $G$ on part of the domain and a solution for $H$ on part of the domain)

Applicability

The factorization method is limited to situations where the expression can be factorized. Further, even if factorization is possible, there is no guarantee that the factor differential equations can be solved. Thus, the method has limited applicability.

However, there are some situations where this method can be used reliably:

• Using the quadratic formula, it can help convert second-degree differential equations to first-degree differential equations. In other words, if the differential equation is quadratic in the highest order derivative, we could factor it as a product of differential equations that are linear in the highest order derivative using the quadratic formula. The quadratic formula involves radicals, so this may come at the expense of introducing radicals into the differential equation.
• In particular, any first-order second-degree autonomous differential equation can be factored to give two first-order first-degree autonomous differential equations, which we have a general method for solving.

Counting degrees of freedom

The total number of degrees of freedom (i.e., freely varying parameters in the general solution) for $F$ is generally taken to be the maximum of the number of degrees of freedom for each of the solution families. In particular, assuming there are no mixed solutions, the number of degrees of freedom for $F$ is the maximum of the number of degrees of freedom for $G$ and the number of degrees of freedom for $H$.

The reason we take the maximum instead of the sum is that the freely varying parameters for the separate solution families are not interacting with each other. Geometrically, this is related to the idea that when we take the union of two one-dimensional things (lines or curves) we still get a one-dimensional thing (line or curve).

Examples

Example without mixed solutions

Consider the differential equation:

$(y')^{2}-3yy'+2y^{2}=0$

This can be factored to give:

$(y'-y)(y'-2y)=0$

Thus, the differential equations we need to solve individually are:

$y'=y,\qquad y'=2y$

The solutions are respectively:

$y=C_{1}e^{x},\qquad y=C_{2}e^{2x},\qquad C_{1},C_{2}\in \mathbb {R}$

The general solution includes both solution families. Let's now consider the possibility of a mixed solution, i.e., a function $y=f(x)$ such that $y'=y$ for part of the domain and $y'=2y$ for the rest of the domain. At a transition point, we must have $y'=y=2y$ forcing $y=0$. But in either solution, setting $y=0$ forces the solution identically to be the zero function. This actually can be realized completely as a solution in either family, hence the only mixed solution actually lies as a pure solution in both families. The general solution is thus a union of two solution families:

• $y=Ce^{x},C\in \mathbb {R}$
• $y=Ce^{2x},C\in \mathbb {R}$

Example with mixed solutions

Consider the differential equation:

$(y'-2x)(y'-3x^{2})=0$

We want to find all possible functional solutions for this on all of $\mathbb {R}$.

The solution families corresponding to the individual factors are:

$y=x^{2}+C_{1},\qquad y=x^{3}+C_{2},\qquad C_{1},C_{2}\in \mathbb {R}$

We now need to consider the possibility of mixed solutions. A mixed solution would be a function that is of the form $x^{2}+C_{1}$ on part of the domain and $x^{3}+C_{2}$ on the rest of the domain, possibly even transitioning multiple times between the two solution forms. Suppose $x=x_{0}$ is a transition point. Then, the $y$-values and $y'$-values for both functions must agree at the point, so:

$x_{0}^{2}+C_{1}=x_{0}^{3}+C_{2},\qquad 2x_{0}=3x_{0}^{2}$

Solving the second equation, we get that $x_{0}=0$ or $x_{0}=2/3$.

• Case of a single transition point $x_{0}=0$: In this case, we get $C_{1}=C_{2}$, so the function description for $y$ could be either of these, with one solution of each type for each $C_{1}\in \mathbb {R}$:

$y=\left\lbrace {\begin{array}{rl}x^{2}+C_{1},&x\leq 0\\x^{3}+C_{1},&x>0\\\end{array}}\right.$

or a transition in the opposite sense:

$y=\left\lbrace {\begin{array}{rl}x^{3}+C_{1},&x\leq 0\\x^{2}+C_{1},&x>0\\\end{array}}\right.$

• Case of a single transition point $x_{0}=2/3$: In this case, $C_{2}=C_{1}+4/27$, so the function description could be either of these, with one solution of each type for each $C_{1}\in \mathbb {R}$:

$y=\left\lbrace {\begin{array}{rl}x^{2}+C_{1},&x\leq 2/3\\x^{3}+C_{1}+4/27,&x>2/3\\\end{array}}\right.$

or a transition in the opposite sense:

$y=\left\lbrace {\begin{array}{rl}x^{3}+C_{1}+4/27,&x\leq 2/3\\x^{2}+C_{1},&x>2/3\\\end{array}}\right.$

• Case of two transition points $x_{0}=0$ and $x_{0}=2/3$: This situation is trickier, The first type of transition is as follows:

$y=\left\lbrace {\begin{array}{rl}x^{2}+C_{1},&x\leq 0\\x^{3}+C_{2},&0<x\leq 2/3\\x^{2}+C_{3},&x>2/3\\\end{array}}\right.$

We get the relations $C_{1}=C_{2}=C_{3}+4/27$, so this becomes:

$y=\left\lbrace {\begin{array}{rl}x^{2}+C_{3}+4/27,&x\leq 0\\x^{3}+C_{3}+4/27,&0<x\leq 2/3\\x^{2}+C_{3},&x>2/3\\\end{array}}\right.$

The other possibility is:

$y=\left\lbrace {\begin{array}{rl}x^{3}+C_{2},&x\leq 0\\x^{2}+C_{1},&0<x\leq 2/3\\x^{3}+C_{3},&x>2/3\\\end{array}}\right.$

In this case, we get $C_{2}=C_{1}$ and $C_{3}=C_{1}+(4/27)$, so plugging in, we get:

$y=\left\lbrace {\begin{array}{rl}x^{3}+C_{1},&x\leq 0\\x^{2}+C_{1},&0<x\leq 2/3\\x^{3}+C_{1}+4/27,&x>2/3\\\end{array}}\right.$

Thus, we overall have eight solution families, each with one freely varying real parameter (which we will now just call $C$ instead of $C_{1},C_{2}C_{3}$):

• $y=x^{2}+C,C\in \mathbb {R}$
• $y=x^{3}+C,C\in \mathbb {R}$
• $y=\left\lbrace {\begin{array}{rl}x^{2}+C,&x\leq 0\\x^{3}+C,&x>0\\\end{array}}\right.,C\in \mathbb {R}$
• $y=\left\lbrace {\begin{array}{rl}x^{3}+C,&x\leq 0\\x^{2}+C,&x>0\\\end{array}}\right.,C\in \mathbb {R}$
• $y=\left\lbrace {\begin{array}{rl}x^{2}+C,&x\leq 2/3\\x^{3}+C+4/27,&x>2/3\\\end{array}}\right.,C\in \mathbb {R}$
• $y=\left\lbrace {\begin{array}{rl}x^{3}+C+4/27,&x\leq 2/3\\x^{2}+C,&x>2/3\\\end{array}}\right.,C\in \mathbb {R}$
• $y=\left\lbrace {\begin{array}{rl}x^{2}+C+4/27,&x\leq 0\\x^{3}+C+4/27,&0<x\leq 2/3\\x^{2}+C,&x>2/3\\\end{array}}\right.,C\in \mathbb {R}$
• $y=\left\lbrace {\begin{array}{rl}x^{3}+C,&x\leq 0\\x^{2}+C,&0<x\leq 2/3\\x^{3}+C+4/27,&x>2/3\\\end{array}}\right.,C\in \mathbb {R}$

Factoring using quadratic formula

Consider the differential equation with independent variable $x$ and dependent variable $y$:

$(y')^{2}+y^{2}=1$

This is a first-order second-degree autonomous differential equation.

Rewrite as:

$(y')^{2}-(1-y^{2})=0$

This does not factor in purely polynomial terms. However, if we allow the use of radicals, we can rewrite it as:

$(y'-{\sqrt {1-y^{2}}})(y'+{\sqrt {1-y^{2}}})=0$

The two factor differential equations are:

$y'={\sqrt {1-y^{2}}},\qquad y'=-{\sqrt {1-y^{2}}}$

Each of them is first-order first-degree autonomous, and in particular separable. The solutions are:

$\arcsin y=x+C_{1},\qquad \arcsin y=-x+C_{2},\qquad C_{1},C_{2}\in \mathbb {R}$

Also, both equations have stationary solutions $y=1$ and $y=-1$.

The following is the general solution and includes all the nowhere stationary functional solutions.

$y=\sin(x+C),C\in \mathbb {R}$

If we want to include somewhere stationary solutions, then the space of solutions becomes complicated. Basically, we could have a function that is an increasing sinusoidal function for part of the domain, then stationary at 1 for a while, then a decreasing sinusoidal function till it reaches -1, then stationary at -1 for a while, and so on.