BC.positive.reg.RST function

Computation of a positive region