We prove existence results of two solutions of the problem [ begin{cases} L(u)+u^{m-1}=lambda u^{p-1} & text{ in $Omega$}, quad u>0 &text{ in $Omega$}, quad u=0 & text{ on $partial Omega$}, end{cases} ] where $L(v)=-{rm div}(M(x) abla v)$ is a linear operator, $pin (2,2^{*}]$ and $lambda$ and $ m$ sufficiently large. Then their asymptotical limit as $mto +infty$ is investigated showing different behaviors.