from sympy import * from sympy import stats
n, p = symbols('n p') for n in [10,12,14,16,18]: X = stats.Binomial('X', n, Rational(1,2)) if stats.E(X**2)==stats.variance(X)+25: print(n)
10