int main() {
array double a[2][2] = {2, 4,
3, 7};
/* a2 is an ill-conditioned matrix */
array double a2[2][2] = {2, 4,
2.001, 4.0001};
array float b[2][2] = {2, 4,
3, 7};
array double complex z[2][2] = {2, complex(4, 3),
3, 7};
/* z2 is an ill-conditioned matrix */
array double complex z2[2][2] = {2, complex(4, 3),
2.001, complex(4.0001,3)};
array complex z3[2][2] = {2, complex(4, 3),
3, 7};
double cond;
cond = condnum(a);
printf("condnum(a) = %g\n", cond);
cond = condnum(a2);
printf("condnum(a2) = %g\n", cond);
cond = condnum(b);
printf("condnum(b) = %g\n", cond);
cond = condnum(z);
printf("condnum(z) = %g\n", cond);
cond = condnum(z2);
printf("condnum(z2) = %g\n", cond);
cond = condnum(z3);
printf("condnum(z3) = %g\n", cond);
}
condnum(a) = 38.9743
condnum(a2) = 10527.6
condnum(b) = 38.9743
condnum(z) = 9.32929
condnum(z2) = 11980.8
condnum(z3) = 9.32929