--- a/liboctave/dbleQRP.cc
+++ b/liboctave/dbleQRP.cc
@@ -53,11 +53,6 @@
 {
   assert (qr_type != QR::raw);
 
-  tau = 0;
-  work = 0;
-  tmp_data = 0;
-  jpvt = 0;
-
   int m = a.rows ();
   int n = a.cols ();
 
@@ -68,33 +63,32 @@
     }
 
   int min_mn = m < n ? m : n;
-  tau = new double[min_mn];
+  Array<double> tau (min_mn);
+  double *ptau = tau.fortran_vec ();
 
   int lwork = 3*n;
-  work = new double[lwork];
+  Array<double> work (lwork);
+  double *pwork = work.fortran_vec ();
 
   int info = 0;
 
+  Matrix A_fact;
   if (m > n)
     {
-      tmp_data = new double [m*m];
-      copy (tmp_data, a.data (), a.length ());
+      A_fact.resize (m, m);
+      A_fact.insert (a, 0, 0);
     }
   else
-    tmp_data = dup (a.data (), a.length ());
+    A_fact = a;
 
-  jpvt = new int[n];
-
-  // Clear Pivot vector (code to enforce a certain permutation would
-  // go here...)
+  double *tmp_data = A_fact.fortran_vec ();
 
-  for (int i = 0; i < n; i++)
-    jpvt[i] = 0;
+  Array<int> jpvt (n, 0);
+  int *pjpvt = jpvt.fortran_vec ();
 
-  F77_XFCN (dgeqpf, DGEQPF, (m, n, tmp_data, m, jpvt, tau, work, info));
+  // Code to enforce a certain permutation could go here...
 
-  delete [] work;
-  work = 0;
+  F77_XFCN (dgeqpf, DGEQPF, (m, n, tmp_data, m, pjpvt, ptau, pwork, info));
 
   if (f77_exception_encountered)
     (*current_liboctave_error_handler) ("unrecoverable error in dgeqpf");
@@ -107,18 +101,15 @@
 	{
 	  p.resize (1, n, 0.0);
 	  for (int j = 0; j < n; j++)
-	    p.elem (0, j) = jpvt[j];
+	    p.elem (0, j) = jpvt.elem (j);
 	}
       else
 	{
 	  p.resize (n, n, 0.0);
 	  for (int j = 0; j < n; j++)
-	    p.elem (jpvt[j]-1, j) = 1.0;
+	    p.elem (jpvt.elem (j) - 1, j) = 1.0;
 	}
 
-      delete [] jpvt;
-      jpvt = 0;
-
       volatile int n2;
 
       if (qr_type == QR::economy && m > n)
@@ -136,36 +127,24 @@
 	{
 	  int limit = j < min_mn-1 ? j : min_mn-1;
 	  for (int i = 0; i <= limit; i++)
-	    r.elem (i, j) = tmp_data[m*j+i];
+	    r.elem (i, j) = A_fact.elem (i, j);
 	}
 
       lwork = 32*m;
-      work = new double[lwork];
+      work.resize (lwork);
+      double *pwork = work.fortran_vec ();
 
-      F77_XFCN (dorgqr, DORGQR, (m, m, min_mn, tmp_data, m, tau, work,
-				 lwork, info));
+      F77_XFCN (dorgqr, DORGQR, (m, m, min_mn, tmp_data, m, ptau,
+				 pwork, lwork, info));
 
       if (f77_exception_encountered)
 	(*current_liboctave_error_handler) ("unrecoverable error in dorgqr");
       else
 	{
 	  q = Matrix (tmp_data, m, m);
-	  tmp_data = 0;
 	  q.resize (m, n2);
 	}
     }
-
-  delete [] tau;
-  tau = 0;
-
-  delete [] work;
-  work = 0;
-
-  delete [] tmp_data;
-  tmp_data = 0;
-
-  delete [] jpvt;
-  jpvt = 0;
 }
 
 /*