How to repeat generating random sample 1000 times without for loop and draw the plots of alpha estimators?

Question

I have to generate one random sample with n= 20, from beta distribution(alpha=5,beta=1) and calculate MLE for alpha. I did this using the function rbeta and ebeta: