Аннотация:Для моделирования колебаний холодной плазмы как в нерелятивистском случае, так и с учетом релятивизма, предложены и обоснованы численные алгоритмы высокой точности. Спецификой подхода является использование лагранжевых переменных для приближенного решения задачи, сформулированной в эйлеровых переменных. Основной результат представлен теоремами о сходимости предложенных алгоритмов относительно малых параметровдискретизации независимых эйлеровых переменных. Численные эксперименты наглядноиллюстрируют полученные теоретические результаты. В частности, проведено моделирование известного эффекта опрокидывания плазменных колебаний и подтверждено, что он имеет характер градиентной катастрофы