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 = sp_T_matrix(ns, ls)
V = V0 .* V_Gaussian.(R, l, ns, transpose(ns))
H = T + V

eigvals(H)