\(y=x\left ( y^{\prime }\right ) ^{2}\) is put in normal form (by replacing \(y^{\prime }\) with \(p\)) and solving for \(y\) gives\begin {align} y & =xp^{2}\tag {1}\\ & =xf\left ( p\right ) \nonumber \end {align}
This is the case when \(f\left ( p\right ) =p^{2}\) and \(g\left ( p\right ) =0\). Writing \(f\equiv f\left ( p\right ) \) and \(g\equiv g\left ( p\right ) \) to make notation simpler but remembering that \(f\) is function of \(p\left ( x\right ) \) which in turn is function of \(x\). Same for \(g\left ( p\right ) \).\[ y=xf \] Taking derivative of the above w.r.t. \(x\) gives\begin {align*} p & =\frac {d}{dx}\left ( xf\right ) \\ p & =f+xf^{\prime }\frac {dp}{dx}\\ p-f & =xf^{\prime }\frac {dp}{dx} \end {align*}
Since \(f=p^{2}\) then the above becomes\begin {equation} p-p^{2}=2xp\frac {dp}{dx} \tag {2} \end {equation} The singular solution is given when \(\frac {dp}{dx}=0\) or \(p-p^{2}=0\). This gives \(p=0\) or \(p=1\). Substituting these values of \(p\) in (1) gives singular solutions\begin {align} y_{s1} & =0\tag {3}\\ y_{s2} & =x \tag {4} \end {align}
General solution is found when \(\frac {dp}{dx}\neq 0\,\). Eq(2) is a first order ode in \(p\). Now we could either solve ode (2) directly as it is for \(p\left ( x\right ) \), or do an inversion and solve for \(x\left ( p\right ) \). We always try the first option first. Since (2) is separable as is, no need to do an inversion. Eq (2) is separable. The solution is\begin {align*} p_{1} & =0\\ p_{2} & =1+\frac {c_{1}}{\sqrt {x}} \end {align*}
For each \(p\), there is a general solution. Substituting each of the above in (1) gives\begin {align*} y_{1}\left ( x\right ) & =0\\ y_{2}\left ( x\right ) & =x\left ( 1+\frac {c_{1}}{\sqrt {x}}\right ) ^{2} \end {align*}
Hence the final solutions are\begin {align*} y & =x\qquad \text {(singular)}\\ y & =0\\ y & =x\left ( 1+\frac {c_{1}}{\sqrt {x}}\right ) ^{2} \end {align*}
But \(y=x\) can be obtained from the general solution when \(c_{1}=0\). Hence it is removed. Therefore the final solutions are\begin {align} y & =0\tag {6}\\ y & =x\left ( 1+\frac {c_{1}}{\sqrt {x}}\right ) ^{2} \tag {7} \end {align}
What will happen if we had done an inversion to \(x\left ( p\right ) \)? Let us find out. ode(5) now becomes\begin {align*} \frac {p-p^{2}}{p}\frac {dx}{dp} & =2x\\ \frac {dx}{2x} & =\frac {p}{p-p^{2}}dp \end {align*}
This is also separable in \(x\). Solving this for \(x\) gives\[ x=\frac {c_{1}}{\left ( p-1\right ) ^{2}}\] Solving for \(p\) from the above gives\begin {align*} p_{1} & =\frac {x+\sqrt {xc_{1}}}{x}\\ p_{2} & =\frac {x-\sqrt {xc_{1}}}{x} \end {align*}
Substituting each of the above in (1) gives\begin {align*} y_{1} & =x\left ( \frac {x+\sqrt {xc_{1}}}{x}\right ) ^{2}\\ & =\frac {\left ( x+\sqrt {xc_{1}}\right ) }{x}^{2}\\ y_{2} & =x\left ( \frac {x-\sqrt {xc_{1}}}{x}\right ) ^{2}\\ & =\frac {\left ( x-\sqrt {xc_{1}}\right ) }{x}^{2} \end {align*}
Now we see that singular solution \(y=x\) can be obtained from the above general solutions from \(c_{1}=0\). But \(y=0\) can not. Hence the final solutions are\begin {align} y & =0\qquad \text {(singular)}\tag {8}\\ y & =\frac {\left ( x+\sqrt {xc_{1}}\right ) }{x}^{2}\tag {9}\\ y & =\frac {\left ( x-\sqrt {xc_{1}}\right ) }{x}^{2} \tag {10} \end {align}
All solutions (6,7,8,9,10) are correct and verified. Maple gives the solutions given in (8,9,10) and not those in (6,7).
Another method to find the singular solutions if it exists is called the p-discriminant. This is used only for first order ode with nonlinear in \(y^{\prime }\). We set up the following two equations\begin {align*} F\left ( x,y,y^{\prime }\right ) & =0\\ \frac {\partial F\left ( x,y,y^{\prime }\right ) }{\partial y^{\prime }} & =0 \end {align*}
We eliminate \(y^{\prime }\) and obtain \(G\left ( x,y\right ) =0\) equation. This is the singular solution. But we still have to check if it satisfies the ode and also if it is true singular solution curve. More on this later. Let us now just find the singular solution found above but using the p-discriminant method. The above two equations are\begin {align*} y-x\left ( y^{\prime }\right ) ^{2} & =0\\ -2xy^{\prime } & =0 \end {align*}
Second equation gives \(y^{\prime }=0\). Hence the first equation now gives the singular solution as \[ y=0 \] Which is the same obtained earlier.