In this work we present a numerical approach for finding positive solutions of the type Du = kf(u) for x 2 X, 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.