clear, clc;
figure, hold on;
NLeaves = 10^4; 
circ = rand(NLeaves,1)*2*pi; 
z = rand(NLeaves,1);
rad = .4*(1-z);
x = rad.*cos(circ); 
y = rad.*sin(circ); 
clr = rand(NLeaves,1);
scatter3(x.*clr,y.*clr,z,[],[zeros(NLeaves,1) (clr.*rad).^.6 zeros(NLeaves,1)],'*')
plot3(0,0,1.05,'rp','markers',20,'markerf','y')
[X,Y,Z] = cylinder(.03,50);
surf(X,Y,-Z*.1)
view(3, 9), axis equal off
title('FROHE WEIHNACHTEN');
