fprintf ( 1, '\n' );
fprintf ( 1, 'multi2_newton_system2:\n' );
fprintf ( 1, '  Use Newton''s method for a root of the multi2 function.\n' );

global print_f_norm
global print_r
global print_x

print_f_norm = true;
print_r = false;
print_x = true;

x = [ 2.0; 2.0 ];
xtol = 0.0000001;
ftol = 0.0000001;
itmax = 50;

[ x, it ] = newton_system2 ( @(x)multi2(x), @(x)multi2_fp(x), x, xtol, ftol, itmax );

fprintf ( 1, '\n' );
fprintf ( 1, '  Number of steps = %d\n', it );
fprintf ( 1, '  x = (%g, %g)\n', x );
fprintf ( 1, '  f(x) = (%g, %g)\n', multi2(x) );

