title: Using scipy.stats.bernoulli to simulate 1/3 probability
body: <p>Our task is to recreate the famous fair coin flip (p=0.5) using only <code>bernoulli.rvs()</code> to generate function with uniform discrete distribution (outcomes = 0,1,2). Is it possible to do using only bernoulli function form <code>scipy.stats</code>?</p>