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