diff --git a/matlab/utilities/dataset/nanautocovariance.m b/matlab/utilities/dataset/nanautocovariance.m
index bc1cf582d17a0ed1f968e8ff94395876b43d023f..eb7a81ab2b54e7bf13672f51db0a37b565bb4441 100644
--- a/matlab/utilities/dataset/nanautocovariance.m
+++ b/matlab/utilities/dataset/nanautocovariance.m
@@ -55,8 +55,9 @@ function autocov = nanautocovariance(data,order)
 
 n = size(data,2);
 missing = isanynan(data);
+autocov = nan(n, n, order);
 order = min(size(data,1)-2,order);
-autocov = zeros(n, n, order);
+autocov(:, :, 1:order)=0;
 
 for lag=1:order
     if missing