Automatic differentiation of a recursive partition function

Computing the energy of a bosonic gas in an harmonic trap using automatic differentiation.