Function stern_gerlach 

pub fn stern_gerlach(
    state: &SpinState,
    measurement_axis: Vec3,
    rng_val: f64,
) -> (SpinState, i8)

Stern-Gerlach measurement along an axis. Returns (collapsed state, outcome +1 or -1).