<![CDATA[
int main() {
array double a[4][3] = {1,2,3,
4,5,6,
7,8,9,
6,3,5};
int n = 4, m = 3, k;
array double t[n][m];
t = triangularmatrix("upper",a);
printf("triangularmatrix(\"upper\", t) =\n%f\n", t);
k = 0;
t = triangularmatrix("upper",a,k);
printf("triangularmatrix(\"upper\", t, 0) =\n%f\n", t);
k = 1;
t = triangularmatrix("upper",a,k);
printf("triangularmatrix(\"upper\", t, 1) =\n%f\n", t);
k = -1;
t = triangularmatrix("upper",a,k);
printf("triangularmatrix(\"upper\", t, -1) =\n%f\n", t);
t = triangularmatrix("lower",a);
printf("triangularmatrix(\"lower\", t) =\n%f\n", t);
k = 0;
t = triangularmatrix("lower",a,k);
printf("triangularmatrix(\"lower\", t, 0) =\n%f\n", t);
k = 1;
t = triangularmatrix("lower",a, k);
printf("triangularmatrix(\"lower\", t, 1) =\n%f\n", t);
k = -1;
t = triangularmatrix("lower",a,k);
printf("triangularmatrix(\"lower\", t, -1) =\n%f\n", t);
}
]]>
<![CDATA[
triangularmatrix("upper", t) =
1.000000 2.000000 3.000000
0.000000 5.000000 6.000000
0.000000 0.000000 9.000000
0.000000 0.000000 0.000000
triangularmatrix("upper", t, 0) =
1.000000 2.000000 3.000000
0.000000 5.000000 6.000000
0.000000 0.000000 9.000000
0.000000 0.000000 0.000000
triangularmatrix("upper", t, 1) =
0.000000 2.000000 3.000000
0.000000 0.000000 6.000000
0.000000 0.000000 0.000000
0.000000 0.000000 0.000000
triangularmatrix("upper", t, -1) =
1.000000 2.000000 3.000000
4.000000 5.000000 6.000000
0.000000 8.000000 9.000000
0.000000 0.000000 5.000000
triangularmatrix("lower", t) =
1.000000 0.000000 0.000000
4.000000 5.000000 0.000000
7.000000 8.000000 9.000000
6.000000 3.000000 5.000000
triangularmatrix("lower", t, 0) =
1.000000 0.000000 0.000000
4.000000 5.000000 0.000000
7.000000 8.000000 9.000000
6.000000 3.000000 5.000000
triangularmatrix("lower", t, 1) =
1.000000 2.000000 0.000000
4.000000 5.000000 6.000000
7.000000 8.000000 9.000000
6.000000 3.000000 5.000000
triangularmatrix("lower", t, -1) =
0.000000 0.000000 0.000000
4.000000 0.000000 0.000000
7.000000 8.000000 0.000000
6.000000 3.000000 5.000000
]]>