Bug fix:  qstsb --> sqrt(qstsb).