For the case n > n0 + 1, I can't comprehend how Iterating \begin{equation} y[n] = \lambda y[n-1] \end{equation} yields \begin{equation} y[n] = \left(\lambda\right)^{n-n_0}\frac{b \;}{\lambda} \;\;\; n > n_0. \end{equation}
For the entire description, please visit this page.