diff --git a/matlab/utilities/general/compare_vectors.m b/matlab/utilities/general/compare_vectors.m
new file mode 100644
index 000000000..971007a78
--- /dev/null
+++ b/matlab/utilities/general/compare_vectors.m
@@ -0,0 +1,63 @@
+function C = compare_vectors(f, A, B)
+ 
+% Performs lexicographical comparison of vectors.
+%
+% INPUTS
+%  o f    function handle: @lt, @gt, @le or @ge
+%  o A    vector of real numbers.
+%  o B    vector of real numbers.
+%
+% OUTPUTS 
+%  o C    integer scalar, 1 or 0.
+%
+% REMARKS 
+%  o It is assumed that vectors A and B have the same number of elements.
+    
+% Copyright (C) 2013 Dynare Team
+%
+% This file is part of Dynare.
+%
+% Dynare is free software: you can redistribute it and/or modify
+% it under the terms of the GNU General Public License as published by
+% the Free Software Foundation, either version 3 of the License, or
+% (at your option) any later version.
+%
+% Dynare is distributed in the hope that it will be useful,
+% but WITHOUT ANY WARRANTY; without even the implied warranty of
+% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+% GNU General Public License for more details.
+%
+% You should have received a copy of the GNU General Public License
+% along with Dynare.  If not, see <http://www.gnu.org/licenses/>.
+
+if ~isvector(A) || ~isvector(B)
+    error('compare_vectors:: Input arguments a and b must be vectors!')
+end
+
+if ~isequal(length(A),length(B))
+    error('compare_vectors_lt:: Input arguments a and b must be of same length!')
+end
+
+if ~ismember(func2str(f),{'lt', 'gt', 'le', 'ge'})
+    error('compare_vectors:: First input argument must be one of the following function handles: @lt, @gt, @le or @ge!')
+end
+
+if feval(f, A(1), B(1))
+    C = 1;
+elseif isequal(A(1),A(1)) && length(A)>1
+    C = compare_vectors(f, A(2:end), B(2:end));
+elseif isequal(A(1),B(1)) && isequal(length(A),1) && ismember(func2str(f),{'le', 'ge'})
+    C = 1;
+else
+    C = 0;
+end
+
+%@test:1
+%$ t(1) = dyn_assert(compare_vectors(@lt, [1990 3], [1991 3]), 1);
+%$ t(2) = dyn_assert(compare_vectors(@lt, [1990 3], [1990 3]), 0);
+%$ t(3) = dyn_assert(compare_vectors(@le, [1990 3], [1990 3]), 1);
+%$ t(4) = dyn_assert(compare_vectors(@ge, [1991 3], [1990 3]), 1);
+%$ t(5) = dyn_assert(compare_vectors(@ge, [1991 3], [1991 3]), 1);
+%$ t(4) = dyn_assert(compare_vectors(@gt, [1991 3], [1990 3]), 1);
+%$ T = all(t);
+%@eof:1
\ No newline at end of file