Pesquisa operacional

sets i produtos/1*4/ t mes /1*4/;

alias (t,j);

parameters c(i) capacidade de producao do item i
/
1 8000
2 5000
3 8000
4 15000/

s(i) custo de preparacao
/
1 12000
2 28000
3 13800
4 19500
/

h(i) custo de estocagem
/
1 35
2 39
3 45
4 85
/
g(i) custo de falta
/
1 30
2 30
3 40
4 80
/;

table d(i,t) demanda 1 2 3 4
1 5000 6000 3000 1000
2 900 1000 4000 5000
3 6000 9000 4000 2000
4 10000 11000 14000 16000;

scalar k capacidade total /28000/;

variables x(i,t) quanto produzir do item i no perido t z fo e(i,t) quanto estocar do produto i em t y(i,t) setup do produto i em t f(i,t) falta;

binary variable y; positive variable e,x,f;

equations custo custo dem(i,t) balancemanto de demanda cap1(i,t) capacidade do item cap2(t) capacidade total lim(i,t) limite de producao;

custo.. z=e= sum((i,t), s(i)*y(i,t) + h(i)*e(i,t) +g(i)*f(i,t)) ; dem(i,t).. e(i,t) - f(i,t)=e= +f(i,t-1) + e(i,t-1) + x(i,t) -d(i,t); cap1(i,t).. x(i,t)=l= c(i); cap2(t).. sum(i, x(i,t))=l= k; lim(i,t).. x(i,t)=l= sum(j,d(i,j)*y(i,t));

model cap /all/

solve cap using mip minimizing z;

display