%% skrypt do calkowania
%% wymaga pliku neur2.m
global I Transm


%% delta =0.008

%% epsilon = 0.1
clf
hold on
%for i=1:10
%    x=-0.01;%rand();
%    y=0.07;%rand();
%% runge-kutty metoda calkowania funkcji polowy3 od t =0 do t=20 
%options = odeset('RelTol',1e-15,'AbsTol',1e-15,'MaxStep',1e-3);
I=0
Transm=0;
[t,ZM]=ode15s(@neur2,[0,2000],[0 0 0 0 0 0 0 0 0 0]);
%I=-0.2e-3;
Transm=1;
[t1,ZM1]=ode15s(@neur2,[2000,2001],ZM(end,:));
%I=0
Transm=0;
[t2,ZM2]=ode15s(@neur2,[2001,4000],ZM1(end,:));
subplot(2,1,1)
plot([t; t1; t2],[ZM(:,1); ZM1(:,1); ZM2(:,1)] ,'r-');
subplot(2,1,2)
plot([t; t1; t2;],[ZM(:,2);ZM1(:,2);ZM2(:,2)],'r-');
%end
%hold off;
%Iinj=0;
%[t,ZM]=ode113(@neur1,[0,100],[0 0 0 0]);
%I=10;
%[t,ZM]=ode45(@neur2,[0,2000],[0 0 0 0 0 0 0 0 0 0]);%ZM(end,:));
%plot(t,ZM(:,1),'b-');
%hold off;