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