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