Rscript -e 'x <- rnorm(1000, mean=50, sd=10); print(summary(x))'

rnorm() draws random numbers from a normal distribution, here 1000 samples with mean 50 and standard deviation 10. summary() then reports its min, quartiles, mean, and max. The sample mean and median should land near 50 for a large draw. Set a seed with set.seed(42) first to make the random draw reproducible.