H=diag(p)-(pp^T),where   (Subscript[p, i])=((Subscript[f, i])/((∑_(i=1)^n) (Subscript[f, i])))