We present the Rothe method for the McKendrick-von Foerster equation with initial and boundary conditions. This method is well known as an abstract Euler scheme in extensive literature, e.g. K. Rektorys, The Method of Discretization in Time and Partial Differential Equations, Reidel, Dordrecht, 1982. Various Banach spaces are exploited, the most popular being the space of bounded and continuous functions. We prove the boundedness of approximate solutions and stability of the Rothe method in $L^\infty $ and $L^1$ norms. Proofs of these results are based on comparison inequalities. Our theory is illustrated by numerical experiments. Our research is motivated by certain models of mathematical biology.