// Question 7 : Méthode de rejet
///////////////////////////////////////////////////////////////////
// Remplacer les ... par la commande appropriée
///////////////////////////////////////////////////////////////////

function [x]=rejetbeta(N)
x=zeros(1,N);
for i=1:N,
  // Mettre une réalisation d'une variable aléatoire
  // suivant une loi beta de paramètre (2,2) dans x(i)
  ...
end;
endfunction

N=10000; // Nombre de réalisations de la variable aléatoire générée

z=rejetbeta(N);

// Tableau de discrétisation des abscisses
// L'intervalle (0,1) est divisé en 24 intervalles de meme longueur
x=linspace(0,1,25); 

xbasc();

histplot(x,z); // Trace l'histogramme de z échantillonné suivant x

// Calcule y, l'image de x par la fonction densité
// de la loi beta de paramètre (2,2)
y=6*x.*(1-x).*(x<1 & x>0);

// Trace y en fonction de x
plot2d(x,y);