function value = bevholt(t)
R = 2; K = 10; N0 = 1;
N(1) = N0;
for k = 1:t
    N(k+1) = R*N(k)/(1 + ((R-1)/K)*N(k));
end
value = N(t+1);