function fsolve_test03 ( )

%*****************************************************************************80
%
%% fsolve_test03 tests fsolve() on 4 nonlinear equations.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    09 April 2021
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'fsolve_test03\n' );
  fprintf ( 1, '  Use fsolve() to solve 4 nonlinear equations.\n' );

  options = optimoptions ( 'fsolve', 'Display', 'off' );

  x0 = zeros ( 4, 1 );
  fx0 = f3 ( x0 );
  r8vec2_print ( x0, fx0, '  Initial X and F(X)' );

  x = fsolve ( @f3, x0, options );
  fx = f3 ( x );
  r8vec2_print ( x, fx, '  Final X and F(X)' );

  return
end

