(Bon D'ieu, c'est pas vraiment clair ce que j'écris.
)
Sur [0,1], on doit pouvoir prendre P_n = (1-x)^n divisé par son intégrale sur [0,1], à vue de nez ça marche.
Ensuite, tu remarques que x |-> int(0,1,P_n(x-y)f(y)dy). Si tu regardes attentivement la tête cette intégrale à paramètre, tu verras que c'est... un polynôme en x.
Ensuite, les trois propriétés P_n > 0, int(0,1,P_n) = 1, et int(delta,1,P_n) -> 0 quand n -> +oo t'assurent que cette intégrale converge uniformément vers f lorsque n tend vers + inf.
C'est mieux, je pense. 