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