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