#############################################################
# Building 'tlsoda' and sample output on the SGIs
#############################################################

einstein% pwd; ls
/usr2/people/phy329/ode/ex1
Makefile   tlsoda.f

#############################################################
#  NOTE: 'lsoda' uses routines from a predecessor of LAPACK
#  called LINPACK as well as BLAS routines.  All three of the
#  libraries: 'odepack', 'linpack', and 'blas' must be 
#  specified at load time.
#############################################################

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

einstein% tlsoda
 usage: tlsoda <y_i> <dy/dx_i>

#############################################################
# Generate a uniform-mesh of x output-values using 'dvmesh'
# and pipe it into tlsoda with initial conditions 
# u(0) = 0, u'(0) = 1  --> u(x) = sin(x)
#############################################################

einstein% dvmesh 0.0 1.0 11 | tlsoda 0.0 1.0
   0.0000000000000000E+00  0.0000000000000000E+00
   0.1000000000000000      9.9833904829367021E-02
   0.2000000000000000      0.1986703509344916    
   0.3000000000000000      0.2955222777237400    
   0.4000000000000000      0.3894213005717221    
   0.5000000000000000      0.4794291455964967    
   0.6000000000000000      0.5646464181480721    
   0.7000000000000000      0.6442215920888248    
   0.7999999999999999      0.7173595216370443    
   0.9000000000000000      0.7833295096871006    
    1.000000000000000      0.8414729721292930

#############################################################
# Use 'nf' to compute error in solution.  Note that, at least
# at the initial time, the error in the solution is roughly
# within the specified tolerance of 1.0d-6.
#############################################################

einstein% dvmesh 0.0 1.0 11 | tlsoda 0.0 1.0 | nf _1 'sin(_1) - _2'
0.0000000000000000E+00   0
0.1000000000000000   -4.88182538865845e-07
0.2000000000000000   -1.02013943034662e-06
0.3000000000000000   -2.07106240041055e-06
0.4000000000000000   -2.95826307156188e-06
0.5000000000000000   -3.606992293681e-06
0.6000000000000000   -3.94475303677755e-06
0.7000000000000000   -3.90485113377359e-06
0.7999999999999999   -3.43073752151479e-06
0.9000000000000000   -2.60005961716381e-06
1.000000000000000   -1.98732139644608e-06