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)