Trova arrotondate ai millesimi le soluzioni del sistema  K(x,y) = U(x,y)  dove:
K(x,y) = 2*(x^4 - 2*x^3 + 3)*(y^4 - 3*y^2 + 4) - 7
U(x,y) = x^2 + y^2 - 4
aiutandoti con del software.
  [a lato i grafici di K e di U]

Il grafico di K è simmetrico rispetto all'asse x in quanto K(x,y) = K(x,−y).  Potrei quindi limitarmi a trovare le intersezioni della parte superiore del grafico di K col cerchio grafico di U.

   Il modo più semplice è ricorrere a WolframAlpha. Introducendo
2*(x^4-2*x^3+3)*(y^4-3*y^2+4)-7=0 & x^2+y^2-4=0
otteniamo il grafico a sinistra e le soluzioni approssimate:
x = 1.30192, y = 1.51823
x = 1.75529, y = 0.958628
x = 1.30192, y = -1.51823
x = 1.75529, y = -0.958628

Ho trovato le soluzioni con molte cifre. Per trovare quanto richiesto dall'esercizio le arrotondo:
(1.755, 0.959),  (1.302, 1.518),
(1.755, -0.959),  (1.302, -1.518).

Come procedere con R.

# Il grafico iniziale è stato ottenuto con:
K = function(x,y) 2*(x^4-2*x^3+3)*(y^4-3*y^2+4)-7
U = function(x,y) x^2+y^2-4
Plane(-2,2, -2,2)
CURVE(K, "blue"); CURVE(U, "brown")
# Quello sopra a sinistra è stato ottenuto con:
Plane(0,2, 0,2)
CURVE(K, "blue"); CURVE(U, "brown")
# Con successivi ingrandimenti:
... arrivo a (grafico al centro):
Plane(1.755285,1.755290, 0.958625,0.958630)
CURVE(K, "brown"); CURVE(U, "blue")
POINT(1.755287, 0.958628, "red")
# Con altri ingrandimenti ...
... arrivo a (grafico a destra):
POINT(1.301916, 1.518228, "red")