Python ODEINT проблемы с арг

Question

Я относительно новым для Python и пытается использовать его, чтобы решить проблему Integrator

x' = - L * x 

Где L ì лапласиан матрицы, то есть матрица представление графа. Это часть моего кода:

def integrate_cons(x, t, l): 
    xdot = -l*x 
    return xdot; 

t = np.linspace(0, 10, 101) 

#laplacian is a 3x3 matrix 
#initial_condition is a vector 
solution = odeint(integrate_cons, initial_conditions, t, args=(laplacian,)) 
print solution 

У меня возникли проблемы с передачей матрицы как аргумента в odeint. Как я могу решить?

Answer 1

import numpy as np 
from scipy.integrate import odeint 

def integrate_cons(x, t, l): 
    # unless you use np.matrix which I never do, you have to use np.dot 
    xdot = -np.dot(l,x) 
    return xdot; 

t = np.linspace(0, 10, 101) 
# just a random matrix 
l = np.random.rand(3,3) 
# initial conditions 
x0 = np.array([1,1,1]) 
#laplacian is a 3x3 matrix 
#initial_condition is a vector 
solution = odeint(integrate_cons, x0, t, args=(l,)) 
print(solution) 

Посмотрите на scipy cookbook для примеров.