
close all;
clear all;
x0 = [0 0 0];

%chaos
%params = [0.456 0.0357];
params = [0.6 0.0357];
tr = 30;

[t, phi] = ode45(@(t,x)blueSky(t,x, params),[0 500],x0);

phi = phi(t>tr,:);
t = t(t>tr);
plot3(phi(:,2),phi(:,3),phi(:,1));
xlabel('x');
ylabel('y');
zlabel('z');
grid on;
figure
plot(t,phi(:,1));
xlabel('t');
ylabel('x');