replaced scipy truncnorm() with kludgy version