substitution for all indexes



On Thu, May 29, 2008 at 10:11 AM, walrand <stephan.walrand at uclouvain.be> wrote:

This can probably be handled by Maxima, but I'm not sure I understand
your notation:

> a  follows a recursive law:
>
> a[n,k_n+1,k_n] = a[n-1,k_n+1,k_n] - a[n-1,k_n+1,k_n-1] *
> a[n-1,k_n,k_n-1] / a[n-1,k_n-1,k_n-1]
>
> I would like to apply this law for a given n, but for all the index k_i,
> in one single command

I think you want k_n to be related to n?  If so, it needs to be
written as k[n].  Does k_n+1 mean k[n+1] or k[n]+1?

Is the result you want a set of simultaneous equations with the above
relationship for n in 1..k?  Or a recursion -- in which case you need
to define the base case a[0,...].  What is the value of a[n,k[i],k[j]]
where i <> j+1?

I'm just not quite clear what you want to do here.

               -s