% Volterra kernel mask 

clear

N=100;
ma=mask5(N);
m5=zeros(1,nchoosek(N+5-1,5));

i=1;
for k1=1:N
    for k2=k1:N
        for k3=k2:N
            for k4=k3:N
                for k5=k4:N
                    m5(i)=ma(k1,k2,k3,k4,k5)/120;
                    i=i+1;
                end
            end
        end
    end
end

save('../mat/mask_P5_N100a m5');

