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