DSolve[-y[x] - x*Derivative[1][y][x] + f[x^2 + a*y[x]^2]*(x + a*y[x]*Derivative[1][y][x]) == 0,y[x],x]