У меня есть система ODE с некоторыми параметрами, назовем их a, b и c. На данный момент они являются только переменными, поэтому я использую этот код для вычисления моих уравнений:
def fct(y, t, a, b, c) :
ydot[1] = ..
ydot[2] = ..
return ydot
Чтобы решить его, я называю это так:
solution = odeint(fct, yinit, t, args = (a,b,c,))
, и он прекрасно работает, но теперь у меня есть та же самая переменная 'a', которая развивается со временем, но у меня нет функции a (t) для нее, просто вектор значения. Я просто хочу реализовать массив (очевидно, размер t), который содержит эволюционирующее значение 'a'
a = np.zeros((len(t),1)) # of course it's not full of zeros in my programm
Возможно ли использовать эти библиотеки?