Function pair_coupling 

fn pair_coupling(
    gap_nm: f64,
    coupling_length_um: f64,
    wavelength_nm: f64,
    core_index: f64,
    cladding_index: f64,
) -> (f64, f64, f64)