OccBin simulation problem

Occbin
2

Hi,
your constraint is wrong.

Instead of

constraint = 'i<1-1/beta';
constraint_relax ='i>1-1/beta';

it is better to use the shadow rate within constraints. E.g.

constraint = 'inot<1-1/beta';
constraint_relax ='inot>1-1/beta';