t = seq( 0, +6, length.out=1000 )
f = rep( 1, length(t) )

f[(2<=t)] = exp(-(t[2<=t]-2))

#postscript("../../WriteUp/Graphics/Chapter6/chap_6_sect_3_prob_9_plot.eps", onefile=FALSE, horizontal=FALSE)
plot( t, f, type='p', cex=0.25, pch=20, col='black' )
points( c(0, 2), c(1, 1), pch=20, col='black', cex=1.5 )
grid()
#dev.off()