In this paper, we show the existence of an S-shaped connected component in the set of radial positive solutions of boundary value problem$$\begin{array}{}
\displaystyle \left\{\,\begin{array}{}
-\text{ div}\big(\phi_N(\nabla y)\big)=\lambda a(|x|)f(y)\, \, \, \, \, \text{in}\, \, \mathcal{A},\\\frac{\partial y}{\partial \nu}=0\, \, \, \,\, \text{ on }\, \, {\it\Gamma}_1,\qquad y=0\, \, \, \, \text{ on}\, \, {\it\Gamma}_2,\\
\end{array}
\right.
\end{array} $$where R2 β (0, β) and R1 β (0, R2) is a given constant, π = {x β βN : R1 < β£xβ£ < R2}, Ξ1 = {x β βN : β£xβ£ = R1}, Ξ2 = {x β βN : β£xβ£ = R2}, $\begin{array}{}
\phi_N(s)=\frac{s}{\sqrt{1-|s|^2}},
\end{array} $s β βN, Ξ» is a positive parameter, a β C[R1, R2], f β C[0, β), $\begin{array}{}
\frac{\partial y}{\partial \nu}
\end{array} $ denotes the outward normal derivative of y and β£β
β£ denotes the Euclidean norm in βN. The proof of main result is based upon bifurcation techniques.