function [] = PollingSystem()
clc

%Начальные данные
mu=[1 1 2   0;
    0 10 0   2;
    2 0 0.1 1;
    0 0 1   1];
lamb=[0.5; 0.5;  0.3; 0.4];
x_beg=[2;  1;  3; 3];
% mu=[0 1  2   0 0 2 0 2;
%     0 10 0   0 0 0 0 0;
%     2 0  0.1 0 0 3 0 0;
%     1 1  1   1 1 1 1 1;
%     0 2  0   2 0 0 0 2;
%     0 0  1   0 2 4 0 0;
%     0 0  0   4 5 0 0 0;
%     0 0  4   0 0 0 1 1;
%     9 0  0   0 2 0 0 0;  
%     0 1  0   0 0 9 0 0];
% lamb=[0.5; 0.5;  0.3; 0.4; 0.7; 0.5; 0.6; 0.7; 0.1; 0.3];
% x_beg=[2;  1;  3; 3; 5; 4; 7; 8; 2; 4];
period2=1;

period1=5;
Tmax=23;
step=0.01;

%Нахождение решений системы (2),(3) и системы (10)
tau=Solver(mu, lamb);
if isequal(tau, [])
    disp("Система неравенств решений не имеет. Система поллинга не стабилизируема.");
else
tau_bk=Work_time (mu, lamb, tau);


%Вычисление освобожденной и пришедшей работы на промежутках [0,(i-1)step] 
amount=Tmax/step+1;
free1=zeros(amount, 1);
in_flow=zeros(amount, 1);
free2=zeros(amount, 1);
for i=1:amount
    free1(i)=Free_work(mu, tau, tau_bk, period1, (i-1)*step);
    in_flow(i)=sum(x_beg)+(i-1)*step*(sum(lamb));
    free2(i)=Free_work(mu, tau, tau_bk, period2, (i-1)*step);
end

%Построение графика изменения суммарной очереди
set(findall(gcf,'type','axes'),'Fontsize',35);
set(findall(gcf,'type','text'),'FontSize',35) ;
t=0:step:Tmax;
Wip1=in_flow-free1;
plot(t, Wip1);
% Wip2=in_flow-free2;
%  plot(t, [Wip1, Wip2]);
% legend('period=5 сек','period=1 сек');

xlabel('Время (сек.)', 'Fontsize', 30);
ylabel('Суммарное содержимое', 'Fontsize', 35);
title('Система поллинга', 'Fontsize', 35);
handle_mas = findobj('type', 'axes');
set(handle_mas, 'FontSize', 35);
end
tau
end