dsolve((y(x)*f(y(x)^2+x^2)-x)*diff(y(x),x)+y(x)+x*f(y(x)^2+x^2) = 0,y(x))