diff --git a/src/MLIAP/mliap_so3.cpp b/src/MLIAP/mliap_so3.cpp
index 2027119d86..97e6abb010 100644
--- a/src/MLIAP/mliap_so3.cpp
+++ b/src/MLIAP/mliap_so3.cpp
@@ -477,8 +477,14 @@ void MLIAP_SO3::compute_W(int nmax, double *arr)
   delete sqrtD;
   delete tempM;
 
-  delete temparr;
-  delete tempvl;
+  for (int iy=0; iy<n; iy++)
+    delete[] temparr[iy];
+  delete[] temparr;
+
+  for (int iy=0; iy<n; iy++)
+    delete[] tempvl[iy];
+  delete[] tempvl;
+
 
 }
 
diff --git a/src/MLIAP/mliap_so3_math.h b/src/MLIAP/mliap_so3_math.h
index 556ba8570f..54f75da1ed 100644
--- a/src/MLIAP/mliap_so3_math.h
+++ b/src/MLIAP/mliap_so3_math.h
@@ -25,7 +25,11 @@ typedef Jacobi<double,double*,double**,double const*const*> Jacobi_v2;
 int SO3Math::jacobin(int n, double const* const* mat, double *eval, double **evec)
 {
   int i,j;
-  double mat_cpy[n][n];
+  double **mat_cpy;
+
+  mat_cpy = new double* [n];
+  for (int iy=0; iy<n; iy++)
+    mat_cpy[iy] = new double[n];
 
   for(i=0;i<n;i++)
     for(j=0;j<n;j++)
@@ -44,6 +48,10 @@ int SO3Math::jacobin(int n, double const* const* mat, double *eval, double **eve
     for (int j=i+1; j<n; j++)
       std::swap(evec[i][j], evec[j][i]);
 
+  for (int iy=0; iy<n; iy++)
+    delete[] mat_cpy[iy];
+  delete[] mat_cpy;
+
   return ierror;
 }