cy = math.sqrt(M[i, i]*M[i, i] + M[j, i]*M[j, i])