close all;

N = 100;
ySS = zeros(1,N);
ySS(1) = 1;

[t,x]=ode45('fischer',[0 50],ySS);
%{
figure
for i=1:N
    plot(t, x(:,(i-1)*2+1));
    ylim([-3.5 3.5])
    pause(0.1);
    %waitforbuttonpress;
end
%}
%{
 %if t>1
   kul=abs(x-0.5);
mi=min(kul);
ind= find(kul<mi+0.000000000000001);
indv(nt,1)=ind;
  %  end
%}
%
figure
for i = 1:length(t)
   plot(1:N,x(i,:));
   ylim([-1 2])
   pause(0.01);
end
%}