# -*- coding: utf-8 -*-
"""
Created on Fri Dec 17 17:01:09 2021

@author: Sim
"""
import numpy as np
n=10
p=0.3
def bern(n,p):
  x = np.random.uniform(0,1,n)
  return list(map(int, (x < p)))

xx = []
for i in range(1000):
    xx = xx + [sum(bern(10,0.3))]
xx = np.array(xx)
print(xx.mean())
print(xx.var(ddof=1))