AbstractSince the problem of waveform distortion in optical communications can be overcome by optical soliton, very wide band transmission becomes possible. It is foreseen that soliton propagation is the candidate for future optical systems. The soliton behavior can be expressed by the coupled nonlinear Schrödinger equations.In this equation, taking the nonlinear oscillatory or loss terms due to the mutual parametric effects of each polarization component into account, the analysis will become very complicated. Concerning this point, in most analytical procedures, such oscillatory terms are neglected. However, the contribution of these terms are large in optical fibers of low birefringence and, therefore, they cannot be neglected.In order to establish a unified and effective approach, taking these perturbation and oscillatory terms into account, an analytical procedure based on the finite element method is proposed. Ultrashort pulse transmission in various fibers and the effects of the birefringence upon soliton pulses are investigated. The correctness of the present method is verified and several new results are obtained. Finally, in case of subpicosecond pulse transmission in optical fibers of low birefringence, it is shown that the contribution of the above mentioned terms is quite remarkable.