Function hc_compartment_freqs 

fn hc_compartment_freqs(theta: f64, b0: f64) -> (f64, f64)

Compartment frequencies (Hz) for fibre-to-B0 angle theta (radians). Returns (Δf_myelin, Δf_axon, Δf_extra); extra-axonal is exactly 0.