int main() {
array double a[2][2] = {2, 4,
3, 7};
/* a2 is an ill-condition matrix */
array double a2[2][2] = {2, 4,
2.001, 4.0001};
array double complex z[2][2] = {2, 4,
3, 7};
array double complex z2[2][2] = {complex(2, 3), complex(4,2),
3, 7};
double rcond;
rcond = rcondnum(a);
printf("rcondnum(a) = %f\n", rcond);
rcond = rcondnum(a2);
printf("rcondnum(a2) = %f\n", rcond);
rcond = rcondnum(z);
printf("rcondnum(z) = %f\n", rcond);
rcond = rcondnum(z2);
printf("rcondnum(z2) = %f\n", rcond);
}
rcondnum(a) = 0.018182
rcondnum(a2) = 0.000079
rcondnum(z) = 0.018182
rcondnum(z2) = 0.131909