% Volterra kernel mask 

clear

N=25;
ma=mask4(N);
m4=zeros(1,nchoosek(N+4-1,4));

i=1;
for k1=1:N
    for k2=k1:N
        for k3=k2:N
            for k4=k3:N
                m4(i)=ma(k1,k2,k3,k4)/24;
                i=i+1;
            end
        end
    end
end

save('../mat/mask_P4_N25a','m4');

