dsolve(diff(y(x),x)^2+2*f(x)*y(x)*diff(y(x),x)+g(x)*y(x)^2-(g(x)-f(x)^2)*exp(-2*int(f(xp),xp = a .. x)) = 0,y(x))