proc findq(ak,d,t,g,wg1,w1); 
   ak2=ak^2; 
   ql=(t+ak2*(wg1-t))./(d+ak2*(w1-d));
   m=rows(g);
   gg=zeros(m,1);
   gg[1:(m-1)]=g[2:m];
   gg[m]=g[m]-1;
   amu=wg1/w1; asd=sqrt(wg1/w1-amu^2);
   qq=selif(ql,ql.>=g .and ql .<=gg);
   retp(qq[1]);
endp;