CMBACCUR

The program optimizes the microwave background radiation.

Features:as described in hubbletension4.pdf we optimize the likelihood of the CMB. The likelihood is calculated with Bayesian statistics. You can optimize 4 parameters :1. the amplitude c_1, 2. the baryon mass,omega_b*h^2, 3.total mass omega_m, 4. spectral density n_s. We do this for a constant angle of the CMB, teta and tau_reionisation=0.06.In effect we optimize ksi=-2*log(likelihood) where ksi is approached from far above 1 to a value near 1. The maximun likelihood=1.

Use the main optm1.

This main starts 'fminserach' for 3 parameters: c_1, the amplitude, the baryon mass, the spectral density. The 4th parameter omega_m is also calculated as it depends only on omega_bh^2 and theta which are then known . This is done by 'iter3' as an iteration of 'findteta = teta' . All is done by the routine 'om'. This routine then apllies 'cmbaccur' which returns ksi and the likelihood. ksi is optimized as a minimum. The optimization depends on start values op the parameters. With different start values of the parameters and h we should find the same output.

'cmbaccur' uses 7 parameters, all single numbers: 1.the amplitude c_1, 2. omega_m, 3. omb_h^2, 4.Hubble factor h, 5. spectral density n_s, 6. tau_reionisation, 7. the accuracy (of the likelihood) with 1 percent for the input value 1 and 0.1% accuray for the value 0. The Hubblefactor h used is 0.67 from the Planck people or 0.73 from the Iae-people.