Подскажите как продолжить решение, популяционная модель Николсона

День добрый, уважаемые, подскажите пожалуйста как решить задачку:

Начинаем с начальной популяции 100 и вводим 100 яиц одновременно (разово). Мы предполагаем, что яйца не были внесены до момента «0» (произвольное начало эксперимента).

Код:

clear all
close all


% Задержка вылупления яйца
a=16;  % дни

% Показатель выживаемости яиц
d=0.5; % как отношение день/день

% Показатель выживаемости куриц
b=0.7; % как отношение день/день

% Начальные условия
D0=100; % число куриц на "нулевой день"

% Время симулиции (всего)
c=49;  % в днях

%
D=zeros(c,1);

S=0*ones(c,1);
S(1,1)=100;

% Цикл 1.
D(1)=D0;% начальное количество популяции
for k=2:a
   D(k) = b*D(k-1);  
   disp("Current D(k):")
   disp(D(k-1))
end

% Цикл 2.
S(1)=D0; % начальное количество яиц
for k=16:a
   S(k) = d*S(k-16);
   disp("Current S(k):")
   disp(S(k-16))
end

% Теперь у нас есть переменные в рабочей области, и нам нужно отобразить результаты

figure,plot(D,S)

Преподаватель сказал, что до 1 цикла включительно - код правильный, со 2 цикла нет и дальше не уверен, в правильности кода. Подскажите, кто знает будьте любезны. Заранее спасибо;

Код приложил.