On intègre par partie en posant

dU = -2ax e^{-ax^2}

V = \frac{x^{n+1}}{n+1}