//////////////////// LoiSing1.sce ////////////////////////
//
// Construction de la fdr F_n de U_n=\sum_{i=1}^n 2^{-i}X_i
// les X_i etant iid de loi Bern(p)
n=input('Rentrer le nombre n : ');
p=input('Rentrer le parametre p de la loi de Bernoulli : ');
q=1-p;
M=[q p];
for i=2:n,
    M=[M.*q M.*p];
end
F=[cumsum(M) 1];
x=0:2^{-n}:1;
xbasc()
plot2d2("onn",x',F',[2])