DSolve[y''[x] + 2 (x^2 - 1)/x^3 y'[x] + (x^2 + 1)/x^6 y[x] == 0, y[x], x]