diff --git a/utils/calctrueforces.m b/utils/calctrueforces.m
index ced4147600e561cba348947448e8caef61303ce3..9d52b9ae365e542492647d8e06122fe8d76514d4 100644
--- a/utils/calctrueforces.m
+++ b/utils/calctrueforces.m
@@ -23,13 +23,14 @@ function [ResData] = calctrueforces(ResData)
                         inF = ResData(i, j).AllPhases(1).(pos).inertiaF;
                         inM = ResData(i, j).AllPhases(1).(pos).inertiaM;
                     else
-                        inF = 0;
-                        inM = 0;
+                        inF = [0, 0, 0];
+                        inM = [0, 0, 0];
                     end
 
                     % Raw force
-                    meanFraw = mean(ResData(i, j).(pos).forces(indexes(iPart, :), :));
-                    meanMraw = mean(ResData(i, j).(pos).moments(indexes(iPart, :), :));
+                    idx = indexes(iPart, 1):indexes(iPart, 2);
+                    meanFraw = mean(ResData(i, j).(pos).forces(idx, :));
+                    meanMraw = mean(ResData(i, j).(pos).moments(idx, :));
 
                     % True value
                     ResData(i, j).AllPhases(1).(pos).meanF(iPart, :) = meanFraw - inF;