#!/bin/sh -x
#-----------------------------------------------------------
# Wave: script which runs 'wave' and then graphs output 
# using gnuplot
#-----------------------------------------------------------
test -f wave || make

# Set the command-line parameters for wave
xlevel=8
olevel=6
tfinal=1.8
dtout=0.02
tol=1.0e-5

# Warn the user re the run time 
echo "This may take a couple of minutes or so.  Please be patient."

# Generate the solution
wave $xlevel $olevel $tfinal $dtout $tol > out8

# Plot the solution
gnuplot<<END
set terminal postscript portrait
set title "Method of lines (LSODA) solution of wave \
equation\nTime-symmetric initial data"
set output "out8.ps"
set ticslevel 0.1
set parametric
set hidden
splot "out8" title "u(x,t)" with lines
quit
END
ls -lt *ps