pub fn wkb_approximation(potential: &[f64], energy: f64, x_grid: &[f64]) -> f64
WKB approximation for tunneling through an arbitrary potential. T ~ exp(-2/hbar * integral sqrt(2m(V(x)-E)) dx) from x1 to x2.