Forward modeling is an essential part of inverse modeling in geophysical problems that is used to calculate the forward response of the physical properties of the model. Through the physics of the problem, it is possible to simulate the ground potential response for a given electrode array. So in the first step, the partial differential equation for a two dimensional field with an arbitrary distribution of conductivity is solved within the framework of the finite difference numerical method. In the next step, based on the existing relationships between conductivity and induced polarization, the apparent induced polarization response is calculated and finally the accuracy of the algorithm in the form of the MATLAB programing language is evaluated by calculating theoretical models.