function ret=mygaussian_hw1(xmean, xsigma)

  ret=[];
  while (isempty(______))
   x=________+(0.5-rand)*8*xsigma; %4 sigma range
   if rand<(__________)
    ret=x;
   end
  end
    
    
endfunction
  
