DSolve[Derivative[2][y][x] == -((((1 - a)*a*x^2 - b*(b + x))*y[x])/x^4),y[x],x]