Le calcul est à moitié fait (je n'ai pas vérifié les formules); il suffit d'introduire, dans l'expression de l'intégrale simple l'expression du carré de l'amplitude de la fonction d'onde radiale, qui donnée par le produit d'un polynôme de degré (n - 1) par une exponentielle décroissante:
Rn, l(r) = Pn-1(r).Exp(-Zr/na0) (a0 désignant le rayon de la première orbite de Bohr).
L'intégrale est toujours convergente, mais l'établissement de son expression est difficile.
As-tu déjà vérifié le résultat particulier correspondant à k = 1 ?