Function velocity_after_impulse 

pub fn velocity_after_impulse(
    mass: f64,
    initial_velocity: f64,
    impulse: f64,
) -> Option<f64>

Computes final velocity after an impulse using v_final = v_initial + J / m.