Function impulse_from_momentum_change 

pub fn impulse_from_momentum_change(
    initial_momentum: f64,
    final_momentum: f64,
) -> Option<f64>

Computes impulse from a change in momentum using J = p_final - p_initial.

Returns None when either input is not finite or when the computed impulse is not finite.