diff --git a/dynare++/kord/first_order.cc b/dynare++/kord/first_order.cc
index 443484f58a04c8552da191ca0adabae28aef29d0..7986481337f058d91f1fb3de205d56788cce3881 100644
--- a/dynare++/kord/first_order.cc
+++ b/dynare++/kord/first_order.cc
@@ -14,7 +14,7 @@ double qz_criterium = 1.000001;
 lapack_int
 order_eigs(const double *alphar, const double *alphai, const double *beta)
 {
-  return (*alphar **alphar + *alphai **alphai < *beta **beta * qz_criterium);
+  return (*alphar **alphar + *alphai **alphai < *beta **beta * qz_criterium * qz_criterium);
 }
 
 /* Here we solve the linear approximation. The result are the matrices