')
disp('xyz: tau = ' + string(tau))


%% Time delay and plot only the x-variable
tau = mdDelay(data(:,1), 'maxLag', 25, 'plottype', 'all');
set(gca,'FontSize',fontSize,'fontWeight','normal')
disp('x: tau = ' + string(tau))

%% Time delay and plot only the y-variable
tau = mdDelay(data(:,2), 'maxLag', 25, 'plottype', 'all');
set(gca,'FontSize',fontSize,'fontWeight','normal')
disp('y: tau = ' + string(tau))


%% Time delay and plot only the z-variable
tau = mdDelay(data(:,3), 'maxLag', 25, 'plottype', 'all')