function d=digse(ra,r)
  // ra e r podem ser escalares, vetores ou matrizes //
  [m,n]=size(r); d=zeros(m,n);
  if ((size(ra,1)~=m) | (size(ra,2)~=n))
     error('dimensoes devem concordar')
  end
  for i=1:m
    for j=1:n
      rij=r(i,j);raij=ra(i,j);
      if (rij==raij) | (abs(rij)<%eps)
        d(i,j)=16.0;
       // usar 32.0 em maquina de 64 bits //
      else
        d(i,j)=max(0,-log10( 2*abs( (raij-rij)/rij  )));
      end
    end
  end
endfunction  // digse //
