%% skrypt do calkowania
%% wymaga pliku neur1.m
global Iinj


%% 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);
Iinj=0
[t,ZM]=ode45(@neur1,[0,100],[0 0 0 0]);
plot(t,ZM(1),'r-');
%end
%hold off;
Iinj=10;
[t,ZM]=ode45(@neur1,[0,100],[0 0 0 0]);
plot(t,ZM(1),'b-');
hold