Using this great post
of J.M. which defines FindAllCrossings2D
f[x_, y_] := (x^2 + y^2 - 4)
g[x_, y_] := (y - x^2 + 2 x - 1)
pts = FindAllCrossings2D[{f[x, y], g[x, y]}, {x, -7/2, 4}, {y, -9/5,
21/5}, Method -> {"Newton", "StepControl" -> "LineSearch"},
PlotPoints -> 85, WorkingPrecision -> 20] // Chop;
pl=ContourPlot[{f[x, y], g[x, y]}, {x, -7/2, 4}, {y, -9/5, 21/5},
Contours -> {0}, ContourShading -> False,
Epilog -> {AbsolutePointSize[6], Red, Point /@ pts}]

we can check it works
NSolve[{(x^2 + y^2 - 4) == 0, (y - x^2 + 2 x - 1) == 0}, {x,y}, Reals]
(* {{x->-0.399864,y->1.95962},{x->1.85894,y->0.737785}} *)
Show[pl,
Graphics[{AbsolutePointSize[12], Purple,
Point[{x, y}] /.
NSolve[{(x^2 + y^2 - 4) == 0, (y - x^2 + 2 x - 1) == 0}, {x, y}, Reals]}]]

NSolve[ {(x^2 + y^2 - 4) == 0 , (y - x^2 + 2 x - 1) == 0}, {x,y}, Reals]
enough for you? $\endgroup$