Using R studio software generate 500 random deviates for X from a B(10, 0.3) distribution and 500 random deviates for Y from a B(5, 0.3) distribution. Add corresponding random deviates from each distribution to form an empirical W=X Y. Then use the theoretical result of Example 3.8.2 and directly generate another 500 random deviates for W, call it W1 from a B(15, 0.3). Order the result of the sum of W=X Y and of W1. As an empirical test of the equivalence of the distributions, you could do a few things. You could create a scatterplot of the random deviates of W1 on the random deviates of W=X Y. It won't be real pretty, but values should follow a line y=x. Or you could plot histograms of W1 and W on the same graph. (In R you call the two histogram commands and on the 2nd add the argument, add=T)