function(progarray = progfile, jump = T, product = T, device = x11(geometry = 
	"640x920", close = T))
{
	device
	par(mfrow = c(2, 1))
	par(mar = c(7, 6, 4, 4) + 0.1)	#
	if(!product) {
		for(i in 1:length(progarray))
			simreact(progarray[i], jump = jump)	#
	}
	else {
		mmm.tau <<- numeric()
		for(i in 1:length(progarray))
			simprod(progarray[i], jump = jump)	#
		xx <<- mmm.tau[, 1]
		yy <<- 1/mmm.tau[, 2]
		if(jump)
			xlabel <- "[E] + [S]"
		else xlabel <- "[S]"
		plot(xx, yy, pch = 1, lab = c(5, 4, 6), las = 1, ylab = "",
			xlab = "")
		axes(axes = F, mgp = c(3, 1, 0), ylab = "", xlab = xlabel)
		axes(axes = F, mgp = c(5, 1, 0), ylab = "1/tau", xlab = "")
		temp <- nlsfit(twostep.fit)
		lines(spline(xx, yy - temp$residuals))
		xxx <- signif(temp$coef, 4)
		title(cex = 0.667, main = paste("k-2 =", xxx[1], "  k2 =",
			xxx[2], "  k-1/k1 =", xxx[3]))
	}
	invisible()
}