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