Program in Matlab or Python if possible

Program conditional Monte Carlo to sample the three-dimensional vector (X1, X2, X3) with probability mass function equal to: P(X1 = X1, X2 = X2, X3 = x3) 2x7x2(533)+x1x2 231=1 32=1 73=1 2y y2 (543)+y1y2 if 0 1