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