joint distribution function