% PKPD
function [p] = method3(r)
k=size(r,1);
C = r + eye(k);
C = 1./C - eye(k);
p = 1./(sum(C')'-(k-2));
p = p/sum(p);