Instructor: Professor Barbara Bailey
Additional Dynamic Links:
symsqrtinv <- function(x) { ee <- eigen(x) ee$vectors%*%diag(sqrt(1/ee$values))%*%t(ee$vectors) }