Add cahcematrix.R