proc gaua(k,amu,h);
   v1=a1*rndn(k,1); 
   theta1=amu+v1';
   theta2=amu-v1';
   retp(" ");
endp;