Finite difference scheme for one system of nonlinear partial integro-differential equations

System of Maxwell equations is considered. Reduction to the integro-differential form is given. Existence, uniqueness and large time behavior of solutions of the initial-boundary value problem for integro-differential model with two-component and one-dimensional case are studied. Finite difference scheme is investigated. Wider class of nonlinearity is studied than one has been investigated before. FreeFem++ realization code and results of numerical experiments are given.