＜BESSEL_YN Example Program 1＞

Program bessel_yn_example_1
  Implicit None
  Integer i
  Real :: x(8) = [ (i*0.5,i=1,8) ]
  Print 1, x
  ! Elemental over X.
  Print 1, Bessel_YN(3,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
  -42.059   -5.822   -2.074   -1.128   -0.756   -0.539   -0.358   -0.182

＜BESSEL_YN Example Program 2＞

Program bessel_yn_example_2
  Implicit None
  Integer i
  Real :: x = 9.5
  Integer :: n(7) = [ 1,3,5,7,11,13,17 ]
  Print 1, (n(i),x,i=1,Size(n))
1 Format(999(2X,'J',I2.2,'(',F3.1,')',:))
  ! Elemental over N.
  Print 2, Bessel_YN(n,x)
2 Format(999F10.5)
End Program

---

■ Execution Results

  J01(9.5)  J03(9.5)  J05(9.5)  J07(9.5)  J11(9.5)  J13(9.5)  J17(9.5)
   0.20318  -0.25726   0.22859   0.11835  -0.64735  -2.01846 -92.43208

＜BESSEL_YN Example Program 3＞

Program bessel_yn_example_3
  Implicit None
  Integer i
  Real :: x = 2.5
  Print 1, (i,x,i=1,7)
1 Format(999(3X,'J',I1,'(',F3.1,')',:))
  ! Transformational.
  Print 2, Bessel_YN(1,7,x)
2 Format(999F10.5)
End Program

---

■ Execution Results

   J1(2.5)   J2(2.5)   J3(2.5)   J4(2.5)   J5(2.5)   J6(2.5)   J7(2.5)
   0.14592  -0.38134  -0.75606  -1.43320  -3.83018 -13.88751 -62.82986