Anyone know how write a function that solves: (1 + c)x1 +x2 +x3 = 5 x1 +(1 + c)x2 +x3 = 5 + 2c x1 +x2 +(1 + c)x3 = 5 + 3c, where c is a small constant, for 1000 equidistant values c = (10^-14, 2*10^-14, ..., 10^-11) by using cholesky decomposition? /P -- View this message in context: http://r.789695.n4.nabble.com/matrix-help-tp2714378p2714378.html Sent from the R help mailing list archive at Nabble.com.