function [z,x,w]=prodf(b)
	c=[5;7;4];
	A=[2,3,1; 1,1,3; 1,3,1];
	LB=[0,0,0];
	UB=[Inf, Inf, Inf];
	ct="UUU";
	vt="CCC";
	[x,z,status,extra]=glpk(c,A,b,LB,UB,ct, vt,-1);
	w=extra.lambda;
end;
	