function vecAmount = simulateBacteriaAmount(vecTime,initAmount) % Find the number of elements in vecTime lengthTime = length(vecTime); % Initialise the vector vecAmount vecAmount = zeros(lengthTime,1); % Amount at vecTime(1) vecAmount(1) = initAmount; % Iteration to find the amount at other time instances for ii = 2:lengthTime vecAmount(ii) = vecAmount(ii-1) + (vecTime(ii)-vecTime(ii-1)) * rate(vecTime(ii-1)); end