fn even_pair(z: Complex, eta: Complex) -> [(Complex, f64); 2]
I_-2 and I_0 in closed form by two complementary error functions, with the magnitude of their two terms as error bounds in units of epsilon.
I_-2
I_0