Bon je reprends car j'avais pas écrit les choses comme toi et donc il y a une remarque que j'ai fait et qui est fausse.
Alors avec tes notations on a f(z)=a+bz+z^2 h(z) holomorphe.
La condition implique que -b + 2 h(z) - z + z h'(z)=3
La résolution de l'équation différentielle donne bien
=Ce^z/z^2-(b+z)(z+1)/z^2)
, avec

solution valable pour tout

Remarque 1. La procédure de résolution sur

est la même que

(rien ne change à la théorie).
Remarque 2. Ce qui change par rapport à ma remarque précédente que j'ai supprimée (perso je n'avais pas le même problème) c'est que parmi les solution obtenues certaine ne sont pas définies en z=0.
Alors avec un développement en z=0 , tu vas trouver une seule valeur de C compatible.
En remplaçant C tu va trouver une famille de solutions qui dépend des paramètres a et b.
Normalement tu vas retrouver le même que moi (peut être à la différence près que mes paramètres alpha et beta sont peut être différent, j'ai un peu oublié ma démarche d'hier).