The Fourier-Legendre series of Bessel functions of the first kind and the summed series involving _2F_3 hypergeometric functions that arise from them

The Bessel function of the first kind J_N(kx) is expanded in a Fourier-Legendre series, as is the modified Bessel functions of the first kind I_N(kx). The purpose of these expansions in Legendre polynomials was not an attempt to rival established numerical methods for calculating Bessel functions, but to provide a form for J_N(kx) useful for analytical work in the area of strong laser fields, where analytical integration over scattering angles is essential. Despite their primary purpose, we can easily truncate the series at 21 terms to provide 33-digit accuracy that matches IEEE extended precision in some compilers. The analytical theme is furthered by showing that infinite series of like-powered contributors (involving _2F_3 hypergeometric functions) extracted from the Fourier-Legendre series may be summed, having values that are inverse powers of the eight primes 1/(2^i3^j5^k7^l11^m13^n17^o19^p) multiplying powers of the coefficient k.