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

hist(unlist(x))

dev.new()
#hist(unlist(medians),xlim=c(0.40,0.60),breaks=50,main="Distribution of Sample Medians")
hist(unlist(medians),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")
hist(unlist(means),breaks=50,main="Distribution of Sample Means")
dev.copy(jpeg,'hist_means.jpg')

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

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