The joint probability distribution function method has been developed in P1; for reflections with rational indices. The positional atomic parameters are considered to be the primitive random variables, uniformly distributed in the interval (0, 1), while the reflection indices are kept fixed. Owing to the rationality of the indices, distributions like P(F(p1), F(p2)) are found to be useful for phasing purposes, where p1 and p2 are any pair of vectorial indices. A variety of conditional distributions like P(|F(p1)| | |F(p2)|), P(|F(p1)| |F(p2)), P(varphi(p1)| |F(p1)|, F(p2)) are derived, which are able to estimate the modulus and phase of F(p1) given the modulus and/or phase of F(p2). The method has been generalized to handle the joint probability distribution of any set of structure factors, i.e. the distributions P(F(1), F(2),ellipsis, F(n+1)), P(|F(1)| |F(2),ellipsis, F(n+1)) and P(varphi(1)| |F|(1), F(2),ellipsis, F(n+1)) have been obtained. Some practical tests prove the efficiency of the method.