Matlab - Grafica Estimada

function [b,tp] = beta34v2(n)

% Escribe un codigo que genere numeros aleatorios con distribucion 

% B(2, 2) ~ (f(x) = 6x(1 - x), 0 >= x >= 1) mediante el metodo de 

% aceptacion-rechazo simple

% Input:

% n :   número de muestras

% Output:

% b :   muestras aleatorias generadas y aceptadas

% tp :  total de puntos generados, aceptados y rechazados 

a = 0.5;

b = []; tp = [];

for i = 1:n,                            % Hasta que tenga todas las muestras

    u1 = rand;                          % Generar u1

    u2 = a * rand;                      % Generar u2

    while (u2 > 6 * u1 * (1-u1) )  % Mientras (U * a * g(x))> f(x) 

        % Bloque para generar x

        tp = [tp ; u1];                 % Recoger todas las muestras

        u1 = rand;                      % Volver a generar u1

        u2 = a * rand;                  % Volver a generar u2

    end

    x = u1;                             % x generado

    b = [b ; x];                        % Recoge las muestras aceptadas

    tp = [tp ; u1];                     % Recoge las muestras aceptadas

end

clear all

figure(gcf)

x = 0:0.01:1;

syms X

a=2;b=2;

y = zeros(size(x));

N=int(X^(a+b)*exp(-X), 0, inf);

D1=int(X^(a)*exp(-X), 0, inf);

D2=int(X^(b)*exp(-X), 0, inf);

for i = 1:size(x,2)

% Aplicacion de la Funcion de Densidad de Probabilidad

y(i)=N/(D1*D2)*(x(i)^(a-1)*(1-x(i))^(b-1));

end

 plot(x,y)

 set(gca,'xtick',0:0.1:1),

 set(gca,'ytick',0:0.2:2.6),

title('Densidad de Probabilidad Teorica')