function [erro, erro2, maxerr] = mde(X,D)

  % erro: Mean Deviation Error (MDE)
  % erro2: Given Distance Mean Error (GDME) from Luo, Wu, JOTA 2011

  n = length(D);
  erro = 0;
  erro2 = 0;
  qtd = 0;
  maxerr=0;
  
  for i=1:n
    for j=(i+1):n
      if (D(i,j)~=0)  
        errij = abs(norm(X(:,i) - X(:,j), 2) - D(i,j));  
        if errij > maxerr
            maxerr = errij;
        end
        erro = erro + errij/D(i,j);
        erro2 = erro2 + errij^2;
        qtd = qtd + 1;
      end
    end
  end
  
  erro = erro/qtd;
  erro2 = sqrt(erro2/qtd);

end