meanboot<-function(object,R){
keep<-rep(0,R)
for(i in 1:R){
newsamp<-sample(object,replace=T)
keep[i]<-mean(newsamp)
}
return(keep)
}

age<-rexp(10,1)
junk<-meanboot(age,1000)
hist(junk)

n<-10
bootmean<-rep(0,100)
bootlb<-rep(0,100)
bootub<-rep(0,100)
normlb<-rep(0,100)
normub<-rep(0,100)
for(j in 1:100){
age<-rexp(n,1)
junk<-meanboot(age,1000)
bootmean[j]<-mean(junk)
bootlb[j]<-sort(junk)[25]
bootub[j]<-sort(junk)[975]
normlb[j]<-mean(age)-qt(.975,(n-1))*sqrt(var(age)/n)
normub[j]<-mean(age)+qt(.975,(n-1))*sqrt(var(age)/n)
print(j)
}

(100-(sum(bootlb>1)+sum(bootub<1)))/100
(100-(sum(normlb>1)+sum(normub<1)))/100