function y=dproducto(x)
j=find(x==0);
l=length(j);
if l==0
    y=prod(x)./x;
elseif l==1
    y=zeros(length(x),1);
    x(j)=[];
    y(j)=prod(x);
else
    y=zeros(length(x),1);
end