Numerical Simulation of the Polarization Inside Lossless Media Characterized by Complex Susceptibilities Using Characteristic-based Method