An enormous assortment of physical situations can be represented by the coupled nonlinear Schrödinger equations (CNLSEs). These equations have been exhibited in a fibre communication system to capture vibration diffusion over perpendicular polarization axes in nonlinear optical fibres and wavelength division of multiplexed systems. CNLSEs also model radiation inside photorefractive matter or crystals by modelling the radiation as water wave collision. Singular waves in these equations are overwhelmingly called vector solitons in the writing as singular waves, in general, consist of two components. In the previous physical situations, clashing vector solitons is a significant case. In the last decade, the system has been examined seriously. Studies have found that, in addition to crossing through the clash, vector solitons can also rebound each other or hinder one another. This article will present (CNLSEs) as the follows  :
where and are the wave amplitudes in pair polarizations, and c is the cross-phase modulation coefficient. The initial conditions are:
and the boundary conditions are
For linearly birefringent fibres , and for elliptically birefringent fibres, c
is a positive number. When , this system enhances two (decoupled) nonlinear Schrödinger equations (NLS), and when , the system is a Manakov equations. In both situations, we get an integrable system. In particular, the solitons of one polarization should pass through the pulses of the opposite interaction’s polarization in the absence of forming shadows, creating an elastic collision. We have a nonintegrable system for other values of c. The clash may be nontrivial, and several complex phenomena like reflection, transmission, trapping, and creation of other solitary waves can happen. Solitons in optical fibres can be described as nonlinear beats that propagate nearly distortion-free, extending into the distance and undergoing (elastic) collision. The sixth-order compact finite difference method is perfect to derive the analytical solution of the CNLS system. We determined the exact solution of the CNLSEs as follows:
Several improved numerical methods have been presented for the nonlinear Schrödinger equations     . Among the numerical techniques for CNLSEs, the majority of the improved methods apply finite difference methods       . Multi-symplectic methods have also been applied for resolving the CNLS systems recently    . Good outcomes have boon obtained in terms of the solutions and the conserved quantities. Wang  presents a numerical solution of the single and coupled nonlinear Schrödinger equation using split-step finite difference method. Galerkin method is used to solve CNLSE in  and , and the outcomes are superior in comparison with those of previous numerical works. Ismail using Collocation method, produced very accurate results. Concerning the finite element method, many articles have been published on nonlinear Schrödinger equations  . Many articles as  and  used Runge-Kutta of order 4 to solve CNLSEs, and they obtained lower accuracy in space direction (second and fourth order) without using the Compact Finite Difference method which gives sixth order.
The conserved quantities of Equations (1) and (2) are follows :
1) Mass conservations:
2) Energy conservations:
2. Numerical Method
We suppose that, outside the interval , the solution of the system (1) and (2) is negligible. Therefore, to study this system numerically, we prefer to consider the system
where and are the wave gage in two polarizations and the cross-phase modulation coefficient is c, with
as initial conditions and
as boundary conditions. To conclude the numerical method to solve 8 and 9, the rectangular mesh point with coordinates
which cover the region . Numerically, we prefer to disjoint the complex functions and into their real and imaginary portions by
where are real functions, and we write,
where we have assumed
and and are all real functions.
By substituting (10) into (8) and (9), the next system is obtained:
The previous system in a matrix-vector is
At the grid point , the exact solution is denoted by while the numerical one is denoted with . The Compact Finite Difference (CFD) formula is derived by :
Now, by Taylor Expansion, we can have the truncation error as the following :
if we solve and , we get
so the truncation error becomes
If then and , which gives the explicit scheme in order four for the second derivative. Furthermore, when , the scheme becomes sixth-order accurate, in this case and .
By using (12)-(15), Equations (17) can be written as
The above we can rewrite as
where , and .
This system in a compact form will be written as
The solution of (26), using the explicit RK4, can be displayed as follows:
The proposed scheme is of sixth order in space and fourth order in time. The previous scheme, which applies the Runge-Kutta of order 4 method, is conditionally stable.
3. Numerical Results
This section presents the numerical results for the aimed scheme. This method’s accuracy is tested with by studying its conservation properties and the truncation error using and norms. Trapezoidal rule is used to calculate the conserved quantities. To reliably estimate the smooth solution of CNLSEs, discrete conservation laws are critical. The accuracy is scaled by using the error norm and error norm.
In the following sections, we will analyze the subsequent problems that arose in studying the properties of the obtained scheme.
3.1. Single Soliton
In this experiment, the initial condition is determined by
where , c and are constants.
To estimate the numerical solution, we use the following parameters:
In Table 1, we demonstrate the conserved quantities recovered from the proposed scheme. It is obvious that all conserved quantities are nearly conserved. Table 2 shows the accuracy of the method where and , . Figure 1 shows the development of a single soliton running to the right with velocity .
3.2. Interaction of Two Solitons
To examine the interaction of two solitons, we exercise the following initial conditions.
Table 1. Single solitons (conserved quantities).
Table 2. The errors of single soliton.
Figure 1. Single soliton.
These initial conditions draw two waves. The complete sum of the masses asymptotically can be shown as .
In this test we choose the parameters
This method is very nearly conserves all conserved quantities (see Table 3). In Figure 2, we demonstrate the interaction of the pair waves moving at the same velocities in opposite directions for (the same for ), which is called head-on interaction. It is recognized that both interacting and moving the interaction zone unaltered in shape and velocities. In Figure 3, we plot the contour of two solitons with the same parameters.
Figure 2. Interaction of two solitons.
Figure 3. Refection of two solitons.
Table 3. Interaction of two solitons (conserved quantities).
We have used a sixth-order finite difference method and Runge-Kutta of order 4 method to solve CNLSEs. The RK4 scheme is fourth order in time and sixth order in space, and it remains conditionally stable. Obviously, from the results, the scheme obtains fourth-order instantly without the need for an extra procedure to improve the accuracy in the temporal direction.
 Ismail, M.S. (2008) Numerical Solution of Coupled Nonlinear Schrödinger Equation by Galerkin Method. Mathematics and Computers in Simulation, 78, 532-547.
 Ismail, M.S. (2010) Collocation Method for Numerical Solution of Coupled Nonlinear Schrödinger Equation. In: AIP Conference Proceedings, American Institute of Physics, College Park, 1429-1432.
 Ismail, M.S. (1996) Finite Difference Method with Cubic Spline for Solving Nonlinear Schrödinger Equation. International Journal of Computer Mathematics, 62, 101-112.
 Muslu, G.M. and Erbay, H.A. (2005) Higher-Order Split-Step Fourier Schemes for the Generalized Nonlinear Schrödinger Equation. Mathematics and Computers in Simulation, 67, 581-595.
 Sheng, Q., Khaliq, A.Q.M. and Al-Said, E.A. (2001) Solving the Generalized Nonlinear Schrödinger Equation via Quartic Spline Approximation. Journal of Computational Physics, 166, 400-417.
 Ismail, M.S. and Taha, T.R. (2007) A Linearly Implicit Conservative Scheme for the Coupled Nonlinear Schrödinger Equation. Mathematics and Computers in Simulation, 74, 302-311.
 Ismail, M.S. and Alamri, S.Z. (2004) Highly Accurate Finite Difference Method for Coupled Nonlinear Schrödinger Equation. International Journal of Computer Mathematics, 81, 333-351.
 Ismail, M.S. and Taha, T.R. (2001) Numerical Simulation of Coupled Nonlinear Schrödinger Equation. Mathematics and Computers in Simulation, 56, 547-562.
 Sonnier, W.J. and Christov, C.I. (2005) Strong Coupling of Schrödinger Equations: Conservative Scheme Approach. Mathematics and Computers in Simulation, 69, 514-525.
 Aydin, A. and Karasozen, B. (2007) Symplectic and Multisymplectic Methods for the Coupled Nonlinear Schrödinger Equations with Periodic Solutions. Computer Physics Communications, 177, 566-583.
 Sun, J.Q. and Qin, M.Z. (2003) Multi-Symplectic Methods for the Coupled 1D Nonlinear Schrödinger System. Computer Physics Communications, 155, 221-235.
 Wang, H. (2005) Numerical Studies on the Split Step Finite Difference Method for the Nonlinear Schrödinger Equations. Applied Mathematics and Computation, 170, 17-35.
 Xu, Q.-B. and Chang, Q.-S. (2010) New Numerical Methods for the Coupled Nonlinear Schrödinger Equations. Acta Mathematicae Applicatae Sinica, English Series, 26, 205-218.