A numerical algorithm for finding solution of multiparameter semipositone Dirichlet problems

In this paper we investigate numerically positive solution of the equation -Δu=λf(x,u)+μg(x,u) with Dirichlet boundary condition in a boundary domain Ω for λ>0, μ∈R. We will show that in which range of λ, this problem achieves a numerical solution.