using LinearAlgebra
include("ho_basis.jl")

l = 0
V0 = -10
R = 1
n_max = 20

ns = collect(0:n_max)
ls = fill(l, n_max + 1)

T = get_sp_T_matrix(ns, ls)

V_l(l, n1, n2) = V0 * V_Gaussian(R, l, n1, n2)
V = get_sp_V_matrix(V_l, ns, ls)

H = T + V

eigvals(collect(H))