powered by
integrand_mtot(d = 10, q, wavelength, epsilon = list(incident = 1.5^2, 1^2), thickness = c(0, 0))
integrand_nr1
integrand_nr2
integrand_nr3
integrand_rad