Function gaussian_wavepacket 

pub fn gaussian_wavepacket(
    x0: f64,
    k0: f64,
    sigma: f64,
    x_grid: &[f64],
) -> Vec<Complex>

Generate a Gaussian wave packet centered at x0 with momentum k0.