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