The coupled nonlinear Schrödinger equation models several interesting physical phenomena. It represents a model equation for optical fiber with linear birefringence. In this paper we introduce a finite difference method for a numerical simulation of this equation. This method is second-order in...