Calcul de la fonction Primorielle P#

rem ============================================================================
rem               La fonction Primorielle
rem ============================================================================
rem La fonction primorielle de p, notée p# (lire primorielle de p), est le produit
rem de tous les nombres premiers inférieurs ou égaux à n.
rem Par exemple, la primorielle de 10 est : 10# = 2*3*5*7 = 210  
rem ============================================================================
rem Le premier nombre premier étant 2, le produit p# commencera toujours par 2.
rem On en déduit que p# sera toujours un nombre pair, quelque soit le nombre p.
rem ============================================================================
rem Exemples :
 print "2# = "  + str$(Primorielle(2))  : ' 2
 print "5# = "  + str$(Primorielle(5))  : ' 2*3*5 = 30
 print "7# = "  + str$(Primorielle(7))  : ' 2*3*5*7 = 210
 print "8# = "  + str$(Primorielle(8))  : ' 2*3*5*7 = 210
 print "9# = "  + str$(Primorielle(9))  : ' 2*3*5*7 = 210
 print "10# = " + str$(Primorielle(10)) : ' 2*3*5*7 = 210
 print "11# = " + str$(Primorielle(11)) : ' 2*3*5*7*11 = 2310
 print "17# = " + str$(Primorielle(17)) : ' 2*3*5*7*11*13*17 =  510510
 print "23# = " + str$(Primorielle(23)) : ' 2*3*5*7*11*13*17*19*23 = 223092870
 print "29# = " + str$(Primorielle(29)) : ' 2*3*5*7*11*13*17*19*23*29 = 6469693230
 print "37# = " + str$(Primorielle(37)) : ' 2*3*5*7*11*13*17*19*23*29*31*37 = 7420738134810

end
rem ============================================================================
' test de primalité de n
' retourne 1 si n est premier, 0 sinon
FNC Est_Premier(n)
    dim_local a, sortir
    a = 3 : sortir = 0
    if n < 1
       message "Nombre entier strictement positif SVP !"
       exit_fnc
    end_if
    if n = 1 then result 0 : exit_fnc :' 1 n'est pas premier
    if n = 2
       result 1 : sortir = 1 : ' 2 est premier
    else
       if mod(n,2) = 0 then result 0 : sortir = 1 : ' n n'est pas premier
    end_if
    while sortir = 0
        if a >= sqr(n) + 1
           result 1  : sortir = 1 : ' n est premier
        else
           if mod(n,a) = 0 : ' si n est multiple de a
              result 0 : sortir = 1 : ' n n'est pas premier
           end_if
        end_if
        a = a + 2
    end_while
END_FNC
rem ============================================================================
' Calcul de la primorielle de p, notée p#
FNC Primorielle(n)
    dim_local i,p
    if n < 2
       message "Nombre entier strictement > 1 SVP !"
       exit_fnc
    end_if
    p = 1    
    for i = 2 to n
        if Est_Premier(i) = 1 then p = p * i
    next i
    result p    
END_FNC
rem ============================================================================