Added a pair of functions for simulating critical values for Chi-squared
with df degrees of freedom