I the file freefemNewmarkEMwave was updated. The program now is solved with initial conditions that are a gausian spatial function and no loads. That is to verify that it is working properly (and so it seems...)
Todo1: Solve a real EM Scattering problem.
Todo2: Absorbing boundary conditions might have to be implemented.