Page 1 of 1

### Equation

Posted: Fri Jul 26, 2019 11:35 pm
In the attachment, 1) constraint 12 has variable x(i,t) but how write that constraint? Do I have to specify three variables x, x1, x0? where x is amount produced.

2) constraint 15, How to write last summation part which is multiplied by M? (k=t+1) to (s-1)

Thank you.

### Re: Equation

Posted: Sat Jul 27, 2019 9:09 am
Hi Omi

1) Only X^0 seems to be an additional variable. X^(-1) is just 1/X.
2) just multiply it with the sum?

Cheers
Renger
BTW: You should start submitting your own code by now, so people can answer your questions more easily building on what you submit.

### Re: Equation

Posted: Sat Jul 27, 2019 3:17 pm
Hello Renger,
Thank you for your response.
This is the code I tried to write for the above constraints 12 to 15.
is this correct?
x is a production amount for product j and I want to split this amount.

Code: Select all

``````ct3(j,t)..x(j,t)=e=x1(j,t)+x0(j,t);

ct4(j,t)..x1(j,t)=l=M*z(j,t);

ct5(j,t)..x0(j,t)=g=mi*(sum((jj),U(jj,j,t)-z(j,t+1)));

ct6(j,t,s)..x0(j,t)+sum(k\$(ord(k)=ord(t)+1),x1(j,k))=g=mi*sum(jj,U(jj,j,t))-M*(sum(k\$(ord(k)=ord(t)+1),R(k))+1-R(s));``````

### Re: Equation

Posted: Mon Jul 29, 2019 8:39 am
Hi Omi

It looks fine to me. If you not sure, if the equation does the right thing, you could initialize the variables in the equations to some values, and define checking parameters to see what they produce.
Here a small example:

Code: Select all

``````set i /1*4/;
variable X(i), OBJ;

equation c1;
c1.. OBJ =E= sum(i, X(i));

* Checking the equations
X.L(i) = 1;

parameter check;
check = sum(i, X.L(i));
display check;
``````
Check should now have 4 as value. You could do the same thing with the equations of your model.

Cheers

Renger

BTW: If you sent code, make sure it runs. In order to check your equations by running the code, I had to make some assumptions and so it is hard to say if the equations are really correct.
I had to add the following to be able to run the code:

Code: Select all

``````set i /1*4/,  t /1*4/,  j /1*2/;

alias(k,s,t);

variables r(k), x(j,t), x1(j,t), x0(j,t), z(j,t);

parameter M /1000/, mi /199/;

parameter u(j,j,t);

alias(j,jj);

equations ct3, ct4, ct5,ct6;
``````