function res = integre4(v1,v2,v3,v4)

l = length(v1);
s = 0;
for k = 1:l
    for kk = 1:l
    s = s+v1(k)*v2(k)*v3(kk)*v4(kk);
    end
end

res = s/(l^2);

end
