The simulation of time-harmonic electromagnetic waves in wide domains with classic methods induces large linear systems which are poorly adapted to iterative resolution. We introduce an iterative quasi-Trefftz method for the Maxwell's equations, where the basis functions are parameterised by a polynomial approximation of an impedance trace space, and computed thanks to a Flux Reconstruction method, allowing to overtake some limitations of the classic choice of plane waves.