% Hindmarsch-Rose neuron modell:
x0 = [0 0 0]; %kezdeti értékek

%params: tüzelési freki, vezérlõ áram, lassú változó ereje, adaptáció, nyugalmi pot, b
params = [1 3 0.01 4 -1.6 5]; %még nem burst
%params = [2 3 0.01 4 -1.6 5]; %fura burst kezdemény
%params = [2.5 3 0.01 4 -1.6 5]; %burst
%params = [3 3 0.01 4 -1.6 5]; %átmenet
%params = [3.5 3 0.01 4 -1.6 5]; %spike
%params = [4.5 3 0.01 4 -1.6 5]; %spike elhal
tr=0; %mekkora tranzienst hagyok le a kirajzolásnál az elejérõl



%{
%%spiking
params = [2.95 5 0.01 4 -1.6 5];
tr = 400;
%}

%{
%%burst
params = [2.6 2.66 0.01 4 -1.6 5];
tr = 400;
%}

%{
%%noidea
params = [1.87 5.24 0.01 4 -1.6 5];
tr = 400;
%}

%{
%%chaos ???
params = [2.96 2.99 0.01 4 -1.6 5];
tr = 100;
%}

%{
%blue sky
params = [2.63 5 0.01 4 -1.6 5];
tr = 200;
%params = [2.63 4.5 0.01 4 -1.6 5];
%tr = 200;
%}

figure;
%for a=1:0.2:6 %tüzelési frekv szabályzása, átkapcs spike és burst közt
  %  params(1)=a;
    % params(1)=3;
 %for I=2.5:0.01:2.8%2.615:0.0001:2.62% %szabálytalan viselkedés az input áram változtatásával
     %params(2)=I;
    
[t, phi] = ode45(@(t,x)hr(t,x, params),[0 1000],x0);

phi = phi(t>tr,:); % a kezdetre nem vagyok kíváncsi, csak a stabil állapotra
t = t(t>tr);
subplot(2,1,1);
plot3(phi(:,1),phi(:,2),phi(:,3));
xlabel('x: membrane potenciál');
ylabel('y: Na, K ion áramlás gyors csatornákon át: spike');
zlabel('z: egyéb ion áramlás (Ca) ,lassú csatornák: burst');
grid on;
subplot(2,1,2);
plot(t,phi(:,1));
xlabel('t');
ylabel('x: membrán potenciál');

title(num2str(params));
waitforbuttonpress;
%end