%%Hodgkin-Huxley: AP terjedés:
close all;
resting = [0.3177    0.0529    0.5962  -59.9995]; %nyugalmi értékek
firstpot = -20;%-47; %kezdeti potenciálváltozás axon elején

N = 60;
%az axon minden pontjához a vektor 4 mezõjét rendelem (n, m, h, V), 
%tehát egy pontot egy 4 hosszú vektor jelöl, 
%ezért léptetek az ySS mentén 4-esével
ySS = zeros(1,N*4);
ySS(1:4:end) = resting(1);
ySS(2:4:end) = resting(2);
ySS(3:4:end) = resting(3);
ySS(4:4:end) = resting(4);
ySS(4) = firstpot;

[t,x]=ode45('hhode_parc',[0, 60],ySS);
%
figure;
for i = 1:length(t)
   plot(linspace(1,N,200),  interp1(x(i,4:4:4*N),linspace(1,N,200),'cubic')  ); %a hhode_parc-ban kiszámolt értékeket (x) finomítom a linspace-szel megadott pontossághoz
  %200 helyett pl 20 ponton interpolálok-> darabosabb fv, a pulzálás a
  %numerikus pontatlanság miatt van
   xlim([0, N]);
  ylim([-80, 80]);
  title(['Time: ', num2str(t(i))]);
  xlabel('Axon');
  ylabel('Activity');
  pause(0.005);
%   waitforbuttonpress;
end
%}