We present the implementation of the Lacuna method, that removes a key difficulty that
currently hampers many existing methods for computing unsteady electromagnetic waves
on unbounded regions. Numerical accuracy and/or stability may deteriorate over long times
due to the treatment of artificial outer boundaries. We describe a developed universal
algorithm and software that correct this problem by employing the Huygens' principle and
lacunae of Maxwell's equations. The algorithm provides a temporally uniform guaranteed …