- matrix h80(80,80);
-
- for (unsigned r=0; r<80; ++r)
- for (unsigned c=0; c<80; ++c)
- h80.set(r,c,numeric(1,r+c+1));
- ex det = h80.determinant();
+ matrix hilbert(n,n);
+
+ for (unsigned r=0; r<n; ++r)
+ for (unsigned c=0; c<n; ++c)
+ hilbert.set(r,c,numeric(1,r+c+1));
+ ex det = hilbert.determinant();
+
+ /*
+ The closed form of the determinant of n x n Hilbert matrices is: