In this paper we will investigate the existence of multiple
solutions for the problem $-\Delta_pu + g(x,u) =\lambda_1h(x)|u|^{p-2}u,\mathrm{in}\Omega,u\in H_0^{1,p}(\Omega)(P)$ where $\Delta_{p}u =\mathrm{div}(|\nabla u|^{p-2}\nabla u)$ is the $p$ -Laplacian operator, $\Omega\subseteq\mathbb{R}^N$ is a bounded domain with smooth boundary, $h$ and $g$ are bounded functions, $N\geq 1$ and $1 < p < \infty$ . Using the Mountain Pass Theorem and the Ekeland Variational Principle, we will show the existence of at least three solutions for (P).