Splitting of any initial field distribution in GRIN media

We show the splitting effect of an electromagnetic field assuming non paraxial propagation in a quadratic GRIN medium. The field distribution at the splitting distance is given by fractional Fourier transforms of the initial field and of a reflected version of it. We apply this result to an initial field given by a Bessel function and show that it splits into two generalized Bessel functions.