proc optmodel;
set Wezly={'z1','z2','mag','p1','p2','q1','q2','o1','o2','o3'};
set Luki={'z1p1','z1p2','z1mag','z2p1','z2p2','z2mag','p1q1','p2q2','q1o1','q1o2','q1o3','q2o1','q2o2','q2o3'};
num c{j in Luki}=[4 4 1 6 7 2 0 0 2 3 9 1 6 5];
num u{j in Luki}=[30 40 10 27 55 27 55 22 15 17];
num b{i in Wezly}=[30 40 10 0 0 0 0 -22 -15 -17];
num A{i in Wezly, j in Luki}=[1 1 1 0 0 0 0 0 0 0 0 0 0 0
                              0 0 0 1 1 1 0 0 0 0 0 0 0 0
							  0 0 -1 0 0 -1 0 0 0 0 0 0 0 0
							  -1 0 0 -1 0 0 1 0 0 0 0 0 0 0
							  0 -1 0 0 -1 0 0 1 0 0 0 0 0 0
							  0 0 0 0 0 0 -1 0 1 1 1 0 0 0
							  0 0 0 0 0 0 0 -1 0 0 0 1 1 1
							  0 0 0 0 0 0 0 0 -1 0 0 -1 0 0
							  0 0 0 0 0 0 0 0 0 -1 0 0 -1 0
							  0 0 0 0 0 0 0 0 0 0 -1 0 0 -1];
var x{j in Luki}>=0 <=u[j];
min KOSZT=sum{j in Luki}c[j]*x[j];
con war{i in Wezly}:sum{j in Luki}A[i,j]*x[j]=b[i];
solve;
print x;
quit;