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