19-Aug-2022 19:53:54 r8but_test(): MATLAB version Test r8but(). R8BUT_DET_TEST R8BUT_DET gets the determinant of an R8BUT matrix. Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.260303 0.00183837 0.109957 0.218418 2 0.912484 0.897504 0.043829 0.956318 3 0.113664 0.350752 0.633966 4 0.351629 0.0945448 5 0.822887 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.267132 0.859097 0.401306 0.0661187 7 0.692066 0.840847 0.754673 0.257578 8 0.561662 0.123104 0.797287 9 0.861216 0.00751236 10 0.453794 Determinant = 0.000317009 R8BUT_INDICATOR_TEST R8BUT_INDICATOR sets up an R8BUT indicator matrix. Matrix order N = 6 Upper bandwidth MU = 2 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 2 22 23 24 3 33 34 35 4 44 45 5 55 6 Col: 6 Row --- 6 66 R8BUT_MTV_TEST R8BUT_MTV computes b=A'*x, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 2 0.0617272 0.257578 0.561695 0.956318 3 0.449539 0.109957 0.415307 4 0.401306 0.043829 5 0.754673 x: 1: 1 2: 2 3: 3 4: 4 5: 5 b=A'*x: 1: 0.633966 2: 0.189573 3: 2.69328 4: 3.2769 5: 7.10724 R8BUT_MV_TEST R8BUT_MV computes b=A*x, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 2 0.0617272 0.257578 0.561695 0.956318 3 0.449539 0.109957 0.415307 4 0.401306 0.043829 5 0.754673 x: 1: 1 2: 2 3: 3 4: 4 5: 5 b=A*x: 1: 4.1284 2: 7.92456 3: 3.86498 4: 1.82437 5: 3.77337 R8BUT_PRINT_TEST R8BUT_PRINT prints an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 2 22 23 24 25 3 33 34 35 4 44 45 5 55 R8BUT_PRINT_SOME_TEST R8BUT_PRINT_SOME prints some of an R8BUT matrix. Matrix order N = 10 Upper bandwidth MU = 3 Rows 1:4, Cols 2:4: Col: 2 3 4 Row --- 2 202 203 204 3 303 304 4 404 R8BUT_RANDOM_TEST R8BUT_RANDOM randomizes an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 2 0.0617272 0.257578 0.561695 0.956318 3 0.449539 0.109957 0.415307 4 0.401306 0.043829 5 0.754673 R8BUT_SL_TEST R8BUT_SL solves A*x=b, where A is an R8BUT matrix; Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.260303 0.00183837 0.109957 0.218418 2 0.912484 0.897504 0.043829 0.956318 3 0.113664 0.350752 0.633966 4 0.351629 0.0945448 5 0.822887 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.267132 0.859097 0.401306 0.0661187 7 0.692066 0.840847 0.754673 0.257578 8 0.561662 0.123104 0.797287 9 0.861216 0.00751236 10 0.453794 x: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 The right hand side: 1: 1.46752 2: 9.47438 3: 9.89089 4: 6.18147 5: 10.6654 6: 11.422 7: 20.9391 8: 13.5741 9: 7.82607 10: 4.53794 Solution x: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 R8BUT_SLT_TEST R8BUT_SLT solves A'*x=b, where A is an R8BUT matrix; Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.260303 0.00183837 0.109957 0.218418 2 0.912484 0.897504 0.043829 0.956318 3 0.113664 0.350752 0.633966 4 0.351629 0.0945448 5 0.822887 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.267132 0.859097 0.401306 0.0661187 7 0.692066 0.840847 0.754673 0.257578 8 0.561662 0.123104 0.797287 9 0.861216 0.00751236 10 0.453794 x: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 b: 1: 0.260303 2: 1.82681 3: 2.24596 4: 2.76485 5: 8.30715 6: 4.40631 7: 14.4935 8: 14.8636 9: 14.4152 10: 12.7869 x: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 R8BUT_TO_R8GE_TEST R8BUT_TO_R8GE converts a matrix from R8BUT to R8GE format. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 2 0.0617272 0.257578 0.561695 0.956318 3 0.449539 0.109957 0.415307 4 0.401306 0.043829 5 0.754673 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 0 2 0 0.0617272 0.257578 0.561695 0.956318 3 0 0 0.449539 0.109957 0.415307 4 0 0 0 0.401306 0.043829 5 0 0 0 0 0.754673 R8BUT_ZEROS_TEST R8BUT_ZEROS zeros an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0 0 0 0 2 0 0 0 0 3 0 0 0 4 0 0 5 0 r8but_test(): Normal end of execution. 19-Aug-2022 19:53:54