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