%% randskripti.m
N=10000;
tulos=zeros(1,1000);
for k=1:N
    if rand =< 0.5
     tulos(k)=1;
    end 
end
P=cumsum(tulos)./(1:N);
plot(P)
grid on
shg
%%
N=1000;
T=rand(1,N);
kruuna=T<=0.5;
krlkm=cumsum(kruuna);
heittojenlkm=1:N;
Pkruunat=krlkm./heittojenlkm;
plot(Pkruunat)
grid on
shg
figure
plot(Pkruunat(end-100:end))
grid on
shg

