diff --git a/matlab/+bvar/irf.m b/matlab/+bvar/irf.m
index 05ec9b31e20e63de4ded3e71c7d2ec8e29b1789c..3889989059c68d65a935b8c7d950d71390c9d260 100644
--- a/matlab/+bvar/irf.m
+++ b/matlab/+bvar/irf.m
@@ -84,8 +84,8 @@ for draw=1:options_.bvar_replic
 
     % Build the IRFs...
     lags_data = zeros(ny,ny*nlags) ;
-    sampled_irfs(:,:,1,draw) = Sigma_lower_chol ;
-    lags_data(:,1:ny) = Sigma_lower_chol ;
+    sampled_irfs(:,:,1,draw) = StructuralMat ;
+    lags_data(:,1:ny) = StructuralMat ;
     for t=2:options_.irf
         sampled_irfs(:,:,t,draw) = lags_data(:,:)*Phi(1:ny*nlags,:) ;
         lags_data(:,ny+1:end) = lags_data(:,1:end-ny) ;