Есть динамические уравнения Эйлера

для тела на орбите,

. Также из модели SGP4 берутся значения

для разных

. Нужно получать решения этой системы в зависимости от

. ode45/ode23, как я понял, использовать не получится, так как

для разных

разное, а некой функции

нету, есть только численные значения для

. Какие есть способы решения, кроме написания тысяч уравнений для разных

?