clear

m=3;                                                % number of locations
n=5;                                                % number of customers

s=[239 225 184];                                    % capacity
d=[92 82 83 69 74];                                 % demand
f=[589 766 886];                                    % fixed cost
c=[14 5 6 24 6; 9 22 26 5 21; 16 11 23 28 24];      % variable cost

x=optimvar('x',[m,n],'Upper',0);
y=optimvar('y',m,'type','integer','Lower',0,'Upper',1);

T=triu(c);
objective=sum(f(:).*y(:))+sum(T(:).*x(:));

constraints.colsum=sum(x,1)==d(:).';
constraints.rowsum=sum(x,2)<=s(:).*y(:);

p=prob2struct(optimproblem('Objective', objective,'Constraints',constraints))