First let's compute the analytical solution to the polynomial fitting problem
python analytic_solution.py
Then let's check that our estimator is indeed unbiased and that the analytic covariance is correct
python montecarlo.py
Let's code up our own mcmc and let's plot the histogram of the chains
python run_mcmc.py
Finally, let's use cobaya mcmc to solve our problem (and getdist to plot the chains)
python run_mcmc_cobaya.py
For part4 the idea is to adapt these functions for our new model and the new data set data_example_precise.txt