Computed the Leung-Prodan magnetoelectric response formula (J. Phys. A: Math. Theor. 46, 085205 (2013)). Modeled and classified a 3D topological insulator model.