import numpy as np
from scipy import stats as st

n, p = 10, 0.5

b = st.binom.pmf(8, n=n, p=p) + st.binom.pmf(9, n=n, p=p) + st.binom.pmf(10, n=n, p=p)

print('La probabilidad de más de 7 caras es:', b)