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