dsolve(diff(y(x),x) = 4*x*(a-1)*(a+1)*(-y(x)^2+a^2*x^2-x^2-2)/(-4*y(x)^3+4*a^2*x^2*y(x)-4*x^2*y(x)-8*y(x)-a^2*y(x)^6+3*a^4*y(x)^4*x^2-6*y(x)^4*a^2*x^2-3*a^6*y(x)^2*x^4+9*y(x)^2*a^4*x^4-9*y(x)^2*a^2*x^4+a^8*x^6-4*a^6*x^6+6*a^4*x^6-4*a^2*x^6+y(x)^6+3*x^2*y(x)^4+3*x^4*y(x)^2+x^6),y(x))