@ -1,4 +1,5 @@
using SpecialFunctions
include("helper.jl")
# Gaussian potentials in HO space
inv_factorial(n) = Iterators.prod(inv.(1:n))