> g = function(w) {
    x = w[1]
    y = w[2]
    x^2 + x * y + y^2 - 3 * x - 3 * y + 4
}

> nlm(g, c(2, 3))
$minimum
[1] 1

$estimate
[1] 0.9999997 0.9999997

$gradient
[1]  0.000000e+00 -1.332268e-09

$code
[1] 1

$iterations
[1] 3