Как найти значения s1 и s2, чтобы значения 1: 5000 и 5001: 10000 были разными?

Question

Это мой код:

  % Initial values
  x0 = 0.5;
  y0 = 0.6;
  dt = 0.01;
  time = 100; % time
  num_iteration = ceil(time / dt);
  num_trajectory = 10000; % Number of trajectories
  c = 1; % constant
  y = zeros(num_trajectory, num_iteration) + y0;
  x = zeros(num_trajectory, num_iteration) + x0;
  s1 = zeros(1, num_trajectory);
  s2 = zeros(num_trajectory);
  for t = 1:5000 % iteration 
     s1 = 2*x(:, t).^2./(y(:, t) + x(:, t) + c);
     s2 = 2*y(:, t);
  end
  for t1 = 5000:10000 % iteration
     s1 = x(:, t1)./y(:, t1);
     s2 = y(:, t1); 
  end

Мне не удалось найти s1 и s2 такие, что: для 1:5000, s1 и s2 будет иметь значение скажем, (0.7, 0.4) соответственно и от 5001:10000 s1 и s2 также будет иметь другое значение. В целом векторы s1 и s2 должны показывать некоторые изменения при изменении времени t.

Как мне найти такие значения?