docs: def gen_psi(self, x):
          """Generates the activity of the basis functions for a given 
          canonical system rollout. 
          
          x float, array: the canonical system state or path
          """
     
          if isinstance(x, np.ndarray):
              x = x[:,None]
          return np.exp(-self.h * (x - self.c)**2)