<![CDATA[
int main() {
array double v[4] = {2,3,4,5}; /* 2x^3 + 3x^2 + 4x + 5 */
int n = 4;
array double a[n-1][n-1];
a = companionmatrix(v);
printf("companionmatrix(a) =\n%f\n", a);
}
]]>
<![CDATA[
companionmatrix(a) =
-1.500000 -2.000000 -2.500000
1.000000 0.000000 0.000000
0.000000 1.000000 0.000000
]]>