of magnitude. Under Eq.  the Higgs field $h$ in Eq.  splits into the familiar MSSM doublets $h_u$ and $h_d$ whose neutral components subsequently develop weak scale VEVs $\langle h_u^0 \rangle = v_2$ and $\langle h_d^0 \rangle = v_1$ with $\tanb = v_2/v_1$. The neutrino fields $\nu^c$ acquires a large mass $M_{RR} \sim \lambda' \langle HH \rangle / M_V$ through the non-renormalizable term in ${\cal W}$ which, together with the Dirac $\nu^c$ – $\nu$ interaction (proportional to $\lambda \langle h_u^0 \rangle$), gives rise to a 2 $\times$ 2 matrix that generates, via see-saw mechanism [@seesaw], a suppressed mass for the left-handed neutrino states. The $D$ field does not develop a VEV but the terms $HHD$ and $\bar H \bar H D$ combine the colour triplet parts of $H$, $\bar H$ and $D$ into acceptable GUT scale mass terms [@AnLe].

The pattern of fermion masses and mixing angles is one of the fundamental problems is particle physics that has not yet been understood. The importance of this unsolved puzzle is demonstrated by the numerous works published in the literature over the past years (see Refs. [@MatrixModels1]-[@BaDo] for a “short” list.) In the standard model (SM) the quark/lepton masses and the CKM matrix 