beta=0|b_tru|(2*b_tru)|(3*b_tru);
n_al=31;

al=seqa(-3,0.2,n_al); 
p_al=zeros(n_al,1);
p_al[1]=cdfn((al[1]+al[2])/2);
i_al=2;
do while i_al.lt n_al;
p_al[i_al]=cdfn((al[i_al]+al[i_al+1])/2)-cdfn((al[i_al-1]+al[i_al])/2);
i_al=i_al+1;
endo;
p_al[n_al]=1-sumc(p_al);

al~p_al;
sumc(al'*p_al);; sumc((al.^2)'*p_al)-sumc(al'*p_al)^2;

p_0=.5*ones(n_al,1); 


{p1,p0}=probs(ga,beta,al);
pdat=p1*(p_0.*p_al)+p0*((1-p_0).*p_al);


botim=bosec;

npoi=201;
dpoi=2/(npoi-1);
sss1=seqa(-1,dpoi,npoi);
npoi=201;
dpoi=2/(npoi-1);
sss2=seqa(-1,dpoi,npoi);
dbe=1/10;
dga=1/5;


res=1000+zeros(1,3);
res[1,1]=-1;
"start";
print /flush " ";

kk=1;

for igam (1,rows(sss1),1);
bt=bosec;

print /flush igam;
g=ga+sss1[igam]*dga;


for ibet2 (1,rows(sss2),1);


bbb=b_tru+sss2[ibet2]*dbe;
be=0|bbb|(2*bbb)|(3*bbb); 
{p1,p0}=probs(g,be,al);
aa=p1~p0;
a=aa|ones(1,cols(aa));

pdata=pdat|1;
/* these are the constrains */

xm=rows(a); xn=cols(a);
aa=(0~ones(1,cols(a)) )| (pdata~-a)|zeros(1,(cols(a)+1));
aa=aa';
xmp=cols(aa); xnp=rows(aa); xm1=0; xm2=0; xm3=rows(pdata);
xica=0; xizro=zeros(1,xn); xiposv=zeros(xm,1);

eps=1.0d-12;
dllcall isimpl(aa,xm,xn,xmp,xnp,xm1,xm2,xm3,xica,eps); 

if (xica.ne -1); res=res|(xica~g~bbb);  endif;

kk=kk+1;


endfor;
print /flush bosec-bt;
endfor;



bt=bosec;

"HEREeeewf timqwe";  

bt-botim;