function r8vec_vector_triple_product_test ( )

%*****************************************************************************80
%
%% R8VEC_VECTOR_TRIPLE_PRODUCT_TEST tests R8VEC_VECTOR_TRIPLE_PRODUCT.
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    14 September 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8VEC_VECTOR_TRIPLE_PRODUCT_TEST:\n' );
  fprintf ( 1, '  R8VEC_VECTOR_TRIPLE_PRODUCT computes the \n' );
  fprintf ( 1, '  vector triple product of three R8VEC''s.\n' );
  fprintf ( 1, '  V4 = V1 x ( V2 x V3 ).\n' );

  v1 = [ 2.0; -1.0;  3.0 ];
  v2 = [ 1.0;  4.0; -5.0 ];
  v3 = [ 3.0; -2.0;  1.0 ];

  r8vec3_print ( v1, v2, v3, '  V1, V2 and V3:' );

  v4 = r8vec_vector_triple_product ( v1, v2, v3 );

  vexact = [ 62.0; 10.0; -38.0 ];
  r8vec2_print ( v4, vexact, '  V4 and Vexact' );

  return
end
