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))