Pure Julia implementation of the finite difference frequency domain (FDFD) method for electromagnetics