int main() { array double complex a[3][3] = {complex(3,4), 0, 6, 0, 2, 1, 1, 0, 1}; array double complex b[3] = {2, 13, 25}; array double complex x[3]; int status; clinsolve(x, a, b); printf("clinsolve(x, a,b) =\n%f\n", x); status = clinsolve(x, a, b); if(status == 0) printf("clinsolve(x,a,b) =\n%f\n", x); else printf("error: numerical error in clinsolve()\n"); }
clinsolve(x, a,b) = complex(17.760000,23.680000) complex(2.880000,11.840000) complex(7.240000,-23.680000) clinsolve(x,a,b) = complex(17.760000,23.680000) complex(2.880000,11.840000) complex(7.240000,-23.680000)