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