# importing numpy and hermite modules
import numpy as np
from numpy.polynomial import hermite
 
# Creating a 3D array of coefficients 'C'
C = np.array([ [[0,1,2], [3,4,5]], [[6,7,8], [9,10,11]] ])
 
# Evaluating the 2 dimensional hermite
# series ant x,y using
# hermval2d() function
print(hermite.hermval2d([1,2],[1,2],C))

[[ 54. 180.]
 [ 63. 205.]
 [ 72. 230.]]

# importing numpy and hermite modules
import numpy as np
from numpy.polynomial import hermite
 
# Creating a 3D array of coefficients 'C'
C = np.arange(36).reshape(2,2,9)
 
# Evaluating the 2 dimensional hermite
# series ant x,y using
# hermval2d() function
x = [2,1]
y = [1,2]
print(hermite.hermval2d(x,y,C))

[[306. 288.]
 [321. 303.]
 [336. 318.]
 [351. 333.]
 [366. 348.]
 [381. 363.]
 [396. 378.]
 [411. 393.]
 [426. 408.]]