This function is to decide the r setting given specific density in each dimension
proportion(m, n_dim)
integer, the number of grid points in each dimension, and we suggest m around 20 for 3 dimension
integer for the dimension, which is equal to the number of sub-population plus 1
matrix of setting the proportion of the population by given specific dimension and density in each dimension