function s = MultiVarNewt(x_0, n)

for i = 1 : n
    v = DI(x_0) \ -I(x_0);
    x_1 = x_0 + v;
    x_0 = x_1;
end

s = x_0;
    
    
