@ -90,7 +90,7 @@ function calculate_Vs(s::system{T}, V_twobody::Function, ϕ::T, n_image::Int)::A
end
Vs[i] *= s.multiplicity[i[1]]
Vs[i] *= s.multiplicity[s.unique_i[i[1]]]
return Vs