fixed bug where wrong mf was set to NAN for t=0

also, currently there is a bug where the range of the
digamma solution is negative, so the mstep must be
calculated incorrectly.