clc
clear
printf("------------task_1----------\n")
c=[1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0]';
A=[0 2.2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 -0.2 0 0 2.2 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 -0.2 0 0 0 -0.2 0 0 2.2 0 0 0 0 0 0 0 0 0 0;
   0 0 -0.2 0 0 0 -0.2 0 0 0 -0.2 0 0 2.2 0 0 0 0 0 0;
   0 0 -0.2 0 0 0 -0.2 0 0 0 -0.2 0 0 0 -0.2 0 0 2.2 0 0;
   0.588 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 -0.08 0.588 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 -0.08 0 0 0 -0.08 0.588 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 -0.08 0 0 0 -0.08 0 0 0 -0.08 0.588 0 0 0 0 0 0 0;
   0 0 0 -0.08 0 0 0 -0.08 0 0 0 -0.08 0 0 0 -0.08 0.588 0 0 0;
   1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 -2.2 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 -2.2 0 0 1 1 1 1 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 -2.2 0 0 1 1 1 1 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 -2.2 0 0 1 1 1 1;
   1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 -0.5 0 0 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 -0.5 0 0 1 0 0 1 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 -0.5 0 0 1 0 0 1 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 -0.5 0 0 1 0 0 1];
b=[5600 5600 5600 5600 5600 800 800 800 800 800 5500 0 0 0 0 2750 0 0 0 0]';
[x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1)
printf("------------task_2_1----------\n")
j=1;
k=1000
b(11)=0;
b(16)=0;
while(b(11)!=[10000])
  b(11)=b(11)+k;
  b(16)=b(16)+k/2;
  [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  y(j)=f_max;
  j=j+1;
 endwhile
disp(y);
x=1000:1000:10000;
plot(x,y);
title('Кол-во станков')
legend('boxoff');
legend('Кол-во стали, потраченное на производство станков',1)
xlabel('Изначальное кол-во станков')
ylabel('Потраченное кол-во стали')
grid on;
printf("------------task_2_2----------\n")
j=1;
k=1000
b(11)=0;
b(16)=0;
while(b(11)!=[10000])
  b(11)=b(11)+k;
  b(16)=b(16)+k/2;
  [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  y(j)=f_max;
  y1(j)=x_max(1)+x_max(5)+x_max(9)+x_max(13)+x_max(17);
  y2(j)=x_max(2)+x_max(6)+x_max(10)+x_max(14)+x_max(18);
  y3(j)=x_max(3)+x_max(7)+x_max(11)+x_max(15)+x_max(19);
  y4(j)=x_max(4)+x_max(8)+x_max(12)+x_max(16)+x_max(20);
  j=j+1;
 endwhile
 x=1000:1000:10000;
 figure()
 plot(x,y1,x,y2,x,y3,x,y4);
 title('Кол-во станков')
legend('boxoff');
legend('Кол-во стали, потраченное на производство станков','Кол-во стали, потраченное на производство стали','Кол-во стали, потраченное на улучшение сталелинейного завода','Кол-во стали, потраченное на улучшение станкостроительного завода',4)
xlabel('Изначальное кол-во станков')
ylabel('Потраченное кол-во стали')
grid on;
printf("------------task_3_1----------\n")
j=1;
k=100
b(6)=0;
b(7)=0;
b(8)=0;
b(9)=0;
b(10)=0;
while(b(6)!=[5000]);
  b(6)=b(6)+k;
  b(7)=b(7)+k;
  b(8)=b(8)+k;
  b(9)=b(9)+k;
  b(10)=b(10)+k;
  [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  y(j)=f_max;
  j=j+1;
 endwhile
 disp(y);
 x=100:100:5000;
 figure()
 plot(x,y);
 title('Кол-во станков')
legend('boxoff');
legend('Кол-во стали, потраченное на производство станков',1)
xlabel('Мощность производства')
ylabel('Кол-во стали')
grid on;
printf("------------task_3_1----------\n")
j=1;
k=100
b(6)=0;
b(7)=0;
b(8)=0;
b(9)=0;
b(10)=0;
while(b(6)!=[5000]);
  b(6)=b(6)+k;
  b(7)=b(7)+k;
  b(8)=b(8)+k;
  b(9)=b(9)+k;
  b(10)=b(10)+k;
  [x_max,f_max,en,extra]=glpk(c,A,b,[],[],"UUUUUUUUUUSSSSSUUUUU","CCCCCCCCCCCCCCCCCCCC",-1);
  y(j)=f_max;
  y1(j)=x_max(1)+x_max(5)+x_max(9)+x_max(13)+x_max(17);
  y2(j)=x_max(2)+x_max(6)+x_max(10)+x_max(14)+x_max(18);
  y3(j)=x_max(3)+x_max(7)+x_max(11)+x_max(15)+x_max(19);
  y4(j)=x_max(4)+x_max(8)+x_max(12)+x_max(16)+x_max(20);
  j=j+1;
 endwhile
 disp(y);
 x=100:100:5000;
 figure()
 plot(x,y1,x,y2,x,y3,x,y4);
 title('Кол-во станков')
legend('boxoff');
legend('Кол-во стали, потраченное на производство станков','Кол-во стали, потраченное на производство стали','Кол-во стали, потраченное на улучшение сталелинейного завода','Кол-во стали, потраченное на улучшение станкостроительного завода',4)
xlabel('Мощность производства')
ylabel('Кол-во стали')
grid on;