Function base_pair 

fn base_pair(odd: bool, z: Complex, eta: Complex) -> [(Complex, f64); 2]

The base pair I_-3, I_-1 (odd, by the series) or I_-2, I_0 (even, in closed form), with error bounds in units of epsilon.