Я не уверен, как использовать scipy.integrate.RK45 или scipy.integrate.BDF. Я пытаюсь решить следующую проблему:
x'(t) = -0.04*x(t) + 10^4*y(t)z(t)
y'(t) = 0.04*x(t) - 10^4*y(t)*z(t) - 3*10^7*(y(t))^2
z'(t) = 3*10^7*(y(t))^2
x(0)= 1
y(0)=z(0)=0
Пока что я не нашел никаких полезных ссылок, чтобы узнать, как использовать эти решатели од. Любая помощь или простые ссылки будут оценены.
Обновлено: z'(t) и x(0)=1