this code is WRONG!

First, it should be called quartile1st and quartile3rd, not quantile*

Second, it already RELIES on being (partitially) sorted since it just looks at the center element and THEN takes the median of the interval... arghh.

Also, in the MEDIAN function we could use nth_element. No need to sort it all.