read.tri <- function(file)

{x <- scan(file)
lx <- length(x)
d <- (sqrt(8*lx+1)-1)/2
m <- matrix(0, d, d)
m[upper.tri(m, T)] <- x
m <- m + t(m) - diag(diag(m))
return(m)}