Python-сообщество

Доброго времени суток!))
Разбирался с 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

Отредактировано mixxx92 (Апрель 16, 2015 11:33:50)

В том что A*x и numpy.dot(A,x) это совершенно разные вещи. Рекомендую смотреть что вы решаете.

я и так и так пробовал, ничего не меняется