dsolve(diff(y(x),x) = -1/216*x/(x^2+1)^4*(-513-432*x-378*y(x)-864*x^4-756*x^3-1134*x^2-594*x^2*y(x)+720*x^3*y(x)-540*y(x)^2-216*y(x)^3-1296*x^2*y(x)^2-648*x^2*y(x)^3-972*x^4*y(x)^2-216*x^6*y(x)^2+432*x^3*y(x)^2-456*x^6-576*x^5-288*y(x)*x^8+288*y(x)*x^7-288*y(x)*x^6+864*y(x)^2*x^5-648*y(x)^3*x^4-216*y(x)*x^4-144*x^7+64*x^9-96*x^8+432*y(x)^2*x^7+1008*x^5*y(x)-216*x^6*y(x)^3),y(x))