\[ y^{\prime }=y\left ( 1-y\right ) \] Since this is linear in \(y^{\prime }\) we can not use the p-discriminant method. By inspection we see this is separable. Hence the candidate singular solutions are obtain when \(y\left ( 1-y\right ) =0\). This is because this is what we have to divide both sides by to integrate. Therefore \begin {align*} y_{s} & =0\\ y_{s} & =1 \end {align*}
The general solution is found from \begin {align*} \int \frac {dy}{y\left ( 1-y\right ) } & =\int dx\\ -\ln \left ( y-1\right ) +\ln y & =x+c\\ \ln \frac {y}{y-1} & =x+c\\ \frac {y}{y-1} & =c_{1}e^{x}\\ y & =\left ( y-1\right ) c_{1}e^{x}\\ y-yc_{1}e^{x} & =-c_{1}e^{x}\\ y\left ( 1-c_{1}e^{x}\right ) & =-c_{1}e^{x}\\ y & =\frac {c_{1}e^{x}}{c_{1}e^{x}-1}\\ y & =\frac {c_{1}}{c_{1}-e^{-x}}\\ y & =\frac {1}{1-c_{2}e^{-x}} \end {align*}
Now we ask, can the singular solutions \(y_{s}=0,y_{s}=1\) be obtained from the above general solution for any value of \(c_{2}\)? We see when \(c_{2}=0\) then \(y=1\). Also when \(c_{2}=\infty \) then \(y=0\). So these are not really singular solutions. Some call these equilibrium solutions. But these should not be called singular solutions. Mathematica generates these when using the option IncludeSingularSolutions. But Maple does not give these when using the option singsol=all.
The following plot shows these equilibrium solutions with the general solution plotted using diļ¬erent values of \(c\).