clear;close all;clc;

lambda = .5;
Ts = 2;
Tf = 10;

N = 10;

jobs = - log(rand(1, N)) / lambda;
job_durations = mod(rand(1, N), Tf) + Ts;

T0 = 0;

t = [];
t_start = [];
t_end = [];
for c = 1:N
    if c == 1
        t(c) = T0;
        t_start(c) = T0;
        t_end(c) = T0 + mod(rand(), Tf) + Ts;
    else
        t(c) = t(c - 1) + jobs(c);
        t_start(c) = max(t(c), t_end(c - 1));
        t_end(c) = t_start(c) + mod(rand(), Tf) + Ts;
    end
end

fprintf('Cреднее время нахождения заявки в системе: %f\n', mean(t_end - t_start));
fprintf('Cреднее время ожидания в очереди: %f\n', mean(t_start - t));
fprintf('Cреднее время простоя обслуживающего прибора: %f\n', mean(t_start(2:end) - t_end(1:end - 1)));

queue = [];
lock = false;

% fprintf('Средняя длина очереди: %f\n', );
% fprintf('Максимальная длина очереди: %f\n', );