%
global a
b = [3 20 100 500];
c = [2 0.1 1 1];

for i=1:4
  a = [-b(i) -c(i); 1 0];
  lambda = eig(a);
  eigenratio(i) = max(abs(lambda))/min(abs(lambda));
  time0 = clock;
  inity = [0 1]';
  %[t, y] = ode45('f508', [0,2], inity);
  [t, y] = ode23s(@f508, [0,2], inity, 0.005);
  et(i) = etime(clock, time0);
%  pause
end;
plot(t, y, 'o')
eigenratio
et
