Hello Sebastian, I think you may need to use tensorA package. You can define different indices and use Einstein or Reimann summation convention.
Grüße, Mehmet On 23 August 2013 15:14, Sebastian Hersberger <sebastian.hersber...@unibas.ch> wrote: > Hi all > > I have a short question relating to the usage of the summation sign in R. > > Let's define A and B as two kxk matrices (which are constant). > C^(n) = Ó_(j=1)^n [(Ó_(i=1)^(j-1) A^i ) B (Ó_(i=1)^(j-1) A^i)’ ] > > My goal is to calculate the matrix C for the periods from 1 to 200 (n=1-200). > > How has that to be implemented in R? > > Thank you for your help!! > > Kind regards, > Sebastian > > [[alternative HTML version deleted]] > > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
