In this paper, the propagation of an electromagnetic pulse due to the interaction of laser pulse and plasma was investigated using a set of fluid relativistic equations and Maxwell’s equations in an unmagnetized collisionless common plasma. Using the multiple scale perturbation approach, the dispersion relation and the group velocity were derived, and finally it was shown that the evaluation of the components of vector potential is governed by two coupled-nonlinear Schrödinger (NLS) equations. It should be noted that using the relations of the phase and group velocities we understated that they are independent from plasma parameters. Then, using analytical methods, the solutions of the amplitude of the coupled Schrodinger equations for both components of the vector potential were obtained. Analytically, the modulation instability conditions in terms of plasma parameters have been investigated. Finally, by plotting the figures, the modulation stability and instability regions are numerically studied.