Fonctions

Fonction Beta.

%--------------------------------------------------------
% utilisation fonction beta
% (c) JPaul Molina
%-------------------------------------------------------
%       fonction gamma
%             infini  x-1  -t
%  Gamma(x)= S       t    e  dt
%             0
%       fonction beta
%
%              gamma(m)gamma(n)
%   beta(n,m)= ----------------
%                 gamma(m+n) 
%
%              1  m-1     n-1
%  beta(n,m)= S  t   (1-t)   dt
%              0
%---------------------------------------------------------
echo on;
%-------------------------------------------------------
%           fonction beta
% évolution en fonction de n et de m
% cf fichier pour explications
%------------------------------------------------------
echo off;
clear all;
N=5; 		% maxi de n
M=6; 		% maxi de m
n=1:0.1:N;  
m=1:0.1:M;
for i=1:M
   plot(n,beta(n,i));
   hold on;
   str=sprintf('m=%d',i);
   y=beta(n,i);
   text(n(1,2), y(1,2),str);
end
	%y1=gamma(n).*gamma(m)/gamma(n+m);
xlabel('n');
ylabel('beta(n,m)');
title('fonction beta');
grid;