do it in Matlab programming

(do not do it in shourt answer ) this is a coding problem :

Suppose that there are five drones parked at five marked positions. All drones take off at the same time, fly around and park at the marked positions randomly (i.e., all parking arrangements are equally likely). What is the expected number of drones parking at their original position? How about 10 drones with 10 positions? Conduct both theoretical analysis and simulation to find out the expectation. Hint: You are asked to find out the expectation only, not the distribution. Let Xi be a Bernoulli random variable. Xi=1 represents drone i parks at its original position and Xi=0 otherwise. Drone i takes one of the five positions with equal probability. Xi and Xj are not independent for i#j. Let X be the number of drones parking at their original position. Then X=X1+...+X5