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