fix in formula for sampleTGivenZT()