#
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
# 
#-----

# temperature measurements: 
Z = c( 200, 202, 208, 204, 204, 207, 207, 204, 202, 199, 201, 198, 200,
       202, 203, 205, 207, 211, 204, 206, 203, 203, 201, 198, 200, 206,
       207, 206, 200, 203, 203, 200, 200, 195, 202, 204 )
Z = as.double(Z)
N = length(Z)

postscript("../../WriteUp/Graphics/Chapter2/prob_1_plot.eps", onefile=FALSE, horizontal=FALSE)
plot( Z, pch=c("o"), col=c("black"), xlab="Time (minutes)", ylab="Temperature" )
dev.off()


postscript("../../WriteUp/Graphics/Chapter2/prob_1_plot_ztp1_vs_zt.eps", onefile=FALSE, horizontal=FALSE)
plot( Z[1:(N-1)], Z[2:N], xlab="z_t", ylab="z_{t+1}" )
dev.off()

postscript("../../WriteUp/Graphics/Chapter2/prob_1_plot_ztp2_vs_zt.eps", onefile=FALSE, horizontal=FALSE)
plot( Z[1:(N-2)], Z[3:N], xlab="z_t", ylab="z_{t+1}" )
dev.off()