Have a working version of the new multivariate normal proposal distribution,
with accompanying tuning procedure. Have also split up pre_tune() as a
totally separate routine to sample() for MCMC objects.