# Section 10; Question 8:
#
y = seq( 0, 1, length.out=100 )
fY2 = 2 * y
fY1p = 10 * y * (1-y^2)^4
fY5p = 10 * y^9
y_min = min( c( fY2, fY1p, fY5p ) )
y_max = max( c( fY2, fY1p, fY5p ) )

#postscript("../../WriteUp/Graphics/Chapter3/chap_3_sect_10_ex_8.eps", onefile=FALSE, horizontal=FALSE)
plot( y, fY2, type='l', ylim=c(y_min, y_max), col='black', xlab='y', ylab='pdfs' )
lines( y, fY1p, col='green' )
lines( y, fY5p, col='red' )
legend('topleft', c('f_Y2', 'f_Y_min', 'f_Y_max'), col=c('black', 'green', 'red'), lty=c(1, 1, 1))
grid()
#dev.off()