1! { dg-do compile }
2!
3! Check whether BESXY functions take scalars and
4! arrays as arguments (PR31760).
5!
6PROGRAM test_erf
7  REAL :: r = 0.0, ra(2) = (/ 0.0, 1.0 /)
8
9  r  = BESJ0(r)
10  r  = BESJ1(r)
11  r  = BESJN(0, r)
12
13  r  = BESY0(r)
14  r  = BESY1(r)
15  r  = BESYN(0, r)
16
17  ra = BESJ0(ra)
18  ra = BESJ1(ra)
19  ra = BESJN(0, ra)
20
21  ra = BESY0(ra)
22  ra = BESY1(ra)
23  ra = BESYN(0, ra)
24
25  r  = BESSEL_J0(r)
26  r  = BESSEL_J1(r)
27  r  = BESSEL_JN(0, r)
28
29  r  = BESSEL_Y0(r)
30  r  = BESSEL_Y1(r)
31  r  = BESSEL_YN(0, r)
32
33  ra = BESSEL_J0(ra)
34  ra = BESSEL_J1(ra)
35  ra = BESSEL_JN(0, ra)
36
37  ra = BESSEL_Y0(ra)
38  ra = BESSEL_Y1(ra)
39  ra = BESSEL_YN(0, ra)
40
41END PROGRAM
42