function asa053_test01 ( )

%*****************************************************************************80
%
%% asa053_test01() generates a random Wishart variate.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    16 April 2014
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'asa053_test01():\n' );
  fprintf ( 1, '  Use wshrt() to generate a single Wishart deviate.\n' );

  d = [ ...
    3.0, ...
    2.0, 4.0, ...
    1.0, 2.0, 5.0 ];
  n = 1;
  np = 3;
 
  fprintf ( 1, '\n' );
  fprintf ( 1, '  The number of variables is %d\n', np );
  fprintf ( 1, '  The number of degrees of freedom is %d\n', n );

  r8utp_print ( np, d, '  The upper Cholesky factor:' );

  sa = wshrt ( d, n, np );

  r8pp_print ( np, sa, '  The sample matrix:' );

  return
end

