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