diff --git a/matlab/my_ordeig.m b/matlab/my_ordeig.m
index 9c2b989ea..5bc93fd0c 100644
--- a/matlab/my_ordeig.m
+++ b/matlab/my_ordeig.m
@@ -1,5 +1,20 @@
 function eval = my_ordeig(t)
-  
+
+% function eval = my_ordeig(t)
+% Computes the eigenvalues of a quasi-triangular matrix
+%
+% INPUTS
+%    t:              quasi-triangular matrix
+%
+% OUTPUTS
+%    eval:           eigenvalues
+%
+% SPECIAL REQUIREMENTS
+%    none
+%  
+% part of DYNARE, copyright Dynare Team (2003-2008)
+% Gnu Public License.
+
   n = size(t,2);
   eval = zeros(n,1);
   for i=1:n-1