# importing numpy module and polyval3d function
import numpy as np
from numpy.polynomial.polynomial import polyval3d
  
# creating an 4d array of coefficiets 'C' 
# using numpy
C = np.arange(24).reshape(2,2,3,2)
  
# Now using polyval3d function we are 
# evaluating the 3D polynomial at points
# (x,y,z)
print(polyval3d([2,1],[1,2],[2,3], C))

[[ 582. 1032.]
 [ 624. 1110.]]

# importing numpy module and polyval3d function
import numpy as np
from numpy.polynomial.polynomial import polyval3d
  
# creating an 4d array of coefficiets 'C'
C = np.arange(72).reshape(3,2,6,2)
  
# Now using polyval3d function evaluate 
# the 3D polynomial at points (x,y,z)
x = [4,1]
y = [1,2]
z = [2,3]
  
print("The result of evaluation by polyval3d function is : \n",
      polyval3d(x, y, z, C))

The result of evaluation by polyval3d function is : 
 [[146412. 134370.]
 [149058. 137646.]]