(1/(k-1))((Subscript[Pe, n]^(k,i))+(Subscript[Pe, n]^(k,j))-(Subscript[Pe, n]^(1,i))-(Subscript[Pe, n]^(1,j))),