close to working analytic_tcprob.out