clear
Npoints = 100
% Test flyers.m
[foo, Nflyers, flyers] = pdf_flyers(Npoints, 2, 0.1, 5);
Nflyers
plot(foo(:,1), foo(:,2), 'bo',... 
     foo(flyers,1), foo(flyers,2), 'r+');
axis('equal')