dsolve(diff(y(x),x) = y(x)/x*(-3*x^3*y(x)-3+y(x)^2*x^7)/(x^3*y(x)+1),y(x))