Membership function for Mahalanobis fuzzy clustering

Input: data set 𝒜 ⊂ (ℝ^n), centers (Subscript[c, 1]),⃛,(Subscript[c, k]) ∈ (ℝ^n), covariance matrices (∑_t)∈(ℝ^(n×n)), t=1⃛k, indexes i,j;
Output: membership function (Subscript[u, ij]);
Submodules: dM[]

q=2;
muGK[A_,cen_,Scov_,i_,j_]:=1/(Sum[(dM[A[[i]],cen[[j]],Scov[[j]]]/dM[A[[i]],cen[[s]],Scov[[s]]])^(1/(q-1)),{s,Length[cen]}])