c===========================================================
c     Demonstrates sample use of SLATEC library to evaluate
c     a Bessel function and an Airy function.
c
c     % maple
c     > Digits := 16:
c     > evalf(BesselJ(0,1));
c       .7651976865579666
c     > evalf(AiryAi(1));
c       .1352924163128814
c===========================================================
      program   tslatec

      implicit  none

      real*8    dbesj0,  dai

c-----------------------------------------------------------
c     Compute and output Bessel function of first kind of 
c     order 0.
c-----------------------------------------------------------
      write(0,1000) dbesj0(1.0d0)
1000  format(' dbesj0(1.0d0) = ',1p,E24.16)
c-----------------------------------------------------------
c     Compute and output Airy function.
c-----------------------------------------------------------
      write(0,1100) dai(1.0d0)
1100  format(' dai(1.0d0)    = ',1p,E24.16)

      stop

      end
