＜BESSEL_Y1 Example Program＞

Program bessel_y1_example
  Implicit None
  Integer i
  Real :: x(8) = [ (i*0.5,i=1,8) ]
  Print 1, x
  Print 1, Bessel_Y1(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
   -1.471   -0.781   -0.412   -0.107    0.146    0.325    0.410    0.398