Sea $J(x)$ y $N(x)$ sean dos soluciones linealmente independientes de la ecuación
$$x^2y''+xy'+(-1+x^2)y=f(x)$$
Se sabe que la solución general de la ecuación viene dada por la fórmula
$$y=-J(x)\int N(x) xe^{-x^2} dx +N(x)\int J(x) xe^{-x^2} dx$$
Entonces
A) $ f(x)=x^2e^{-x^2} $
B) $ f(x)=xe^{-x^2} $
C) $ f(x)=e^{-x^2} $
D) $ f(x)=x^{-1}e^{-x^2} $
E) Ninguna de las anteriores
No estoy seguro de cómo encontrar el Determinante Wronskiano de $J(x)$ y $N(x)$ porque por lo que pienso la solución sería en forma de serie de potencias. Si $W(x)$ entonces puedo sustituirlo por $\frac{f(x)}{W(x)}= xe^{-x^2}$ pero no creo que sea así como se supone que hay que hacerlo.
¿Puede sugerir alguna forma de abordar este problema? Muchas gracias.