Математическое моделирование систем и явлений / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1986.

Вкачественачальногораспределениячастицпоколебательнымуровням бралось больцмановскоесколебательнойтемпературойТ , существеннобольшейначаль нойпоступательнойтемпературы газаТ0. ЗадачаКошидлясистемыуравнений(I), (2) решаласьчисленонаЭВМ ЕС-1045. ПриэтомварьировалисьзначенияТ0 , Тѵ , ѵ*. Длярешенияприменялся методГира/5/, разработанныйдляинтегрирования"жестких" системдифферен циальныхуравнений. Точностьвычисленийвыбираласьравной10”^, изменение шагаинтегрированияповремениосуществлялосьавтоматическиизависелоот ошибки, накопленнойнашаге. Результатырасчетовпредставленынарис.1-5. Видно, чтообщийхарактер эволюци функциираспределенияобусловленсуществованиемдвухстадийколеба тельнойрелаксации: стади формированияквазистационарногораспределенияи РисЛ. Т0=75СЖ, ТѴ-5500К. 1-8 - 0, 1.8. 6.7, 71, 431, 1770, 5390, ь490 мкс. Рис. 2. То =500К, Тѵ=5000К. 1-7 - 0.4, 3.8, 9.3, 31, 285, 6094, 15050 мкс. о• Рис.З. Т0=300К, ТѴ=5800К. І-ІО - 8.6, 12.7, 16.4 , 23 , 36 , 50 , 56 , 61, 72, 98 мкс. 69