function f_vec = quad_fun ( n, xy_vec )

%*****************************************************************************80
%
%% quad_fun() is a sample integrand function for triangulation_quad().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    22 January 2007
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer N, the number of evaluation points.
%
%    real XY_VEC(2,N), the evaluation points.
%
%  Output:
%
%    real F_VEC(N), the value of the integrand
%    function at the evaluation points.
%
  f_vec(1:n) = exp ( xy_vec(1,1:n).^2 + xy_vec(2,1:n).^2 );

  return
end
