Problem passing arrays to linsolve( )



I think you can do this with
f(p):=solve( makelist(sum(binomial(j+1,i)*a[j],j,i,p)=binomial(p,i),i,0,p),
          makelist(a[i],i,0,p))$