[Help-glpk] Out of domain error

[Juan Pablo Caballero Villalobos]

Hello for everybody, I have a problem with a little model I wrote in mathprog. 
When I try to execute it, I get the following error message:

Reading model section from planeacion_agregada.mod...
Reading data section from planeacion_agregada.mod...
50 lines were read
Generating costo...
planeacion_agregada.mod:24: Cw[1] out of domain
Model processing error


The complete model is this:

set T; /* Periodos de la planeacion */

param D{t in T};  /*demanda de productos en periodo t*/
param n{t in T};  /*numero de unidades que puede producir un operario en 
periodo t*/
param Cp{t in T}; /*costo unitario de produccion en periodo t*/
param Cw{t in T}; /*costo de un operario en periodo t*/
param Ch{t in T}; /*costo unitario de contratación en periodo t*/
param Cl{t in T}; /*costo unitario de despido en periodo t*/
param Ci{t in T}; /*costo unitario de inventario en periodo t*/
param Cb{t in T}; /*costo unitario de faltantes en periodo t*/
param i_inicial;  /*inventario inicial*/
param b_inicial;  /*faltantes inicial*/
param w_inicial;  /*número de operarios inicial*/


var p{t in T} >=0; /*unidades producidas en periodo t*/
var w{t in T} >=0; /*operarios disponibles en periodo t*/
var h{t in T} >=0; /*operarios contratados en periodo t*/
var l{t in T} >=0; /*operarios despedidos en periodo t*/
var i{t in T} >=0; /*unidades en inventario en periodo t*/
var b{t in T} >=0; /*unidades faltantes en periodo t*/


minimize costo: sum{t in T}(Cp[t]*p[t] + Cw[t]*w[t] + Ch[t]*h[t] +
                            Cl[t]*l[t] + Ci[t]*i[t] + Cb[t]*b[t]);


s.t. capacidad_prod{t in T}: p[t] <= n[t]*w[t];
s.t. balance_trabaj{t in T}: w[t] = if(t==1) then w_inicial else w[t-1] + h[t] 
- l[t];
s.t. invent_faltant{t in T}: i[t]-b[t] = if(t==1) then i_inicial else i[t-1] -
                                         if(t==1) then b_inicial else b[t-1] +
                                         p[t] - D[t];

data;

set T:= 1..6; /*Periodos*/

param D:= 1 2760 2 3320 3 3970 4 3540 5 3180 6 2900;  /*demanda de productos en 
periodo t*/
param n:= 1 84 2 80 3 92 4 84 5 88 6 88;  /*numero de unidades que puede 
producir un operario en periodo t*/
param Cp default 0; /*costo unitario de produccion en periodo t*/
param Cw:= 1 2520 2 2400 3 2760 4 2520 5 2640 6 2640; /*costo de un operario en 
periodo t, n*120 */
param Ch default 450; /*costo unitario de contratación en periodo t*/
param Cl default 600; /*costo unitario de despido en periodo t*/
param Ci default 5; /*costo unitario de inventario en periodo t*/
param Cb default 0; /*costo unitario de faltantes en periodo t*/
param i_inicial:= 0;  /*inventario inicial*/
param b_inicial:= 0;  /*faltantes inicial*/
param w_inicial:= 35;  /*número de operarios inicial*/

end;


I really will apreciate your colaboration to solve my model problem,



Ing. Juan Pablo Caballero Villalobos
Jefe Sección de Producción
Departamento de Procesos Productivos
Pontificia Universidad Javeriana
320 8320 ext 5362
address@hidden