Matlab

Write a program to produce the probability mass function of the Bernouilli distribution, i.e. a function that only returns 2 possible values(0 and 1). The probability of returning 1 is p, the probability returning 0 is q. where p + q = 1. Using p = 0.4 and 100 samples verify that the expected value of such a distribution is p and the variance is p