Starlikeness of the generalized Bessel function

Abstract: For a fixed $a \in {1, 2, 3, \ldots},$ the radius of starlikeness of positive order is obtained for each of the normalized analytic functions \begin{align*} \mathtt{f}{a, \nu}(z)&:= \bigg(2^{a \nu-a+1} a^{-\frac{a(a\nu-a+1)}{2}} \Gamma(a \nu+1) {}_a\mathtt{B}{2a-1, a \nu-a+1, 1}(a^{a/2} z)\bigg)^{\tfrac{1}{a \nu-a+1}},\ \mathtt{g}{a, \nu}(z)&:= 2^{a \nu-a+1} a^{-\frac{a}{2}(a\nu-a+1)} \Gamma(a \nu+1) z^{a-a\nu} {}_a\mathtt{B}{2a-1, a \nu-a+1, 1}(a^{a/2} z),\ \mathtt{h}{a, \nu}(z)&:= 2^{a \nu-a+1} a^{-\frac{a}{2}(a\nu-a+1)} \Gamma(a \nu+1) z^{\frac{1}{2}(1+a-a\nu)} {}_a\mathtt{B}{2a-1, a \nu-a+1, 1}(a^{a/2} \sqrt{z}) \end{align*} in the unit disk, where ${}a\mathtt{B}{b, p, c}$ is the generalized Bessel function \begin{align*} {}a\mathtt{B}{b, p, c}(z):= \sum_{k=0}^\infty \frac{(-c)^k}{k! \; \mathrm{\Gamma}{\left( a k +p+\frac{b+1}{2}\right)} } \left(\frac{z}{2}\right)^{2k+p}. \end{align*} The best range on $\nu$ is also obtained for a fixed $a$ to ensure the functions $\mathtt{f}{a, \nu}$ and $\mathtt{g}{a, \nu}$ are starlike of positive order in the unit disk. When $a=1,$ the results obtained reduced to earlier known results.