Got compute_c_q0 working, now need to write a comparison with simulations