int main() {
array double complex a[2][2] = {complex(1,2), 4,
3, 7};
/* a2 is an ill-detition matrix */
array double a2[2][2] = {2, 4,
2.001, 4.0001};
/* a3 is singular */
array double a3[2][2] = {2, 4,
4, 8};
array complex b[2][2] = {2, 4,
3, 7};
array double c[3][3] = {-1,5,6,
3,-6,1,
6,8, 9} ; /* n-by-n matrix */
double complex det;
det = cdeterminant(a);
printf("cdeterminant(a) = %g\n", det);
det = cdeterminant(a2);
printf("cdeterminant(a2) = %g\n", det);
det = cdeterminant(a3);
printf("cdeterminant(a3) = %g\n", det);
det = cdeterminant(b);
printf("cdeterminant(b) = %g\n", det);
det = cdeterminant(c);
printf("cdeterminant(c) = %g\n", det);
}
cdeterminant(a) = complex(-5
,14
)cdeterminant(a2) = complex(-0.0038
,-0
)cdeterminant(a3) = complex(-0
,-0
)cdeterminant(b) = complex(-2
,0
)cdeterminant(c) = complex(317
,-0
)