Plusieurs choses :
-le calcul de Luc est juste et, cherchant une solution sous forme de série entière, on tombe sur cette classe de fonctions ;
-même si la démonstration de Cauchy-Lipschitz te semble ardue, c'est un théorème extrêmement important qu'il te faut savoir maîtriser. Ce que dit le théorème de Cauchy-Lipschitz linéaire, c'est que si tu as une EDO linéaire homogène d'ordre n, réglée en
})
(c'est-à-dire que le coeff devant est 1 ; ici, ce sera donc le cas sur chacun de tes sous-intervalles quitte à diviser par le coefficient devant

), l'espace de tes solutions forme un espace vectoriel de dimension

. Et que par ailleurs, si tu fixes
, f'(y_0), ..., f^{(n-1)}(y_0))
, tu as une unique solution qui prend les valeurs que tu as imposées.
Autrement dit, c'est un théorème bien pratique : tu peux chercher un peu comme tu veux les solutions (en trouver une au détour d'un buisson), si tu en trouves suffisamment, Cauchy-Lipschitz te dit que tu n'as plus besoin de chercher autre chose.
En pratique, ici : l'équation est de degré 2, l'espace des solutions doit donc être de dimension 2. La fonction
=\frac{x}{(1-x)^2})
que Luc te donne (et tous ses multiples scalaires) te donne une solution. Il manque une solution linéairement indépendante pour couvrir tout l'espace (ensuite, l'ensemble des solutions sera l'ensemble des combinaisons linéaires de ces deux fonctions). Une façon classique de chercher une solution linéairement indépendante une fois qu'on en a une est de la chercher comme multiple de l'autre, i.e. de chercher
=\lambda(x)\times y_1(x))
. Cherche-la sous cette forme, écris ce que ça veut dire pour
)
et conclus.