title: From for loops to numpy dot implementation
body: <pre><code>N=100 
  
  reliab=zeros((N,N))
  
  
  for i in range(N):    
  
      for j in range(N):
          if random() &lt; 0.6:
              reliab[i,j] = 1
          else:
              reliab[i,j] = 0
  </code></pre>
  
  <p>As in my code this matrix filling is recalled a huge number of times, these for loops should be changed with a dot product...but I dunno how to do it...
  Anyone willing to help me?</p>