dsolve(diff(y(x),x) = (y(x)-exp(x))^2+exp(x),y(x))