while l < 150:
    while X < 150:
        if X >= l:
            P1 = (math.factorial(X)) / (math.factorial(l) * math.factorial(X-l))
            P2 = (u**l) * ((1-u) ** (X - l))
            P = (P1 * P2) * number2_diag[X]
            nlist.append(P)
            X = X + 1

        else:
            P = 0 * X
            nlist.append(P)
            X = X + 1