The fastest way to create a random network?

I am trying to create a square matrix of 1's and zeros, where 1's represents a communication link and 0's represents a lack of the same. In this matrix, the probability of observing 1 at each location is given by beta.

What is the most efficient way to create the stated network? I have tried the following, but there are no doubt better ways:

n = 100 beta = 0.25

A = rand(n,n) B = find(A<=beta) C = find(A>beta) A(B) = 1 A(C) = 0

Answer 1

[Roger Stafford]

   A = 1*(beta>=rand(n));