Nimm doch den Ansatz und leite ihn zweimal ab
$$\begin{aligned}f(x) &= \sum_{k=0}^{\infty} c_k x^k \\ f'(x) &= \sum_{k=1}^{\infty} k c_k x^{k-1} \\ f''(x) &= \sum_{k=2}^{\infty} k(k-1) c_k x^{k-2} \end{aligned} $$
dann setzte die Funktionen in die DGL ein
$$ (1-x^2) f''(x) - 2xf'(x) + 2f(x) = 0 $$ $$ (1-x^2) \left( \sum_{k=2}^{\infty} k(k-1) c_k x^{k-2} \right) - 2x \left( \sum_{k=1}^{\infty} k c_k x^{k-1}\right) + 2\sum_{k=0}^{\infty} c_k x^k = 0 $$ das Ganze wird jetzt so sortiert, dass die Koeffizienten vor gleich großen Exponenten von \(x\) zusammen gefasst werden: $$ \left( \sum_{k=2}^{\infty} k(k-1) c_k x^{k-2} \right) - \left( \sum_{k=2}^{\infty} k(k-1) c_k x^{k} \right) - 2 \left( \sum_{k=1}^{\infty} k c_k x^{k} - \sum_{k=1}^{\infty} c_k x^k\right) + 2c_0= 0 $$ $$ \left( \sum_{k=0}^{\infty} (k+2)(k+1) c_{k+2} x^{k} \right) - \left( \sum_{k=2}^{\infty} k(k-1) c_k x^{k} \right) - 2 \left( \sum_{k=2}^{\infty} (k-1) c_k x^{k}\right) + 2c_0= 0 $$ $$ 2c_2 + 6c_3 x + \left( \sum_{k=2}^{\infty} (k+2)(k+1) c_{k+2} x^{k} \right) - \left( \sum_{k=2}^{\infty} (k+2)(k-1) c_k x^{k} \right)+ 2c_0= 0 $$ $$ 2(c_0 + c_2) + 6c_3 x + \sum_{k=2}^{\infty} (k+2)\left[ (k+1)c_{k+2} - (k-1) c_k\right] x^{k} = 0 $$ damit obige Gleichung erfüllt ist, muss jeder Koeffizient =0 sein. Weiter ist als Anfangsbedingung \(f(0)=f'(0)=1\) gegeben. Daraus folgt, dass \(c_0=1\) und \(c_1=1\) ist. Die weiteren \(c_k\) sind dann: $$ \begin{aligned} \Rightarrow c_2 &= -c_0 = -1 \\ c_3 &= 0 \\ c_{k+2} &= \frac{k-1}{k+1} c_k\end{aligned}$$ d.h. alle \(c_k\) mit geradem \(k\) sind ungleich 0 ohne Beschränkung. Also ist die Funktion beliebig oft differenzierbar. (nein deshalb nicht s.Kommentar unten)
Nachtrag: Das Polynom lautet demnach:
$$f(x)= 1 + x -\sum_{n=1}^{\infty} \frac{1}{2n-1}x^{2n}$$
Ich habe die Differenzialgleichung nochmal bei Wolfram Alpha eingegeben und als Lösung hat er mir eine Legendre's Gleichung ausgespuckt
$$f(x)= 1 + x + \frac12 x \big( \ln(1 - x) - \ln(1 + x) \big)$$
Zusätzlich habe ich einen nummerischen Algorithmus (nach Runge Kutta) auf die DGL angesetzt und alle drei Funktionen sehen als Graph so aus:
Die nummerische Lösung (f(x) blau; f'(x) orange) und das Polynom (grau; bis 8.Grades) stimmen sehr gut überein. In der Nähe von 1 werden sie aus naheliegenden Gründen 'unscharf'. Warum aber die Lösung von Wolfram Alpha (gelb) so daneben liegt, dafür habe ich keine Erklärung :-/
Gruß Werner