--- a/liboctave/dMatrix.cc	Sat Sep 19 21:38:30 2009 +0200
+++ b/liboctave/dMatrix.cc	Sat Sep 19 22:17:12 2009 +0200
@@ -3414,13 +3414,13 @@
     retval(i, 0) = x1(i);
 
   // The last column is not needed while using delta.
-  double *delta = &retval(0, 1); 
+  double *delta = &retval(0, n-1); 
   for (octave_idx_type i = 0; i < m; i++)
     delta[i] = (x2(i) - x1(i)) / (n - 1);
 
   for (octave_idx_type j = 1; j < n-1; j++)
     for (octave_idx_type i = 0; i < m; i++)
-      retval(i, j) = retval(i, j-1) + delta[i];
+      retval(i, j) = x1(i) + j*delta[i];
 
   for (octave_idx_type i = 0; i < m; i++)
     retval(i, n-1) = x2(i);