ConsensusFnChange <- function(old,new,pos,clusters,amat,alpha) {
    n = length(clusters)
    result = 0

    for (j in 1:n) 
     {
        if (j!=pos)
        {
            if (clusters[j]==old)
            {
                if (j<pos)
                    result = result - (amat[j,pos] - alpha) 
                else
                    result = result - (amat[pos,j] - alpha) 
            }
            if (clusters[j]==new)
            {
                if (j<pos)
                    result = result + (amat[j,pos] - alpha)
                else
                    result = result + (amat[pos,j] - alpha) 
            }
        }
    }
    ConsensusFnChange = result
}