function matrix_column_sum_test ( )

%*****************************************************************************80
%
%% matrix_column_sum_test() computes matrix column sums using einsum().
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    08 August 2026
%
%  Author:
%
%    John Burkardt
%  
  fprintf ( 1, '\n' );
  fprintf ( 1, '  matrix_column_sum_test():\n' );

  A = [ ...
    11, 12, 13; ...
    21, 22, 23; ...
    31, 32, 33 ];

  fprintf ( 1, '\n' );
  fprintf ( 1, '  A:\n' );
  disp ( A );

  Ac = einsum ( 'ij -> j', A );
  fprintf ( 1, '  Ac = einsum ( "ij -> j", A ):\n' );
  disp ( Ac );

  return
end

