A generalization of di-Bessel function of Exton (Q788873)

If the exponential generating function of the Bessel function is generalized to the product of two hypergeometric functions \({}_ 0F_{n-1}\), in which all the parameters are put equal to unity, the function generated is \[ A_{\nu,n}(x)=\frac{(x/2)^{\nu}}{\{\Gamma(\nu +1)\}^ n}_ 0F_{2n-1}(-;\nu +1,...,\nu +1,1,...,1;-\frac{x^ 2}{4}). \] If \(n=1\), we have the Bessel function and for \(n=2\), the so- called di-Bessel function, first introduced by the reviewer [Indian. J. Pure Appl. Math. 11, 856-862 (1980; Zbl 0443.33006)] is obtained. This paper gives many interesting results which effect considerable generalization and extension of the existing work on the di-Bessel function. This includes orthogonality relations of Fourier-Bessel type arising from a Sturm-Liouville system of order 2n.