// EXO 4 // 1 function [y]=suivant(x,p) rand("uniform"); // On initialise la loi uniforme sur [0,1] comme rand(). if rand()
0 then t=t+1, s=suivant(s,p); end endfunction // 4 function [r]=ET(p) r=0; for i=1:1000 r=r+T(p); end r=r/1000; endfunction // La réponse semble être toujours autour de 10 // 5 // Conditionnellement à T>1, T-1 est la somme de 9 v.a. géométriques de // paramètre 1-p, donc E[T] = P(T=1) + P(T>1) * [1 + 9 *(1-p)^{-1}] // = p + (1-p) + 9 = 10