title: binom.pmf only returning zero
body: <p>code:</p>
  <pre><code>def expected_profit(n):
      total = 0
      X = np.arange(0,n+1)
      p = np.arange(0,n+1)
      profit = np.arange(0,n+1)
      for i in list(range(1,n+1)):
          print(&quot;X_i:&quot;, X[i])
          p[i] = binom.pmf(X[i],n,19/20)
          print(p[i])
          if X[i] &gt; 100:
              profit[i] = 50*n-60*(X[i]-100)
          else:
              profit[i] = 50*n
          total += profit[i]*p[i]
      return total
  expected_profit(10)
  &gt;&gt;&gt;0
  </code></pre>
  <p>For some reason, after each iteration, <code>p[i]</code> is equal to zero. Yet when I manually type out (for example) <code>binom.pmf(10,10,19/20)</code> I get a non zero answer. What is the problem here?</p>
  <p>This seems to happen with any call to <code>binom.pmf</code> within the function call.</p>