ps.options(horizontal=FALSE,onefile=FALSE)
ps.options(pagecentre=TRUE)
ps.options(paper="letter")
ps.options(width=7.0,height=5.0)

for (name in c("smooth", "filter") ) {

  particle <- read.table(paste(name,".csv",sep=""),header=T,sep=",")
  
  mean <- particle[,1]
  ulim <- mean + 2.0*particle[,2]
  llim <- mean - 2.0*particle[,2]
  true <- particle[,3]
  data <- particle[,4]
  
  ybot <- min(c(mean,ulim,llim,true))
  ytop <- max(c(mean,ulim,llim,true))
  
  xbot <- 1
  xtop <- length(mean)
  
  time <- seq(xbot,xtop)
  
  postscript(file=paste(name,"1.eps",sep=""))
  
  ybot <- min(c(mean,ulim,llim,true))
  ytop <- max(c(mean,ulim,llim,true))
  
  plot(c(xbot,xtop),c(ybot,ytop),type="n",ylab="",xlab="")
    lines(time,mean,lty="solid",col="red")
    lines(time,ulim,lty="dashed")
    lines(time,llim,lty="dashed")
  
  dev.off()
  
  postscript(file=paste(name,"2.eps",sep=""))
  
  ybot <- min(c(mean,true))
  ytop <- max(c(mean,true))
  
  plot(c(xbot,xtop),c(ybot,ytop),type="n",ylab="",xlab="")
    lines(time,mean,lty="solid",col="red")
    lines(time,true,lty="solid",col="blue")
  
  dev.off()
}