/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        E N C O D A G E   D E S   D E C I M A L E S   D E   ' P I '   G R A C E   A U X   N O M B R E S   P R E M I E R S          */
/*        " A   L A   G O D E L "  :                                                                                                 */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtb/SuitesGodel.01.SubmissionSloane.pi' :                                                                      */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 20220320113103).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
scale = 1000;

print     "\n";

pn[0]=2;pn[1]=3;pn[2]=5;pn[3]=7;pn[4]=11;pn[5]=13;pn[6]=17;pn[7]=19;pn[8]=23;pn[9]=29;pn[10]=31;pn[11]=37;pn[12]=41;
pn[13]=-1;
                                        /* Definition des premiers nombres premiers ("Prime Numbers").                               */

pi[0]=1;pi[1]=4;pi[2]=1;pi[3]=5;pi[4]=9;pi[5]=2;pi[6]=6;pi[7]=5;
pi[8]=-1;
                                        /* Definition des premieres decimales de 'PI'...                                             */

define    pow(x,n)  {
                    produit = 1;
                    i = 1;
                    while    (i <= n)
                             {
                             produit = produit*x;
                             i=i+1;
                                        /* L'utilisation de 'i++' provoquerait l'edition de 'i', d'ou 'i=i+1'...                     */
                             }
                    return(produit);
                    }

iterer = 0;
i1 = 0;
i2 = 0;
while     (iterer == 0)
          {
          if        ((pn[i1] > 0) && (pi[i1] >= 0))
                    {
                    term = 1;
                    for       (i2 = 0 ; i2 <= i1 ; i2 = i2+1)
                              {
                              term = term*pow(pn[i2],pi[i2]);
                              if        (i2 == i1) print     "a(",i2+1,")=",term,"\n";
                              }
                    } else
                    {
                    if        (pn[i1] <= 0) print     "OVERFLOW\n";
                    iterer = 1;
                    }
          i1 = i1+1;
                                        /* L'utilisation de 'i1++' provoquerait l'edition de 'i1', d'ou 'i1=i1+1'...                 */
          }
quit;



Copyright © Jean-François Colonna, 2022-2022.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2022-2022.