This paper investigates finite difference schemes for solving a system of the nonlinear Schrödinger (NLS) equations. Several types of schemes, including explicit, implicit, Hopscotch-type and Crank-Nicholson-type are defined. Cubic spline interpolation is used for solving time-shifting part of equations. The numerical results of the different solution methods are compared using two analytical invariant properties.