Find the fixed points of the implicit Euler scheme

\begin{equation} y_{n+1}-y_{n}= hf(t_{n+1},y_{n+1})

\end{equation}

when applied to the differential equation $y'=y(1-y)$ and investigate their stability?

=>

implicit Euler scheme

\begin{equation} y_{n+1}-y_{n}= hf(t_{n+1},y_{n+1})

\end{equation}

$y'=y(1-y)$

\begin{equation} y_{n+1}=y_{n}+hy_{n+1}(1-y_{n+1})

\end{equation}

\begin{equation} y_{n+1}=y_{n}+hy_{n+1}-hy^2_{n+1}

\end{equation}

For fixed points

$y_{n+1}=y_{n}$

\begin{equation} y_{n}=y_{n}+hy_{n}-hy^2_{n}

\end{equation}

$y_{n}=0$ or $1$

I got problem with stability but this is what I have done

$y_{n}= \alpha +\epsilon^n$, $y_{n+1}= \alpha +\epsilon^{n+1}$,

\begin{equation} \alpha +\epsilon^{n+1}= \alpha +\epsilon^n + h (\alpha +\epsilon^{n+1})(1-\alpha -\epsilon^{n+1}) \end{equation}

\begin{equation} \epsilon^{n+1}= \epsilon^n + h (\alpha +\epsilon^{n+1})(1-\alpha -\epsilon^{n+1}) \end{equation}

When $y_{n}=0=\alpha$

\begin{equation} \epsilon^{n+1}= \epsilon^n + h \epsilon^{n+1}(1-\epsilon^{n+1}) \end{equation}

I don't what to say or do after that to determine the stability.

When $y_{n}=1=\alpha$

\begin{equation} \epsilon^{n+1}= \epsilon^n - h \epsilon^{n+1}(1+\epsilon^{n+1}) \end{equation}

same again what can say about with my answer to investigate the stability.