span = 1000;
[t, x, y] = rk45('blackbox', [0, span]);
subplot(2,1,1);
plot(t, x, '.');
xlabel('Time (second)');
title('Time Series');
drawnow

% Delete the initial unstable part
x(1:200) = [];

% Collect training data for your model
% ...
