> a = np.floor(10*np.random.random((2,2)))
> a
array([[ 8., 8.],
    [ 0., 0.]])
> b = np.floor(10*np.random.random((2,2)))
> b
array([[ 1., 8.],
    [ 0., 4.]])
> np.vstack((a,b))
array([[ 8., 8.],
    [ 0., 0.],
    [ 1., 8.],
    [ 0., 4.]])
> np.hstack((a,b))
array([[ 8., 8., 1., 8.],
    [ 0., 0., 0., 4.]])

> from numpy import newaxis
> np.column_stack((a,b))  # With 2D arrays
array([[ 8., 8., 1., 8.],
    [ 0., 0., 0., 4.]])
> a = np.array([4.,2.])
> b = np.array([2.,8.])
> a[:,newaxis] # This allows to have a 2D columns vector
array([[ 4.],
    [ 2.]])
> np.column_stack((a[:,newaxis],b[:,newaxis]))
array([[ 4., 2.],
    [ 2., 8.]])
> np.vstack((a[:,newaxis],b[:,newaxis])) # The behavior of vstack is different
array([[ 4.],
    [ 2.],
    [ 2.],
    [ 8.]])