Я пытаюсь решить систему дифференциальных уравнений, используя solve_idp. Мои уравнения содержат параметры с временными единицами секунд, и я хочу интегрировать их в течение 20 лет с одинаковым временным шагом, чтобы изменение времени имело смысл.
Я попробовал следующее:
year_to_sec = 60*60*24*365.25
t_span = [0, 200*year_to_sec]
times = (np.linspace(0, 200, 20000))*year_to_sec
TS_0 = np.array([0, 0, 0, 0, 0, 0,
35, 35, 35, 35, 35, 35])
sol = solve_ivp(ode, t_span, TS_0,
method="RK23", t_eval=times, rtol=1e3)
Идея заключалась в том, чтобы сделать световой допуск, чтобы временной интервал оставался прежним, и определить times для указанных c точек, которые я хочу ( 100 баллов в год). Тем не менее, результаты для времени sol.t идут от 0 до 1 года.
Чего мне не хватает?