dsolve(diff(y(x),x) = ((y(x)-ln(x)-Ci(x))^2+cos(x))/x,y(x))