# parse cmd line args
args = commandArgs(trailingOnly=TRUE);
src = args[1];
dst = args[2];

# read data
data = read.table(src, colClasses="double")

# setup output
format = substr(dst, nchar(dst)-3, nchar(dst));
if (format == ".pdf" || format == "pdf~") {
   pdf(file = dst)
}

# plot data
require(gplots);
hist2d(data, nbins=200, col=c("white", heat.colors(12)));
#fac=0.05;
#rfac=1/fac;
#par(pty="s", pch=22, cex=fac, cex.axis=rfac, cex.lab=rfac, cex.main=rfac, cex.sub=rfac);
#plot.default(data, xlab="time [s]", ylab="neuron #");