n = 100
Q = Matrix(toeplitz(c(1, -0.5, rep(0, n-2))))
map <- contourmap(mu = seq(-5, 5, length=n),Q,n.levels = 2)
cols = contourmap.colors(map, col=heat.colors(100, 1),
credible.col = grey(0.5, 1))Run the code above in your browser using DataLab