<![CDATA[
#define N1 3
#define N2 4
int main() {
array double a[N1][3] = {3, 0, 6,
0, 2, 1,
1, 0, 1};
array double a2[N2][3] = {3, 0, 6,
0, 2, 1,
1, 0, 1,
4, 5, 2};
array double b[N1] = {2,
13,
25};
array double b2[N2] = {2,
13,
25,
1};
array double x[N1], x2[N2];
int status;
llsqsolve(x, a, b);
printf("llsqsolve(x,a,b) = \n%f\n", x);
llsqsolve(x2, a2, b2);
printf("llsqsolve(x2, a2,b2) = \n%f\n", x2);
status = llsqsolve(x, a, b);
if(status == 0)
printf("llsqsolve(x, a,b) =\n%f\n", x);
else
printf("error: numerical error in llsqsolve()\n");
}
]]>
<![CDATA[
llsqsolve(x,a,b) =
49.333333 18.666667 -24.333333
llsqsolve(x2, a2,b2) =
-2.256881 1.715596 2.198777 -26.654661
llsqsolve(x, a,b) =
49.333333 18.666667 -24.333333
]]>