Mahalanobis fuzzy Xie-Beni index

Input: data set 𝒜 ⊂ (ℝ^n), centers (Subscript[z, 1]),⃛,(Subscript[z, k]) ∈ (ℝ^n), membership matrix U, covariance matrices (∑_j)∈(ℝ^(n×n)), j=1⃛k;
Output: MFXB index;
Submodules: dM[]

Clear[MFXB]
MFXB[A_,cen_,U_,cov_]:=Module[{m,k},m=Length[A];k=Length[cen];Sum[U[[i,j]]^q dM[cen[[j]],A[[i]],cov[[j]]],{i,m},{j,k}]/(m*Min[Flatten[Table[dM[cen[[i]],cen[[j]],cov[[i]]+cov[[j]]],{i,k-1},{j,i+1,k}],1]])
]