DSolve[-((-1 + a^2*x^(2*n))*y[x]) + Derivative[2][y][x] == 0,y[x],x]