Numerical model of long fiber Raman laser is proposed. The model based on the equations, describing the propagation of pump and Stokes waves, linear coupling of oppositely running waves due to scattering and its nonlinear interaction. The derivation of equations for slowly varying pulse envelopes uses the field decomposition in terms of spatial harmonics rather then commonly used temporal harmonics, which allows to avoid the twopoint boundary conditions, and to employ the numerical scheme of Courant–Isaakson–Rees.
