Function calculate_limit_cycle_yaw_with_inertias 

pub fn calculate_limit_cycle_yaw_with_inertias(
    velocity_mps: f64,
    spin_rate_rad_s: f64,
    stability_factor: f64,
    spin_inertia: f64,
    transverse_inertia: f64,
) -> f64

Calculate the first-order flat-fire yaw-of-repose magnitude from the projectile inertias.

This is the classical gravity/gyroscopic reduction yaw = 4 * Iy * Sg * g / (Ix * |p| * V). Sg = 1 returns its finite stable-side limit; Sg < 1 returns 0.0 because no stable equilibrium exists. Crosswind is absent because it produces a damped transient, not persistent equilibrium yaw.