ons. This script was used to
% produce Figure 2 in the article.
%
% 


%% Load the data and set font size for plots

data = load('lorenz_3d_timeseries.txt');
fontSize = 18;


%% Estimate time delay and plot AMI using all variables (Figure 2a in article)
tau = mdDelay(data, 'maxLag', 25, 'plottype', 'all');
set(gca,'FontSize',fontSize,'fontWeight','normal')
disp('xyz: tau = ' + string(tau))
print('Figure2a','-dpng')


%% Estimate the embedding dimension (Figure 2b in the article)
figure()
[fnnPercent, embedd