Разбирался с scipy и numpy, понял как решать системы дифференциальных уравнений , если у нас есть фиксированное количество уравнений(например,2)
Когда попытался решить уравнение n*n(задавая матрицу коэффициентов А), у меня выскочила ошибка
Вот код:
from scipy import integrate import numpy as np A=np.array([[3,2],[1,-1]]) def func(t,x): f=np.array(A*x) return f y0=np.array([0,1]) t=np.linspace(0,1,10) result = integrate.odeint(func,y0,t) print result
Вот ошибка:
ValueError: object too deep for desired array
odepack.error: Result from function call is not a proper array of floats.
ValueError: object too deep for desired array
odepack.error: Result from function call is not a proper array of floats.
Может, кто-нибудь знает в чем дело?