dsolve(diff(y(x),x) = (y(x)-a*ln(y(x))*x+x^2)/(-y(x)*ln(y(x))-y(x)*ln(x)-y(x)+a*x)*y(x)/x,y(x))