% by product!!

function [p] = method5(r)
k=size(r,1);
rr = r - diag(diag(r))+ eye(k);
p = prod(rr,2);
p = p/sum(p);
