I have several random variables X = [x1, x2,...... xn] represented by columns in a Data matrix and rows are representing random samples. I also have the marginal Probability Density Functions as f(x1), f(x2), ... f(xn) for individual random variables. I would like to calculate their joint PDF as f(x1,x2,....xn).