Code to compute the magnetoelectric response

Implemented Leung-Prodan formula to model and classify a disordered 3D
topological insulator model.