function res = integre2(v1,v2,v3)

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

res = s;

end

