In this work we present a numerical approach for finding positive solutions with Dirichlet boundary condition, where f is a superlinear function of u. We will show in which range of k, this problem achieves multiple numerical solutions and what is the behavior of the branches of solutions.
