This ode is an example where \(y\) does not appear explicitly in the ode so not possible to directly solve for \(y\). It
is given here to show possible problems with this method.
\begin{equation} y^{\prime }=\sqrt {1+x+y} \tag {1A}\end{equation}
This ode is squared to first solve for
\(y\)
which gives
\begin{equation} \left ( y^{\prime }\right ) ^{2}=1+x+y \tag {2A}\end{equation}
However, here care is needed. To get back to original ode (1A) then (2A) means
two possible equations
\[ y^{\prime }=\pm \sqrt {1+x+y}\]
Hence the solutions obtained using (2A) can be the solution to one of
these
\begin{align} y^{\prime } & =+\sqrt {1+x+y}\tag {B1}\\ y^{\prime } & =-\sqrt {1+x+y} \tag {B2}\end{align}
Therefore the solution obtained by squaring both sides of (1A), which is done in order to solve for \(y\), must
be checked to see if it satisfies the original ode, else it will be extraneous solution resulting from squaring
both sides of the ode.
Starting from (2A), in normal form (by replacing \(y^{\prime }\) with \(p\)) it becomes
\begin{align} y & =-x-1+p^{2}\tag {1}\\ & =xf+g\nonumber \end{align}
Where \(f=-1,g=-1+p^{2}\). Taking derivative w.r.t. \(x\) gives
\begin{align} p & =f+\left ( xf^{\prime }+g^{\prime }\right ) \frac {dp}{dx}\nonumber \\ p+1 & =2p\frac {dp}{dx} \tag {2}\end{align}
Since \(\frac {\partial \phi }{\partial x}=-1\neq p\) then this is d’Alembert ode. The singular solution is found by setting \(\frac {dp}{dx}=0\) which results in
\(p=-1\). Substituting this in (1) gives the singular solution
\begin{equation} y\left ( x\right ) =-x \tag {3}\end{equation}
But this solution does not satisfy the ode, hence it is
extraneous. The
general solution is found by finding
\(p\) from (2). Since (2) is nonlinear, then it is inverted
which gives
\begin{align*} \frac {p+1}{2p} & =\frac {dp}{dx}\\ \frac {dx}{dp} & =\frac {2p}{p+1}\end{align*}
Which is linear in \(x\). Solving gives
\begin{equation} x=2p-2\ln \left ( p+1\right ) +c_{1} \tag {4}\end{equation}
Instead of inverting this to find
\(p\) in terms of
\(x\),
\(p\) is found from (1) which
gives
\begin{align*} y+x+1 & =p^{2}\\ p & =\pm \sqrt {y+x+1}\end{align*}
Substituting these solutions in (4) gives implicit solutions as
\begin{align*} x & =2\sqrt {y+x+1}-2\ln \left ( 1+\sqrt {y+x+1}\right ) +c_{1}\\ x & =-2\sqrt {y+x+1}-2\ln \left ( 1-\sqrt {y+x+1}\right ) +c_{1}\end{align*}
But only the first one above satisfies the ode. The second is extraneous. Therefore the final solution is
\[ x=2\sqrt {y+x+1}-2\ln \left ( 1+\sqrt {y+x+1}\right ) +c_{1}\]
And
no singular solutions exist. If instead of doing the above,
\(p\) was found from (4) using inversion, then it will
be
\[ p=-\operatorname *{LambertW}\left ( -c_{1}e^{\frac {-x}{2}-1}\right ) -1 \]
Substituting this in (1) gives
\[ y=-x-1+\left ( -\operatorname *{LambertW}\left ( -c_{1}e^{\frac {-x}{2}-1}\right ) -1\right ) ^{2}\]
But this general solution does not satisfy the original ode. In general, it
is best to avoid squaring both side of the ode in order to solve for
\(y\) as this can generate extraneous
solutions. Only use this method if the original ode is already given in the form where
\(y\) shows
explicitly.