tmpar <- par(mfcol=c(1,2))
on.exit( par( tmpar ))

plot(catalyst)
attach(catalyst)
interaction.plot(Conc, Cat, Yield)
detach("catalyst")