1) Generate a sample using Bernoulli distribution

2) Define a function to calculate the negative maximum likelihood estimate for parameter p.

3) Apply either the nlm or mle function to find the estimate for the parameter of the distribution.

4) Repeat steps 1-4 for a sample of different size and for different initial parameter of p (you can create new chuck for this question)

Compare the results to the MLE for p derived analytically.