calculErlangQNAP.py

# pour acceder aux fonctions maths evoluees
import math

# pour RHO/Erlang de 1 a 80
for p in range(1,80+1):
    # pour C/capacite de RHO a RHO*2
    for c in range(p,p*2+1):
        # calcul numerateur
        numerateur = pow(p,c) * 1.0 / math.factorial(c)
        # init denominateur a 0
        denominateur = 0
        # Somme dans le denominateur
        for i in range(0,c):
            denominateur += pow(p,i) * 1.0 / math.factorial(i)
        # affichage resultat (redirection vers un fichier donne un fichier utilisable par gnuplot)
        print( str(c) + " " + str(p) + " " + str(numerateur / denominateur))