Avec ce genre de somme, perso j'ai envie de poser
f(x) = somme(k=1, n, C(n,k)x^k/k²).
On remarque au passage que ta somme vaut f(1).
On remarque que f''(x) = somme(k=2, n, C(n,k)x^(k-2) (k-1)/k)
= somme(k=2, n, C(n,k)x^(k-2)) - somme(k=2, n, C(n,k)x^(k-2)/k)
= [(1+x)^n - 1 - x]/x² - (f'(x) - n)/x.
(J'espère pas m'être planté sur le dernier terme, j'en suis pas très sûr.
)
Bon bah maintenant que t'as une équa diff du premier ordre sur g := f', t'as plus qu'à la résoudre. Posons
b(x) = [(1+x)^n - 1 - x]/x² + n/x, on a
g'(x) = -g(x)/x + b(x).
Solution de l'équation homogène : C/x.
Perso, ça m'incite à faire encore un autre changement de fonction. Posons
h(x) = x*g(x).
Notre nouvelle fonction h est solution de
h'(x) = x*b(x) = [(1+x)^n - 1 - x]/x + n.
On a donc
h(x) = int_0^x ([(1+t)^n - 1 - t]/t + n) dt,
puis
g(x) = int_0^x ([(1+t)^n - 1 - t]/t + n) dt/x,
puis
f(x) = int_0^x int_0^t ([(1+s)^n - 1 - s]/s + n) ds/t dt.,
car f(0) = 0 si on se reporte à la définition de f.
Bon, c'est pas beaucoup plus explicite que ce que tu as obtenu, m'enfin ça a le mérite de ne faire intervenir que des plynômes, si on examine soigneusement le résultat. 