Problem passing arrays to linsolve( )

I changed your binom to binomial. Maybe you had defined a function binom
somewhere else?

Barton

-----maxima-bounces at math.utexas.edu wrote: -----


>Try?using?lists,?not?arrays:
>
>solution(p):=?block([eq?:?[],?v?:?[]],
>????for?i:0?thru?p?do?(
>???????eq?:?cons?(sum(binomial(j+1,i)*a[j],j,i,p)?=?binomial(p,i),?eq),
>???????v?:?cons(a[i],?v)),
>????linsolve(eq,v));
>