diff --git a/vbr/vbrCore/functions/Q_xfit_mxw.m b/vbr/vbrCore/functions/Q_xfit_mxw.m
index 4c51b42..e895cfa 100644
--- a/vbr/vbrCore/functions/Q_xfit_mxw.m
+++ b/vbr/vbrCore/functions/Q_xfit_mxw.m
@@ -51,7 +51,6 @@
   f_norm_glob=J1; tau_norm_glob=J1;
 
   % vectorized rho and Vave
-  rho_vec = reshape(rho_in,size(Mu_in(1:n_th)));
   Vave=reshape(zeros(sz),size(Mu_in(1:n_th)));
 
   % ====================================================