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

einstein% pwd; ls
/usr2/people/phy329/fd/ex2
Makefile     rexbvp1d.f

einstein% make
        f77 -g -c rexbvp1d.f
        f77 -g -L/usr/local/lib rexbvp1d.o -lp329f -llapack -lblas -o rexbvp1d

einstein% rexbvp1d
 usage: rexbvp1d <level> [<option>]

        Specify option .ne. 0 for output
        of error instead of solution

einstein% rexbvp1d 5
   0.0000000000000000E+00  9.6589403142388603E-15
   6.2500000000000000E-02  0.3826738794434802    
   0.1250000000000000      0.7070891296890379    
   0.1875000000000000      0.9238564697183627    
   0.2500000000000000      0.9999750370128206    
   0.3125000000000000      0.9238564697183612    
   0.3750000000000000      0.7070891296890351    
   0.4375000000000000      0.3826738794434761    
   0.5000000000000000      4.1818400594214227E-15
   0.5625000000000000     -0.3826738794434678    
   0.6250000000000000     -0.7070891296890267    
   0.6875000000000000     -0.9238564697183529    
   0.7500000000000000     -0.9999750370128120    
   0.8125000000000000     -0.9238564697183538    
   0.8750000000000000     -0.7070891296890289    
   0.9375000000000000     -0.3826738794434710    
    1.000000000000000      0.0000000000000000E+00
 rmserr =   1.7124467931210358E-05

#############################################################
# Convergence test: RMS error goes down by a factor of 
# close to 16 at each level indicating O(h^4) convergence.
# Overall results very similar to those from 'bvp1d4'
#############################################################
einstein% foreach level (4 5 6 7 8 9 10)
foreach? rexbvp1d $level > /dev/null
foreach? end
 rmserr =   2.7251125733154054E-04
 rmserr =   1.7124467931210358E-05
 rmserr =   1.0801454943036140E-06
 rmserr =   6.7928843322335184E-08
 rmserr =   4.2604421714816242E-09
 rmserr =   2.6653456044743611E-10
 rmserr =   1.7013500115372792E-11

#############################################################
# Generate some sample output for plotting
#############################################################
einstein% rexbvp1d 4 > out4
 rmserr =   2.7251125733154054E-04