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.hersberger at 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 at 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.
>