code: def doublewell_eigs(n_grid, lag_time=1):
      
      return _brownian_eigs(n_grid, lag_time, DOUBLEWELL_GRAD_POTENTIAL,
                            -np.pi, np.pi, reflect_bc=True)
docs: Analytic eigenvalues/eigenvectors for the doublwell system
  
      TODO: DOCUMENT ME