diff --git a/matlab/var_sample_moments.m b/matlab/var_sample_moments.m
index 834eb81522a1945f0cce60ab5c054f52308a0d2c..ab7b4b74b6575b2be7d6d04ba19593ae560b3495 100644
--- a/matlab/var_sample_moments.m
+++ b/matlab/var_sample_moments.m
@@ -78,7 +78,7 @@ XtX = [];
 data = read_variables(datafile,varobs,[],xls_sheet,xls_range);
 
 if qlag > FirstObservation
-    disp('VarSampleMoments :: not enough data to initialize! Try to increase FirstObservation.')
+    error('VarSampleMoments :: not enough data to initialize! Try to increase FirstObservation.')
     return
 end
 
@@ -93,7 +93,7 @@ elseif var_trend_order == 1;% Constant and linear trend case.
 X = ones(NumberOfObservations,NumberOfVariables*qlag+2);
 indx = NumberOfVariables*qlag+1:NumberOfVariables*qlag+2;
 else
-    disp('var_sample_moments :: trend must be equal to -1,0 or 1!')
+    error('var_sample_moments :: trend must be equal to -1,0 or 1!')
     return
 end