Quantum Tunneling

Numerical Domain
xmin = -4
xmax =  4
Nx = 4000
CFL = 0.125
Nt = 150000

Initial Conditions
p0 = 50.0
a = 0.25
psi = (2.0d0/pi/a**2)**0.25*exp(-((x+2.0d0)/a)**2)*exp(cmplx(0.0d0,p0*x))

Potential
if (|x| < 0.05):
    V = 1500
else:
    V = 0