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