diff --git a/matlab/CheckPath.m b/matlab/CheckPath.m
new file mode 100644
index 0000000000000000000000000000000000000000..87ad27875cf80749216106a09d3e306e0ab97653
--- /dev/null
+++ b/matlab/CheckPath.m
@@ -0,0 +1,9 @@
+function DirectoryName = CheckPath(type)
+% 06-03-2005
+global M_
+
+DirectoryName = [ M_.dname '\' type ];
+
+if ~isdir(DirectoryName) 
+    system([ 'md ' DirectoryName]); 
+end
\ No newline at end of file
diff --git a/matlab/irf.m b/matlab/irf.m
index 12916b69b9ffbaaba0326737c0343e5f1e3f564e..53480efcaefd39141b597c661ce6fc2a862d8fdb 100644
--- a/matlab/irf.m
+++ b/matlab/irf.m
@@ -17,7 +17,7 @@ function y_=irf(dr, e1, long_, drop_, replic, iorder)
   
   else
     % eliminate shocks with 0 variance
-    i_exo_var = setdiff([1:M_.exo_nbr],find(diag(M_.Sigma_e) == 0));
+    i_exo_var = setdiff([1:M_.exo_nbr],find(diag(M_.Sigma_e) == 0 ));
     nxs = length(i_exo_var);
     ex1_ = zeros(long_+drop_+M_.maximum_lag,M_.exo_nbr);
     ex2_ = ex1_;
@@ -35,9 +35,7 @@ function y_=irf(dr, e1, long_, drop_, replic, iorder)
     y_=y_/replic;
   end
   options_.periods = old_iter;
-
 % 01/18/02 MJ corrected for many lags
-% 03/11/22 MJ input is now entire shock vector e1 (for orthogonalized
-%              IRFs)
+% 03/11/22 MJ input is now entire shock vector e1 (for orthogonalized IRFs)