Let's indicate with :
- $A(k)$ the number of adult cow at year $k$
- $C_0(k)$ the number of calf born at year $k$
- $C_1(k)$ the number of one-year-old calves at year $k$
- $C_2(k)$ the number of two-years-old calves at year $k$
- $C_3(k)$ the number of three-years-old calves at year $k$
Of course, the total population is $C(k) = A(k) + C_0(k)+C_1(k)+C_2(k)+C_3(k).$
The dynamics of the subpopulation is the following:
$$\begin{cases}
A(k+1) = A(k) + C_3(k)\\
C_0(k+1) = A(k)\\
C_1(k+1) = C_0(k)\\
C_2(k+1) = C_1(k)\\
C_3(k+1) = C_2(k)\\
\end{cases}$$
By iterating this, with initial conditions $A(0) = 1$, $C_0(0) = C_1(0) = C_2(0) = C_3(0) = 0$, you get that $C(17) = 185.$
Here you can find the Matlab code (if you are familiar with it) to iterate this:
clear all
close all
A = zeros(18, 1);
C0 = zeros(18, 1);
C1 = zeros(18, 1);
C2 = zeros(18, 1);
C3 = zeros(18, 1);
A(1) = 1;
for k=1:17
A(k+1) = A(k) + C3(k);
C0(k+1) = A(k);
C1(k+1) = C0(k);
C2(k+1) = C1(k);
C3(k+1) = C2(k);
end
C = A + C0 + C1 + C2 + C3;
figure
plot(0:17, A, 'r--', 'LineWidth', 2)
hold on
plot(0:17, C0, 'b--', 'LineWidth', 2)
plot(0:17, C1, 'g--', 'LineWidth', 2)
plot(0:17, C2, 'c--', 'LineWidth', 2)
plot(0:17, C3, 'm--', 'LineWidth', 2)
plot(0:17, C, 'k', 'LineWidth', 2)
legend('Adults', 'Brand new calves', 'One-year-old calves', 'Two-years-old calves', 'Three-years-old calves', 'Total cows')
xlim([0 17])
The (graphical) results you get is the following:
