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