diff --git a/.gitignore b/.gitignore
index 89bb900f537a44a3c57cd12d33d0f9766b8fe8df..1f51cd2e6ca1e884711e577058fe44142431b697 100644
--- a/.gitignore
+++ b/.gitignore
@@ -4,6 +4,7 @@
 *.eps
 
 # Folders to ignore
+figures/
 data/wtFiles/*/*
 data/arduinoFiles/*
 data/ORIGINAL_RAW_DATA
diff --git a/utils/comparefrequencies.m b/utils/comparefrequencies.m
new file mode 100644
index 0000000000000000000000000000000000000000..f67219768f650fcaa8d72b74e9a28ac1523dde35
--- /dev/null
+++ b/utils/comparefrequencies.m
@@ -0,0 +1,44 @@
+function [allRelDiff, meanDiff, stdDiff] = comparefrequencies(ResData)
+    % COMAPREFREQUENCIES Compare the frequencies of the two sets (should be identical)
+    %   Todo
+
+    % ----------------------------------------------------------------------------------------------
+    % (c) Copyright 2022 University of Liege
+    % Author: Thomas Lambert <t.lambert@uliege.be>
+    % ULiege - Aeroelasticity and Experimental Aerodynamics
+    % MIT License
+    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+    allRelDiff = zeros(size(ResData));
+
+    for i = 1:size(ResData, 1)
+        for j = 1:size(ResData, 2)
+
+            allFreqF(i, j) = ResData(i, j).Front.trueFreq;
+            allFreqA(i, j) = ResData(i, j).Aft.trueFreq;
+
+            allRelDiff(i, j) = abs(ResData(i, j).Front.trueFreq - ResData(i, j).Aft.trueFreq) / ...
+                ResData(i, j).Front.trueFreq;
+
+        end
+    end
+
+    allRelDiff = abs(allFreqF - allFreqA) ./ allFreqF;
+
+    meanDiff = mean(mean(allRelDiff));
+    stdDiff = mean(std(allRelDiff));
+
+    figure('Name', 'Frequency differences');
+    hold on;
+    plot(allFreqF, allFreqA, 'ok', 'MarkerFaceColor', 'k');
+    xlim([1.5, 3.5]);
+    ylim(xlim);
+    plot(xlim, xlim, '--');
+
+    hold off;
+    grid on;
+
+    xlabel('f_{aft}  [Hz]');
+    ylabel('f_{front}  [Hz]');
+
+end