Numerical Solution of Multiscale Electromagnetic Systems