p <- seq(0,1,.01)
n <- 100

mse1 <- p*(1-p)/n

mse2 <- p*(1-p)/n + .0009

mse3 <- n*p*(1-p)/106^2 + ((n*p-106*p+3)/106)^2
 
plot(p, mse2, type="l", col="blue", lwd=3, ylim=c(0,.0035), ylab="MSE")
points(p, mse1, type="l", col="green", lwd=3)
points(p, mse3, type="l", col="purple", lwd=3)