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

%%spiking
%params = [2.95 5 0.01 4];
%tr = 400;

%%burst
%params = [2.6 2.66 0.01 4];
%tr = 400;

%%noidea
%params = [1.87 5.24 0.01 4];
%tr = 400;

%%chaos
%params = [2.96 2.99 0.01 4];
%tr = 100;

%blue sky
%params = [2.63 5 0.01 4];
%tr = 200;
%params = [2.63 4.5 0.01 4];
%tr = 200;

[t, phi] = ode45(@(t,x)hr(t,x, params),[0 3000],x0);

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