gama fonksiyonu

 clc

clear

k=5;

L=10;

n=1000;

for j=1:n

for i=1:k

    u=rand;

    x(i)=-L*log(1-u);

end

y(j)=sum(x);

end

mean(y)

hist(y)