master
calculate_analytical_cross_sections
def calculate_analytical_cross_sections(eps_cyl, eps_bg, wavelength, radius, N=50):