Here, for convenience, the vector in 2D is represented by a complex number: the real part represents the x-component and the imaginary part represents the y-component. How to get the Fourier transform is the key of this program. Here I used the a

this program calculates and plots the photonic bands for a 2D photonic crystal consisting of cylinders with circular cross-section and infinite height, arranged in a triangular lattice