#############################################################
# 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.8818253886584451e-07
0.2000000000000000   -1.0201394303466227e-06
0.3000000000000000   -2.0710624004105505e-06
0.4000000000000000   -2.9582630715618841e-06
0.5000000000000000   -3.6069922936809995e-06
0.6000000000000000   -3.9447530367775485e-06
0.7000000000000000   -3.9048511337735903e-06
0.7999999999999999   -3.4307375215147928e-06
0.9000000000000000   -2.600059617163808e-06
1.000000000000000   -1.9873213964460845e-06