<p style='text-indent:20px;'>We consider the following semilinear problem with a gradient term in the nonlinearity</p><p style='text-indent:20px;'><disp-formula> <label/> <tex-math id="FE1"> \begin{document}$ \begin{align*} -\Delta u = \lambda \frac{(1+|\nabla u|^q)}{(1-u)^p}\quad\text{in}\quad\Omega,\quad u>0\quad \text{in}\quad \Omega, \quad u = 0\quad\text{on}\quad \partial \Omega. \end{align*} $\end{document} </tex-math></disp-formula></p><p style='text-indent:20px;'>where <inline-formula><tex-math id="M1">\begin{document}$ \lambda,p,q>0 $\end{document}</tex-math></inline-formula> and <inline-formula><tex-math id="M2">\begin{document}$ \Omega $\end{document}</tex-math></inline-formula> be a bounded, smooth domain in <inline-formula><tex-math id="M3">\begin{document}$ {\mathbb R}^N $\end{document}</tex-math></inline-formula>. We prove that when <inline-formula><tex-math id="M4">\begin{document}$ \Omega $\end{document}</tex-math></inline-formula> is a unit ball and <inline-formula><tex-math id="M5">\begin{document}$ p = 1 $\end{document}</tex-math></inline-formula> for <inline-formula><tex-math id="M6">\begin{document}$ q\in (0,q^*(N)) $\end{document}</tex-math></inline-formula> with <inline-formula><tex-math id="M7">\begin{document}$ q^*(N)\in (1,2) $\end{document}</tex-math></inline-formula>, we have infinitely many radial solutions for <inline-formula><tex-math id="M8">\begin{document}$ 2\leq N<2\frac{6-q+2\sqrt{8-2q}}{(2-q)^2}+1 $\end{document}</tex-math></inline-formula> and <inline-formula><tex-math id="M9">\begin{document}$ \lambda = \tilde \lambda $\end{document}</tex-math></inline-formula>. On the other hand, for <inline-formula><tex-math id="M10">\begin{document}$ N>2\frac{6-q+2\sqrt{8-2q}}{(2-q)^2}+1 $\end{document}</tex-math></inline-formula> there exists a unique radial solution for <inline-formula><tex-math id="M11">\begin{document}$ 0<\lambda<\tilde \lambda $\end{document}</tex-math></inline-formula>.</p>