Note, however, that your equation is defined over J and K and then you try to use J again for summing. THis won't work (how will GAMS keep track of which j to use for the equation itself and the sum?).

The equation makes more sense as being defined only over the set K (EMISSION(K)).

i had tried using K (EMISSION(K)) before ,but more error appeared
8 ')' expected
10 ',' expected
119 Number (primary) expected
121 Set expected
171 Domain violation for set
257 Solve statement not checked because of previous errors

Maybe i should say i propose of those,
i would mutiply the p(j.k) with a coeifficent f which less than a limit value,
is it possible to do that? Or i need to review my model?