diff --git a/src/lyapunov.c b/src/lyapunov.c
index 1cef444..b041a83 100644
--- a/src/lyapunov.c
+++ b/src/lyapunov.c
@@ -137,7 +137,8 @@ int lyapunov(
     }
 
     // QR decomposition
-    if((lyapunov_trigger==LYAPUNOV_TRIGGER_TIME && time>(n+1)*lyapunov_reset) || (lyapunov_trigger==LYAPUNOV_TRIGGER_SIZE && norm>lyapunov_reset)){
+    // Do it also if it is the last step
+    if((lyapunov_trigger==LYAPUNOV_TRIGGER_TIME && time>(n+1)*lyapunov_reset) || (lyapunov_trigger==LYAPUNOV_TRIGGER_SIZE && norm>lyapunov_reset) || time>=final_time){
       n++;
 
       // QR decomposition