program f08_BESSEL
    implicit none
    integer::i,n
    real(8)::a,b,x

    n=100
    a=0.0D0
    b=10.0D0
    do i=0,n
        x=a+(b-a)/dble(n)*dble(i)
        write(6,'(5f10.5)') x,bessel_jn(0,x),bessel_jn(1,x),bessel_jn(2,x),bessel_jn(3,x)
    end do

end program f08_BESSEL