from scipy.special import binom

def grade_of_service(n, p, c):
    prob = 0
    k = c + 1
    while k <= n:
        prob += binom(n, k)*(p**k)*((1-p)**(n-k))
        k += 1
    return prob