Write a script to simulate a sampling procedure from 1,000 identical lots using the binominal distribution with probability of a part being defective of 0.01. Create a histogram showing the frequencies of the number of observed defective parts in 1,000 lots using 10,000 trials. On the same set of axes plot the scaled probability density function for the binominal distribution with the appropriate parameters. The density function can be plotted over the range from 0 to 25. Compare the empirical and the theoretical distributions. Note that the binomial distribution function binopdf returns non-zero values only when the argument is an integer.