R help - Nov 2010 - matrix limit

Dear experts,

I need to write R code for lim(I+X/n)^n with "repeat loop". (note:
limit is
n from 1 to infinity, I is identity matrix, X is square symmetric matrix)
Can anyone please provide help with my question.

Thank you. 
-- 
View this message in context:
http://r.789695.n4.nabble.com/matrix-limit-tp3035508p3035508.html
Sent from the R help mailing list archive at Nabble.com.

zhiji19 <zhiji19 <at> gmail.com> writes:
> I need to write R code for lim(I+X/n)^n with "repeat loop".
(note: limit is
> n from 1 to infinity, I is identity matrix, X is square symmetric matrix)
> Can anyone please provide help with my question.
> 
 Here's a more efficient solution. If you *need* to use a repeat { }
loop, then I am guessing this would be a homework problem ... see
?"repeat" (although that does not actually provide very much guidance
...)

library(expm)  ## needs to be installed first

z <- matrix(1:9,nrow=3)

tmpf <- function(M,n) {
  (diag(nrow(M))+M/n)%^%n
}
tmpf(z,1)
tmpf(z,20)
tmpf(z,100)