MAT #CONVOLUTION(C,A,k)
computes k first coefficients of the convolution
of the columns of matrix A and saves them as
a column vector C.
It is assumed that elements a0,a1,a2,... of each A column are
coefficients of a polynomial a0+a1*x+a2*x^2+...
Default for k is k=(m-1)*n+1 where m,n are dimensions of A.
Alternatively:
MAT #CONVOLUTION(C,A,k,N)
when A has only one column computes the N-fold convolution of this column.
MAT #CONVOLUTION(C,A,B)
computes the convolution C of vectors (1st columns) of A and B.
..................................................................
Example: Probabilities of binomial distribution
as convolution of n Bernoulli distributions
p=1/3 n=13
MATRIX P ///
1-p
p
MAT SAVE P / Save probabilities of Bernoulli distribution
MAT #CONVOLUTION(C,P,n+1,n)
/MATSHOW C / See binomial probabilities
..................................................................
A = More about additional MAT #operations
M = More about MAT operations