function H=hproducto(x)
H=[];
n=length(x);
n1=n-1;
for j=1:n
    y=x;
    y(j)=[];
    y=dproducto(y);
    H=[H [y(1:j-1);0;y(j:n1)]];
end