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