
function[x]=aanormalpdf(m,s,N)
%function[x]=aanormalpdf(m,s,N);
%s=standard deviation;m=mean value;
for i=1:N
    r1=rand;
    r2=rand;
    z(i)=sqrt(-2*log(r1))*cos(2*pi*r2);
end;
x=s*z+m;
