equations checked again
perturbations fixed. C_l, P_k and background agree well with CMBeasy