Предложена пространственно-временная численная модель для расчета распространения лазерного излучения, модулированного по частоте, в двухуровневой нелинейно-оптической среде с насыщением поглощения и дисперсии. Используемый численный метод основан на разложении поля по поперечной координате по модам Гаусса - Лагерра. В качестве исследуемых параметров рассмотрены зависимости выходной интенсивности пучка от времени и частоты при условиях, когда период модуляции сравним с атомными временами релаксации.