lambda = - 2(1+(2q - 1)βe^(-lambda tau)) for tau=0