Fortran Wiki
bessel_j1