%% EM_constants
mu0 = 4e-7*pi;
c0 = 299792458;
e0 = 1/(mu0*c0^2);