par 関数で mfrow のパラメータを与えてあげればよい。
par(mfrow=c(2,3))
この場合、縦に 2 つ、横に 3 つのグラフが描画される。
例えば、こんな感じ。
par(mfrow=c(2,2))
plot(sin,-pi,2*pi)
plot(-10:10,-10:10,col="blue")
plot(-10:10,(-10:10)^2,type="l",col="red")
plot(-10:10,(-10:10)^3,type="l",col="green",lwd=4)