with(LinearAlgebra);
CAb := <1, 1, 2;1, 2, 3;1, 3, 1>;
# Get the QR decomposition of CAb:
Q, R := QRDecomposition(CAb); R;
Error in R
Maple Matrix(3, 3, {(1, 1) = 1.7321, (1, 2) = 3.4642, (1, 3) = 3.4642, (2, 1) = 0., (2, 2) = 1.4142, (2, 3) = -.70710, (3, 1) = 0., (3, 2) = 0., (3, 3) = 1.2248})
MathCad Matrix(3, 3, {(1, 1) = 1.7321, (1, 2) = 3.4642, (1, 3) = 3.4642, (2, 1) = 0., (2, 2) = 1.4142, (2, 3) = -.70710, (3, 1) = 0., (3, 2) = 0., (3, 3) = -1.2248})
# Extract the lower-right submatrix:
M := SubMatrix(R, [2, 3], [2, 3]);
# Get
SingularValues(M, output = ['U', 'S', 'Vt']);
#
Error, (in LA_Main:-SingularValues) Matrix of type complex(extended_numeric) required for output object(s) [U, Vt]
