Mahalanobis fuzzy  Calinski–Harabasz 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: MFCH index;
Submodules: dM[]

Clear[MFCH]
MFCH[A_,cen_,U_,cov_]:=Module[{m,k},m=Length[A];k=Length[cen];
cc=Mean[A];
covcc=Sum[KroneckerProduct[(A[[i]]-cc),(A[[i]]-cc)],{i,m}]/m;
Return[(Sum[Sum[U[[ii,j]]^q,{ii,m}] *dM[cc,cen[[j]],cov[[j]]+covcc],{j,k}]/(k-1))/(Sum[U[[i,j]]^q dM[cen[[j]],A[[i]],cov[[j]]],{i,m},{j,k}]/(m-k))]
]