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