Matlab - Grafica Matlab

clear all

 figure(gcf);

 x=0:0.01:1;

 colores={'r','g','b','m','k'};

 alfa=[0.5 5 1 2 2];

 betha=[0.5 1 3 2 5];

 hold on

 for i=1:length(alfa)

 y=betapdf(x,alfa(i),betha(i));

h= plot(x,y,colores{i});

 end

 hold off

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

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

 title('Beta')

axis([0 1 0 2.6])

x = (0:0.1:1);

a=2;b=2;

y = zeros(size(x));

%r(n)=(n-1)!

%r(n+1)=nr(n)

for i = 1:size(x,2)

    % Aplicacion de la Funcion de Densidad de Probabilidad

    y(i)= (a*(factorial(a-1))/(factorial(a-1)*(factorial(b-1))))*(i.^(a-1)*(1-i).^(b-1));

end

subplot(212),plot(x,y), hold on

title('Densidad de Probabilidad')

clear all

figure(gcf)

x = 0:0.01:1;

syms X

a=2;b=2;

y = zeros(size(x));

%r(n)=(n-1)!

%r(n+1)=nr(n)

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)= (a*(factorial(a-1))/(factorial(a-1)*(factorial(b-1))))*(x(i)^(a-1)*(1-x(i))^(b-1));

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')