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