c=[2; -1; -3];

A = [ 1  -1  -1;
     -1   0   6;
      2   0   1];
b=[-2; -4; 6];
ci=[0; 0; 0];
k = 1;
[x,kl,f] = linpro(c, A, b, ci, [], k);
disp(x, "Значення невідомих x=");
disp(abs(f), "Максимальне значення функціі f= ");