srs.mean=function(data,x,n,N){
  data=as.name(data); x=x; n=n; N=N; fpc=(1-n/N)
  muhat=mean(x); s2=var(x)
  vhat.muhat=(s2/n)*fpc
  bound=round(2*sqrt(vhat.muhat),2)
  lower=muhat-bound; upper=muhat+bound
  cat("","\n","Results from SRS: Data =",data,"\n","N =",N,"n =",n,"\n","FPC =",fpc,"\n","Muhat =",muhat,"\n","Vhat muhat =",
      vhat.muhat,"\n","Bound =",bound,"\n","Lower Bound =",lower,"Upper Bound =",upper,"\n","")
  results=list(muhat=muhat,data=data,n=n,N=N,fpc=fpc,vhat.muhat=vhat.muhat,bound=bound,lower=lower,upper=upper,s2=s2)
}

# to use the function with its call:
# srs.mean(data,x,n,N)

# data: name of dataset, in quotes
# x: vector (data)
# n: sample size
# N: population size