# HG changeset patch
# User Kai T. Ohlhus <k.ohlhus@gmail.com>
# Date 1593587062 -32400
# Node ID 23a33db2bdb3f4ec68c7927ad6882deac4c9251a
# Parent  d4563c5d4060eeee77686bbda0b3ce638f337a74# Parent  627da618dcc4b7e84732c9c132c35ea5738070b3
maint: merge stable to default.

diff -r d4563c5d4060 -r 23a33db2bdb3 scripts/special-matrix/hilb.m
--- a/scripts/special-matrix/hilb.m	Mon Jun 29 23:45:19 2020 -0400
+++ b/scripts/special-matrix/hilb.m	Wed Jul 01 16:04:22 2020 +0900
@@ -66,12 +66,10 @@
     error ("hilb: N must be a scalar integer");
   endif
 
-  retval = zeros (n);
-  tmp = 1:n;
-  for i = 1:n
-    retval(i, :) = 1.0 ./ tmp;
-    tmp += 1;
-  endfor
+  ## Very elegant solution by N. Higham
+  ## https://nhigham.com/2020/06/30/what-is-the-hilbert-matrix/
+  j = 1:n;
+  retval = 1 ./ (j' + j - 1);
 
 endfunction