function [x,y]=boxmuller(u,v)
% the box muller transformation.
r=sqrt(-2*log(u));
theta=2*pi*v;
x=r.*cos(theta);
y=r.*sin(theta);

