Magnetoelectric Responses of a Topological Insulator

Implemented Leung-Prodan formula to model and classify the disordered
3D Fu-Kane-Mele model in EMV6.f90. Computed the energy gap of the
model in EMVGap.f90.

MEresponse.pdf shows the numerical results for various disorder
simutlations.