#
# 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.
# 
#-----

N=20

rhos = c(1,2/3)
for( ii in 2:(N+1) ){
  newRho = rhos[ii-1] - 0.5*rhos[ii-2]
  rhos = c(rhos,newRho)
}

rhos 

#postscript("../../WriteUp/Graphics/Chapter3/prob_6_plot.eps", onefile=FALSE, horizontal=FALSE)
plot( 0:N, rhos, pch=c("o"), col=c("black"), xlab="index (k)", ylab="rho_k" )
#dev.off()