<![CDATA[
int main() {
array double a[3][3] = {3, 0, 6,
0, 2, 1,
1, 0, 1};
array double b[3] = {2,
13,
25};
array float a1[3][3] = {3, 0, 6,
0, 2, 1,
1, 0, 1};
array float b1[3] = {2,
13,
25};
array double x[3];
int status;
linsolve(x, a, b);
printf("linsolve(a,b) =\n%f\n", x);
linsolve(x, a1, b1);
printf("linsolve(a1,b1) =\n%f\n", x);
status = linsolve(x, a, b);
if(status == 0)
printf("linsolve(x,a,b) =\n%f\n", x);
else
printf("error: numerical error in linsolve()\n");
}
]]>
<![CDATA[
linsolve(a,b) =
49.333333 18.666667 -24.333333
linsolve(a1,b1) =
49.333333 18.666667 -24.333333
linsolve(x,a,b) =
49.333333 18.666667 -24.333333
]]>