Find $$f(x) = x + \frac{2}{3}x^3 + \frac{2\cdot4}{3\cdot5}x^5 + \frac{2\cdot4\cdot6}{3\cdot5\cdot7}x^7+\cdots +\infty\,,\quad|x|<1$$
My solution:
$$f'(x) = 1 + 2x^2 + \frac{2}{3}\cdot4x^4 + \frac{2\cdot4}{3\cdot5}\cdot6x^6 + \cdots + \infty $$ $$ = 1+x\left(\frac{d}{dx} xf(x)\right) $$ (By observation)
$$ = 1+x^2f'(x)+xf(x)$$
$\implies (1-x^2)\frac{dy}{dx} = 1 + xy$ (where $y=f(x), \frac{dy}{dx}=f'(x)$)
This reduces to a linear first order differential equation, which along with $f(0)=0$, gives $$f(x) = \frac{\sin^{-1}(x)}{\sqrt{1-x^2}}$$
Is there any other way to solve this question? Preferably along the lines of Taylor series, infinite GP, fubini theorem, etc.?