In this work we presented the results of numerical analysis of higher order nonlinear Schrodinger equation (NLSE) performed by the method of lines (MOL). In the higher order NLSE we take into account the imaginary Raman term. It is studied the impact of the dissipative Raman term on the femtosecond first order soliton stability. According to the results of numerical simulations, when real part of Raman term is equal to zero in the wavelength range, where second and third order dispersion negative, there exist analytic solitary-wave solutions. Our numerical simulations closely follow the published analytical data.© (2010) COPYRIGHT SPIE--The International Society for Optical Engineering. Downloading of the abstract is permitted for personal use only.