Рассмотрено применение метода конечных разностей во временной области (FDTD) для моделирования сверхбыстрой динамики излучения в двухуровневой активной среде. В основе метода лежат уравнения Максвелла и уравнение Блоха для диагонального элемента матрицы плотности. Проведено моделирование процесса образования n-импулъсов в инвертированной усиливаю¬щей среде, которое иллюстрирует перспективность данного подхода для применения в физике лазеров ультракоротких импульсов.