S <- mat.or.vec(101,2);
for (i in 1:101) {
  x <- i-50.0;
  S[i,1] <- x;
  S[i,2] <- exp(x)/(1.0 + exp(x));
# S[i,2] <- (exp(x)-1.0)/(2.0 + 2.0*exp(x));
}

source("psopts.r");
ps.options(width=5.0,height=5.0);
postscript(file="squash.eps");

par(mar=c(2,2,0,0));

plot(S,xlab="",ylab="",type="l");

dev.off();