Custom median / normal

file1.r

# set.seed(42)
X = 25
Y = 26
N = 100
set.sd = 5

set <- rnorm(N, Y, set.sd)
set.left <- set[set < X]
set <- c(set.left, X + (X - set.left))

set <- c(set, 
         X + ((set.sd / 2) * sign(X - Y)),
         ((length(set)+ 2) * Y) 
         - sum(set, X + ((set.sd / 2) * sign(X - Y)))
         )
median(set)
mean(set)
hist(set)
abline(v= median(set), col = 'green', lwd = 3)
abline(v = mean(set), col = 'red', lwd = 3)