|
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)) |