%swerling
clear all;
clc;
N=1024;
u1=rand(1,N);
u2=rand(1,N);
for i=1:N;
    x(i)=sqrt(-2*log(u1(i)))*cos(2*pi*u2(i));
    y(i)=sqrt(-2*log(u1(i)))*sin(2*pi*u2(i));
end;
r1=5;
r=r1*(x.^2+y.^2);
subplot(211);
plot(r);
subplot(212);
pyulear(r,8,'twosided');