function value = smith1_f ( m, n, x )

%*****************************************************************************80
%
%% smith1_f() returns the function in the smith1 problem.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    21 January 2016
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Andrew Smith,
%    Fast construction of constant bound functions for sparse polynomials,
%    Journal of Global Optimization,
%    Volume 43, 2009, pages 445-458.
%
%  Input:
%
%    integer M, the number of variables.
%
%    integer N, the number of points.
%
%    real X(M,N), the points.
%
%  Output:
%
%    real VALUE(N), the value of the function at X.
%
  value(1:n) = ( ...
      3.0E+00 * x(1,1:n) .^2 .* x(2,1:n) .^3 .* x(3,1:n) .^4 ...
    +           x(1,1:n) .^3 .* x(2,1:n)     .* x(3,1:n) .^3 ...
    - 5.0E+00 * x(1,1:n)     .* x(2,1:n)     .* x(4,1:n) .^5 ...
    +           x(3,1:n)     .* x(4,1:n)     .* x(5,1:n) .^3 )';
 
  return
end


