`xch` <-
function(p,rho)
{
if (p <= 0) {stop('n must be at least 1')}
else if (p==1) {r = 1}
else if (p > 1)
{
c1 = rep(rho,p-1)
r = toeplitz(c(1,c1))
}
return(r)
}
`xch` <-
function(p,rho)
{
if (p <= 0) {stop('n must be at least 1')}
else if (p==1) {r = 1}
else if (p > 1)
{
c1 = rep(rho,p-1)
r = toeplitz(c(1,c1))
}
return(r)
}