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