Total decay rate of a dipole near a multilayer interface
integrand_mtot(d = 10, q, wavelength, epsilon = list(incident = 1.5^2, 1^2), thickness = c(0, 0))
Integrand without transformation of variables
Other integrands dipole: integrand_nr1,
integrand_nr2, integrand_nr3,
integrand_rad