При численном моделировании динамики околозвездных дисков лагранжевыми методами возникает задача расчета траекторий твердых частиц. Возникает вопрос о выборе метода интегрирования и способа расчета коэффициента трения.

В данной работе исследовалась применимость разностных схем для решения жесткой системы обыкновенных дифференциальных уравнений, описывающей движение частицы по орбите под действием силы трения о газ и притяжения центральной звезды, так же исследовалась чувствительность решения к коэффициенту трения. В работе сравнивались полунеявная схема первого порядка и схема Рунге-Кутты четвёртого порядка, сравнивались коэффициенты трения Хендерсона и стандартный астрофизический.

В результате было установлено, что полунеявная схема первого порядка позволяет получать решения с приемлемой точностью, требуя меньший объём вычислительных затрат и памяти. Был определен диапазон параметров в котором решение чувствительно к коэффиценту трения.