dsolve(diff(y(x),x) = -(-1/x*ln(y(x))+1/x/ln(x)*ln(y(x))-_F1(x))*y(x),y(x))