//calculate gammas
Gamma1 = ((S12Rare * t21) / (float)((C12 * t10 * t11)) - 1);
Gamma2 = ((S12Rare * t12) / (float)((C12 * t01 * t11)) - 1);
Gamma3 = ((S12Rare / C12) * (S12Rare / C12)) * ( t22 / (float)(t10 * t01 * t11));
Gamma3 = Gamma3 - ((S12Rare * t11) / (float)(C12 * t01 * t10)) - Gamma1 - Gamma2;
//calculate gammas
Gamma1 = ((S12Rare * t21) / (float)((C12 * t10 * t11)) - 1);
Gamma2 = ((S12Rare * t12) / (float)((C12 * t01 * t11)) - 1);
Gamma3 = ((S12Rare / C12) * (S12Rare / C12)) * ( t22 / (float)(t10 * t01 * t11));
Gamma3 = Gamma3 - ((S12Rare * t11) / (float)(C12 * t01 * t10)) - Gamma1 - Gamma2;