# chlsoi.s                         J F Monahan, June 2001
 # Inverse of positive definite matrix -- uses 'solve'
A<- c(4, 2, 2, 4, 2, 5, 7, 0, 2, 7, 19, 11, 4, 0, 11, 25 )
 A <- matrix(A,4,4)
 "matrix"
 A
 Ainv <- solve(A)
 "inverse"
 Ainv
 rm(list=ls())