08-Oct-2025 14:15:26 r8ut_test(): MATLAB/Octave version 6.4.0 Test r8ut(). r8ge_to_r8ut_test(): r8ge_to_r8ut() converts an R8GE matrix to R8UT format. The random R8GE matrix: Col: 1 2 3 4 Row --- 1 0.0954606 0.0279161 0.458938 0.291044 2 0.654761 0.860739 0.332972 0.967934 3 0.458925 0.251329 0.239728 0.492033 4 0.855824 0.023408 0.766155 0.926891 5 0.199325 0.530789 0.0806359 0.156327 The R8UT matrix: Col: 1 2 3 4 Row --- 1 0.0954606 0.0279161 0.458938 0.291044 2 0.860739 0.332972 0.967934 3 0.239728 0.492033 4 0.926891 5 r8ut_det_test(): r8ut_det() computes the determinant of an R8UT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.64632 0.883287 0.15621 0.249378 0.884457 2 0.133648 0.0378193 0.404395 0.688721 3 0.861032 0.672516 0.267055 4 0.217977 0.543825 5 0.477712 Determinant is 0.00774471 r8ut_indicator_test(): r8ut_indicator() sets up an indicator matrix in R8UT format Matrix rows M = 5 Matrix columns N = 4 The indicator matrix: Col: 1 2 3 4 Row --- 1 11 12 13 14 2 22 23 24 3 33 34 4 44 5 r8ut_inverse_test(): r8ut_inverse() computes the inverse of an R8UT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.164606 0.622652 0.509675 0.10143 0.643015 2 0.404028 0.609919 0.800203 0.234717 3 0.955959 0.473298 0.477913 4 0.0392809 0.807345 5 0.438412 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 6.07512 -9.36244 2.73442 142.091 -268.542 2 2.47507 -1.57914 -31.3933 58.2078 3 1.04607 -12.6041 22.0705 4 25.4576 -46.8808 5 2.28096 The product A * B: Col: 1 2 3 4 5 Row --- 1 1 0 0 -8.88178e-16 6.21725e-15 2 1 0 0 4.66294e-15 3 1 0 -2.22045e-16 4 1 0 5 1 r8ut_mm_test(): r8ut_mm() computes C = A * B for R8UT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 1 1 1 1 1 2 1 1 1 1 3 1 1 1 4 1 1 5 1 The product C = A * A Col: 1 2 3 4 5 Row --- 1 1 2 3 4 5 2 1 2 3 4 3 1 2 3 4 1 2 5 1 r8ut_mtm_test(): r8ut_mtm() computes C = A' * B for R8UT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 1 1 1 1 1 2 1 1 1 1 3 1 1 1 4 1 1 5 1 The product C = A' * A Col: 1 2 3 4 5 Row --- 1 1 1 1 1 1 2 1 2 2 2 2 3 1 2 3 3 3 4 1 2 3 4 4 5 1 2 3 4 5 r8ut_mtv_test(): r8ut_mtv() computes a matrix product b=A'*x for an R8UT matrix. The matrix A: Col: 1 2 3 4 Row --- 1 11 12 13 14 2 22 23 24 3 33 34 4 44 5 The vector X: 1: 1 2: 2 3: 3 4: 4 5: 5 The vector b=A'*x: 1: 11 2: 56 3: 158 4: 340 r8ut_mv_test(): r8ut_mv() computes a product b=A*x for an R8UT matrix. The R8UT matrix A: Col: 1 2 3 4 Row --- 1 11 12 13 14 2 22 23 24 3 33 34 4 44 5 Vector x: 1: 1 2: 2 3: 3 4: 4 Vector b = A*x: 1: 130 2: 209 3: 235 4: 176 5: 0 r8ut_print_test(): r8ut_print() prints an R8UT matrix. Matrix rows M = 6 Matrix columns N = 4 The R8UT matrix: Col: 1 2 3 4 Row --- 1 11 12 13 14 2 22 23 24 3 33 34 4 44 5 6 r8ut_print_some_test(): r8ut_print_some() prints some of an R8UT matrix. Matrix rows M = 4 Matrix columns N = 6 Some of the matrix: Col: 4 5 6 Row --- 1 14 15 16 2 24 25 26 3 34 35 36 r8ut_random_test(): r8ut_random() randomizes an R8UT matrix. Matrix order M, N = 5, 4 Matrix A: Col: 1 2 3 4 Row --- 1 0.0541768 0.179889 0.307603 0.839813 2 0.185606 0.349781 0.0572231 3 0.0789334 0.524889 4 0.502525 5 r8ut_sl_test(): r8ut_sl() solves a linear system A*x=b with R8UT matrix; Matrix order N = 5 The upper triangular matrix: Col: 1 2 3 4 5 Row --- 1 1 2 3 4 5 2 2 3 4 5 3 3 4 5 4 4 5 5 5 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 r8ut_slt_test(): r8ut_slt() solves a linear system A'x=b with R8UT matrix; Matrix order N = 5 The upper triangular matrix: Col: 1 2 3 4 5 Row --- 1 1 2 3 4 5 2 2 3 4 5 3 3 4 5 4 4 5 5 5 Solution to transposed system: 1: 1 2: 2 3: 3 4: 4 5: 5 r8ut_to_r8ge_test(): r8ut_to_r8ge() converts an R8UT matrix to R8GE format. The random R8UT matrix: Col: 1 2 3 4 Row --- 1 0.32259 0.818884 0.759875 0.160324 2 0.769207 0.599746 0.411228 3 0.0738303 0.534609 4 0.405603 5 The R8GE matrix: Col: 1 2 3 4 Row --- 1 0.32259 0.818884 0.759875 0.160324 2 0 0.769207 0.599746 0.411228 3 0 0 0.0738303 0.534609 4 0 0 0 0.405603 5 0 0 0 0 r8ut_zeros_test(): r8ut_zeros() zeros out space for an R8UT matrix. Matrix order M, N = 5, 4 Matrix A: Col: 1 2 3 4 Row --- 1 0 0 0 0 2 0 0 0 3 0 0 4 0 5 r8ut_test(): Normal end of execution. 08-Oct-2025 14:15:26