macro triangle [1]
vec/dele * 
vec/cre babar(1)
vec/inp babar 0. 
vec/cre rho(100)
vec/cre eta(100)
vec/cre rho_amb(4)
vec/cre eta_amb(4)
vec/cre eta_up
vec/inp eta_up 1.

*----  Initialisation 
call contrainte2.for(0.)
zone
set ymgu 5
set ymgl 7
null -1. 1. 0 1.



*----  Dmd
call contrainte2.for(1.)
igset fais 3
igset faci 2
igset plci 1
igset fasi 345
farea 100 rho eta
pline 100 rho eta
igset plci 2 
if [1] = 'old' goto jump1 
if [1] = 'OLD' goto jump1 
dbox -1. -.8 1.5 1.6
text -.75 1.5 '[D]m?d! constraint' .3
jump1:
*----  Epsilon K 
call contrainte2.for(2.)
igset fais 3
igset faci 3
igset plci 1
igset fasi 395
farea 100 rho eta
pline 100 rho eta
if [1] = 'old' goto jump2 
if [1] = 'OLD' goto jump2 
dbox -1. -.8 1.3 1.4
text -.75 1.3 '[e]?K! constraint' .3
jump2:
*----  Vbu/Vcb
call contrainte2.for(3.)
igset fais 3
igset faci 7
igset plci 1 
igset fasi 353
farea 100 rho eta
pline 100 rho eta
if [1] = 'old' goto jump3 
if [1] = 'OLD' goto jump3 
dbox -1. -.8 1.1 1.2
text -.75 1.1 'V?ub!/V?cb! constraint' .3
jump3:
if $sigma(babar(1)) = 0 goto nobabar

2dh 2122 'chi2' 100 -1. 1. 100 0. 1.

call contrainte2.for(4.)
igset plci 4 
set hcol 4 
vec/cre cl(1) r 6. 
contour 2122 1 s cl 	

igset pmci 4 
set mtyp 20
igset mscf 1.5
pmarker 4 rho_amb eta_amb
igset mscf 1.
if [1] = 'old' goto jump4 
if [1] = 'OLD' goto jump4 
igset mscf 1.5
key .05 1.5 20 '[r,c] ambiguities from sin(2[a]) sin(2[b])'
igset mscf 1.
dline -.1 .05 1.3  1.3 
text .1 1.3 '95"Y#CL zones from sin(2[a]) sin(2[b])' .3
jump4:
nobabar:

atit '[r]' '[c]'
return