function [theta1, theta2]=gaua(k,amu,h);
   global a1;
   v1=a1*randn(k,1);
   theta1=amu+v1';
   theta2=amu-v1';   
end;