2D uniform distribution sample

I am looking for a way to take a random sample of a uniform distribution in 2D square. Is there a single function for this? Or do I have to take the x and y axis separately uniformly distributed each like:

x=numpy.random.uniform(0,1,1000)
y=numpy.random.uniform(0,1,1000)
sample=numpy.column_stack((x, y))