% Volterra kernel mask 

clear

N=100;
ma=mask3(N);
m3=zeros(1,nchoosek(N+3-1,3));

i=1;
for k1=1:N
  for k2=k1:N
    for k3=k2:N
      m3(i)=ma(k1,k2,k3)/6;
      i=i+1;
    end
  end
end

save('../mat/mask_P3_N100a', 'm3');

