############################################################
# Building 'integral' and sample usage on SGIs
############################################################
einstein% ls  
Makefile     fcn.f        fcn.inc      integral.f   makefile

einstein% make
        make -f Makefile
        f77 -g -n32 -c integral.f
        f77 -g -n32 -c fcn.f
        f77 -g -n32 -L/usr/localn32/lib -n32 integral.o fcn.o \
            -lp329f -lodepack -llinpack  -lblas -o integral

einstein% integral
 integral: Insufficient initial values  specified on command line.
 usage: integral <iv 1> <iv 2> ... <iv neq> [<tol>]
  
        Output times read from standard input

############################################################
# Use a fairly stringent error tolerance of 1.0d-12
############################################################
einstein% integral 0.0 1.0d-12
-5.0
5.0
^D
   -5.000000000000000      0.0000000000000000E+00
    5.000000000000000       1.772453850933237    

############################################################
# The value 1.772 4538 5093 3237 is the estimate of the 
# integral.  The exact value is 1.772 4538 5090 5516
############################################################