code: def H_n(self, n, x):
          
          n_array = np.zeros(n+1)
          n_array[n] = 1
          return hermite.hermval(x, n_array, tensor=False)
docs: constructs the Hermite polynomial of order n at position x (dimensionless)
  
          :param n: The n'the basis function.
          :type name: int.
          :param x: 1-dim position (dimensionless)
          :type state: float or numpy array.
          :returns:  array-- H_n(x).
          :raises: AttributeError, KeyError