code: def H_n(self, n, x):
          
          if not self._interpolation:
              n_array = np.zeros(n+1)
              n_array[n] = 1
              return self.hermval(x, n_array, tensor=False)  # attention, this routine calculates every single hermite polynomial and multiplies it with zero (exept the right one)
          else:
              return np.interp(x, self.x_grid, self.H_interp[n])
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