1. Napisać funkcję generującą wartości losowe z rozkładu gamma τkxk−1 exp(−xτ)Γ(k), dla parametrów k = 9,τ = 0.5 (2 pkt.)
funkcja.m
function [ gamma ] = funkcja( x,t,k )
gamma=(t.^k).*(x.^(k-1)).*(exp(-x*t))
end

while liczba <= 1000
x=rand()*zakres;
y=rand();
if y < tau^k*x^(k-1)*exp(-x*tau)/gamma(k)
wyniki(liczba)=x;
liczba=liczba+1;
end
end





2. Narysować krzywą gęstości (1 pkt.)


krzywa.m
x=(1:60);
t=0.5;
k=9;
funkcja(x,t,k);
plot(x,ans);






3. Stworzyć histogram z wylosowanych 1000 próbek (2 pkt.)

* Stworzyć animację z losowań do histogramu (2 pkt.)

histogram.m
clear
for i=1:1000
x=rand(1,1000);
t=0.5;
k=9;
funkcja(x,t,k);
hist(ans);
drawnow
pause(0.5);
end

http://www.mathworks.com/help/stats/makedist.html