cards <- 1:13
S <- cards
S[1] <- 11
S[11:13] <- 10
P <- function(s) sum(S==s)/13
p <- sapply(1:11, P)
sum((1:11)*p)