Let X1, X2, . . . , X20 be independent and identically distributed random variables with pmf P(X = x) = (1 p)^x-1*p , x = 1, 2, . .

By simulating 50000 sets of twenty outcomes when p = 1/3, compute the corresponding maximum likelihood estimates and plot a histogram of the outcomes. Need to do this as a matlab code.