matrix_diagonal_test <- function ( )

#*****************************************************************************80
#
## matrix_diagonal_test() computes a matrix diagonal using einsum().
#
#  Licensing:
#
#    This code is distributed under the MIT license. 
#
#  Modified:
#
#    07 August 2026
#
#  Author:
#
#    John Burkardt
#  
cat ( '\n' )
cat ( '  matrix_diagonal_test():\n' )

A = matrix ( c (
    11, 12, 13,
    21, 22, 23,
    31, 32, 33 ), nrow = 3, ncol = 3, byrow = TRUE )
cat ( '\n' )
cat ( '  A:\n' )
print ( A )

D = einsum ( 'ii -> i', A )
cat ( '\n' )
cat ( '  D = einsum ( \'ii -> i\', A ):\n' )
print ( D )