real*8 function besj0(x) C C used as an interface to the bes0 function written in ndp c C real*8 x real xx,y xx=x call bes0(xx,y) besj0 = y return end real*8 function besj1(x) C C used as an interface to the bes1 function written in ndp c C real*8 x real xx,y xx=x call bes1(xx,y) besj1 = y return end