joint probability function