diff --git a/matlab/compute_model_moments.m b/matlab/compute_model_moments.m
index 5b09d82de..4cee165ef 100644
--- a/matlab/compute_model_moments.m
+++ b/matlab/compute_model_moments.m
@@ -1,4 +1,4 @@
-function moments=compute_model_moments(dr,M_,options,_)
+function moments=compute_model_moments(dr,M_,options_)
 %
 % INPUTS
 %    dr:        structure describing model solution