# Script to simulate tunnelling in 1-D (animation)
# Isaac Lenton
# August 2017

set terminal gif animate delay 10
set output 'output.gif'
set xrange [-10:10]
set x2range [-10:10]
set yrange [0.0:25.0]
set y2range [0.0:1.0]
set xlabel "Position"

do for [i=0:1000-1] {
    plot "potential.dat" with lines title "Potential" axes x1y1, \
      sprintf("output_%06d.dat", i) with lines \
        title sprintf("|Wavefunction|^2 (at %d)", i) axes x2y2
}