Added to defs.py the joint probability calculating function for Yt, Yhist, XXhist.