t = seq( 0, +3, length.out=1000 ) f = rep( 1, length(t) ) f[ (0 <= t) & (t < 2) ] = t[ (0 <= t) & (t < 2) ]^2 #postscript("../../WriteUp/Graphics/Chapter6/chap_6_sect_3_prob_10_plot.eps", onefile=FALSE, horizontal=FALSE) plot( t, f, type='p', cex=0.25, pch=20, col='black' ) points( c(0, 2), c(0, 1), pch=20, col='black', cex=1.5 ) grid() #dev.off()