Hi,
ich sehe keinen Fehler, aber Du must weiterrechnen. Für \( k = 4 \) eribt sich z.B. \( 5 a_5 + a_2^2 = 0 \), also \( a_5 = -\frac{1}{20} \).
Der nächste Term ungleich Null wäre bei \( n = 7 \) und es ergibt sich \( a_8 = \frac{1}{160} \). Dann der nächste bei \( n = 10 \) usw.
Eine Systematik habe ich nicht erkannt. Laut Wolfram hat es aber vielleicht was mit den Besselfunktionen zu tun.
Die Konvergenz ist auch sehr langsam, wie man an dem Bild sieht. Die angenäherte Lösung in blau lautet
$$ p(x) = \frac{1}{2}x^2 -\frac{1}{20}x^5 + \frac{1}{160}x^8 - \frac{7}{8800}x^{11} $$ Der rote Graph ist die korrekte Lösung, numerisch ermittelt.