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