[nr nc] = BestArrayDims(n)
Given n elements, BestArrayDims returns the dimensions of an array in which to store them. The function also accepts user preferences as to whether they require less wasted space or a more square configuration.

I have found this function particularly well-suited for finding subplot dimensions if it is not known in advance the number of plots that will be needed.