powered by
#An matrices dicomposition function
gmedc(A, B = diag(nrow(A)))
a square matrix
a diagonal matrix