n=length(A)-1; m=length(x); B=zeros(m,n-1); wun=ones(m,1); B(:,1)=A(1)*wun; for j=2:n B(:,j)=A(j)*wun+x.'.*B(:,j-1); end y=A(n+1)*wun+x.'.*B(:,n); y=y.' B