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

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

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