proc iml;
reset storage=mlib.dea;
c = { 38 50 35 };
A = {  2  2  1, 0  2  3, 3  1  2, 0.5 0.5 0.5 };
b = { 132, 140, 144, 40 };
rel = j(4,1,'<=');
call mlinprog(c, 'max', A, rel, b, obj_val, primal, dual);
print obj_val primal dual;