model selection_sw_comps-lin.mod; data simple.dat; option solver cplex; solve; display delivery.body; display exp(reliability.body); display {i in S, j in J[i] : x[i,j] = 1} sum{k in 1..Nmax} k*nu[i,j,k]; for {i in S, j in I[i] union J[i] : x[i,j] = 1} { printf "slot %d: component %d\n", i,j; }