c----------------------------------------------------------------------
c     Computes definite integral of 
c
c          exp(-x^2)
c
c----------------------------------------------------------------------
c     neq = 1
c----------------------------------------------------------------------

c----------------------------------------------------------------------
c     Implements differential equations.
c     Called by ODEPACK routine LSODA.
c----------------------------------------------------------------------
      subroutine fcn(neq,t,y,yprime)
         implicit    none

         include    'fcn.inc'

         integer     neq
         real*8      t,     y(neq),    yprime(neq)

         yprime(1) = exp(-t**2)

         return
      end 

c----------------------------------------------------------------------
c     Implements Jacobian  (optional) ...
c----------------------------------------------------------------------
      subroutine jac
         implicit    none

         include    'fcn.inc'

         return
      end