/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   D E   L ' E N T R O P I E   D E   C O N F I G U R A T I O N  :                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtb/EntropieConfiguration.01' :                                                                                */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 20210821115356).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
define    factorielle(n)
          {
          if   (n<=1)         {
                              return(1);
                              }                   else                {
                                                                      return(n*factorielle(n-1));
                                                                      }
          }
                                        /* Definition de la fonction factorielle de 'n'.                                             */

define    combinaison(n,k)
          {
          return(factorielle(n)/(factorielle(k)*factorielle(n-k)));
          }
                                        /* Definition de la fonction combinaison de {n,k}.                                           */

scale=0;

combinaison(6,0);
combinaison(6,1);
combinaison(6,2);
combinaison(6,3);
combinaison(6,4);
combinaison(6,5);
combinaison(6,6);
print     "\n";
                                        /* Valeurs de test...                                                                        */

combinaison(100,50);
print     "\n";
                                        /* Valeur apparaissant dans "Jusqu'a la fin des temps" de Brian Greene (page 42)...          */

n=64;

for       (k=0 ; k<=n ; k++)  {
                              cnk=combinaison(n,k);
                              log=l(cnk);

                              print     "C(",n,",",k,")=",cnk;
                              scale=2;
                              print     "     log(C(",n,",",k,"))=",l(cnk),"\n";
                              scale=0;
                              }

quit;