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