Initial script to compare optimization and MCMC sampling