＜BESSEL_J1 Example Program＞

Program bessel_j1_example
  Implicit None
  Integer i
  Real :: x(8) = [ (i*0.5,i=1,8) ]
  Print 1, x
  Print 1, Bessel_J1(x)
1 Format(999F9.3)
End Program

---

■ Execution Results

    0.500    1.000    1.500    2.000    2.500    3.000    3.500    4.000
    0.242    0.440    0.558    0.577    0.497    0.339    0.137   -0.066