#!/bin/csh -f
rm schrodinger
mv wavefunctions/wave* old_data/
mv *.txt old_data/
cc -o schrodinger schrod.c /usr/local/lib/librecipes_c.a
schrodinger

#Plot a sampling of wavefunctions:
gnuplot -persist schrod.gp >schrod.ps
gnuplot -persist schrod1.gp >schrod1.ps

#Plot energy levels vs. wall fuzziness
gnuplot schrod2.gp >schrod2.ps
gs -q schrod2.ps # -q means quiet (less text to stdout).

#foreach thick (0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0)
#   foreach eigenvalue (`ls wavefunc_T"$thick"_n*.txt`)
#      
#   end
#end

