min <- 0                                 
max <- 1
n <- 2000; #sample size
num_trials <- 10000; #num samples
means <- list();
medians <-list();
for(tn in 1:num_trials) {
	min <- 0                                 
	max <- 1
	x <- runif(n,min,max); #draw from uniform distribution
    means[[tn]] <- mean(x) # record sample mean for each sample
	medians[[tn]] <- median(x) # record sample median for each sample
}

dev.new()
hist(unlist(medians),xlim=c(0.40,0.60),breaks=50,main="Distribution of Sample Medians")
dev.copy(jpeg,'hist_medians.jpg')

dev.new()
hist(unlist(means),xlim=c(0.40,0.60),breaks=50,main="Distribution of Sample Means")
dev.copy(jpeg,'hist_means.jpg')

print("compare mean of sample means and medians")
mean(unlist(means))
mean(unlist(medians))

print("compare sd of sample means and medians")
sd(unlist(means))
sd(unlist(medians))