function fp = triangle_integrand_04 ( p_num, p, f_num )

%*****************************************************************************80
%
%% TRIANGLE_INTEGRAND_04 evaluates 4 integrand functions.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    15 August 2009
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Input, integer P_NUM, the number of points.
%
%    Input, real P(2,P_NUM), the evaluation points.
%
%    Input, integer F_NUM, the number of integrands.
%
%    Output, real FP(F_NUM,P_NUM), the integrand values.
%
  fp = zeros ( 4, p_num );

  fp(1,1:p_num) = p(1,1:p_num) .* p(1,1:p_num) .* p(1,1:p_num);
  fp(2,1:p_num) = p(1,1:p_num) .* p(1,1:p_num) .* p(2,1:p_num);
  fp(3,1:p_num) = p(1,1:p_num) .* p(2,1:p_num) .* p(2,1:p_num);
  fp(4,1:p_num) = p(2,1:p_num) .* p(2,1:p_num) .* p(2,1:p_num);

  return
end