%% Parc

clc
close
close all;
resting = [-1.2 -0.625];
%resting = [-3 -3];

N = 60;
ySS = zeros(1,N*2);
ySS(1:2:end) = resting(1);
ySS(2:2:end) = resting(2);
ySS(2) = -3;
ySS(1) = -3;

[t,x]=ode45('fnode_parc',[0 150],ySS);
%
figure
for i=1:N
    plot(t, x(:,(i-1)*2+1));
    ylim([-3.5 3.5])
    pause(0.1);
    %waitforbuttonpress;
end
%}
%
%figure
for i = 1:length(t)
   plot(linspace(1,N,200),interp1(1:N,x(i,1:2:2*N-1),linspace(1,N,200),'cubic'));
   ylim([-3.5 3.5])
   pause(0.005);
end
%}