diff --git a/matlab/swz/c-code/sbvar/var/VARio.c b/matlab/swz/c-code/sbvar/var/VARio.c
index 1de0581cd..3e000128b 100644
--- a/matlab/swz/c-code/sbvar/var/VARio.c
+++ b/matlab/swz/c-code/sbvar/var/VARio.c
@@ -701,6 +701,8 @@ void WriteAllParameters(FILE *f, char *filename, char *id, TStateModel *model)
 
   WriteTransitionMatrices(f_in,(char*)NULL,id,model);
   Write_VAR_Parameters(f_in,(char*)NULL,id,model);
+
+  if(!f) fclose(f_in);
 }
 
 /*******************************************************************************/
diff --git a/matlab/swz/c-code/sbvar/var/mhm_VAR_main_1.c b/matlab/swz/c-code/sbvar/var/mhm_VAR_main_1.c
index 997cf4b11..c0defd12a 100644
--- a/matlab/swz/c-code/sbvar/var/mhm_VAR_main_1.c
+++ b/matlab/swz/c-code/sbvar/var/mhm_VAR_main_1.c
@@ -513,6 +513,7 @@ int main(int nargs, char **args)
       printf("Number inconsistent normalizations: %d\n",((T_VAR_Parameters*)mhm->model->theta)->WZ_inconsistancies);
       printf("Number singular inverse variances: %d\n\n",Get_VAR_Improper_Distribution_Counter());
 
+      fclose(mhm->f_out_regime_counts);
       fclose(mhm->f_out);
 
 /*        // After modified harmonic mean draws   ansi-c*/
diff --git a/matlab/swz/c-code/sbvar/var/probabilities.c b/matlab/swz/c-code/sbvar/var/probabilities.c
index cb7c7cbd9..17cdc287e 100644
--- a/matlab/swz/c-code/sbvar/var/probabilities.c
+++ b/matlab/swz/c-code/sbvar/var/probabilities.c
@@ -130,6 +130,7 @@ int main(int nargs, char **args)
       fprintf(f_out,"\n");
     }
 
+  fclose(f_out);
   swzFree(spec);
   swzFree(out);
   swzFree(head);