# import hermite_e libraries
from numpy.polynomial import hermite_e as h
 
# create a multidimensional array
# 'arr' of coefficients
arr = [[0, 1, 2],[3, 4, 5]]
 
# integrate a Hermite_e series using
# hermite_e.hermeint() function
print(h.hermeint(arr, m=2, k=[1, 2], lbnd=-1, axis=0))

[[2.         2.66666667 3.33333333]
 [1.         2.         3.        ]
 [0.         0.5        1.        ]
 [0.5        0.66666667 0.83333333]]

# import numpy and hermite_e libraries
import numpy as np
from numpy.polynomial import hermite_e
 
# create a multidimensional array
# 'arr' of coefficients
arr = np.arange(6).reshape(2,3)
 
# integrate a Hermite_e series using
# hermite_e.hermeint() function
print(hermite_e.hermeint(arr, axis = 0))

[[1.5 2.  2.5]
 [0.  1.  2. ]
 [1.5 2.  2.5]]