Function root_locus 

pub fn root_locus(
    tf: &TransferFunction,
    gains: &[f64],
) -> Vec<(f64, Vec<Complex>)>

Compute root locus: poles of 1 + K*H(s) for varying gain K.

Returns a vector of (gain, poles) pairs.