pub fn double_well_1d(x: f64) -> f64
1D double well potential: V(x) = (x^2 - 1)^2 Symmetric, minima at +/- 1, barrier at x=0