c===========================================================
c     Implements ODE for computation of definite integral of 
c
c          exp(-x^2)
c===========================================================
      subroutine fcn(neq,x,y,yprime)
         implicit    none

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

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

         return
      end 

c===========================================================
c     Dummy Jacobian routine.
c===========================================================
      subroutine jac
         implicit    none

         return
      end
