Transient M/M/1 queue variance computation using generalized Q functions