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