# HG changeset patch
# User jwe
# Date 775160563 0
# Node ID b805b12f3b6616bcb432653e3ace869b6e3beba3
# Parent  50cf0f8ebde13b1a2a127ffba4db3d151dbc8f34
[project @ 1994-07-25 18:22:43 by jwe]

diff --git a/scripts/special-matrix/hilb.m b/scripts/special-matrix/hilb.m
--- a/scripts/special-matrix/hilb.m
+++ b/scripts/special-matrix/hilb.m
@@ -1,4 +1,4 @@
-# Copyright (C) 1993 John W. Eaton
+# Copyright (C) 1993, 1994 John W. Eaton
 # 
 # This file is part of Octave.
 # 
@@ -35,10 +35,9 @@
   nmax = length (n);
   if (nmax == 1)
     retval = zeros (n);
-    for j = 1:n
-      for i = 1:n
-        retval (i, j) = 1 / (i + j - 1);
-      endfor
+    tmp = 1:n;
+    for i = 1:n
+      retval (i, :) = 1.0 ./ (tmp + (i - 1));
     endfor
   else
     error ("hilb: expecting scalar argument, found something else");