x=zeros(3,4);
y=x;
C=4;
CO=3;
vi=[480 650 580 390];
p=[310; 380; 350; 285];

b1=[18 15 23 12];
b3=[6800 8700 5300 0];
b2=[10 16 8 0];
for w=1:1:C
    a1(w)=sum(y(:,w));
end
for r=1:1:CO
    a2(r)=sum(y(r,:));
    a2(1,4)=0;
end
for h=1:1:CO
    a3(h)=sum(vi(h).*y(h,:));
    a3(1,4)=0;
end
for u=1:1:CO
    a4(u)=1./b2(u)*sum(y(u,:));
end
for o=1:1:C
    for j=1:1:CO
fobj(j,o) = p(o)*y(j,o);
    end
end
f=fobj;
A=[a1';a2';a3'];
b=[b1;b2;b3];
Aeq=a4;
beq=z;
lb = zeros(CO,C);
ub=[];
[x,fvar,exitflag,output,lambda] = linprog(f,A,b,Aeq,beq,lb,ub);