let alpha = 0.05
let m = 3
let v = 38
let zeta = .5*(1 - exp(ln(1-alpha)/m))
let TSTAR = tppf(zeta, v)
.end of calculations