calc_mc_tk function

Conduct multiple TK simulations using Monte Carlo