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