Я использую quad из scipy.integrate, чтобы получить интеграл в ограниченном диапазоне от объекта. предположим, что целевой объект находится под ударом:
Код: Выделить всё
∫expm(A*X).expm(B*X)dx
To solve this I have used blow code:
Код: Выделить всё
from scipy.integrate import quad
from scipy.linalg import expm
import numpy as np
def integrand(X, A, B):
return np.dot(expm(A*X),expm(B*X))
A = np.array([[1, 2], [3, 4]])
B = np.array([[1, 2], [3, 4]])
I= quad(integrand, 0, 1, args=(A,B))
Код: Выделить всё
TypeError: only length-1 arrays can be converted to Python scalars
Источник: https://stackoverflow.com/questions/547 ... hen-gettin