a1 *= ((d1 *(e1 + f1)+ e1 * f1) + c1 *((e1 + f1) + d1))
             + b1 *(((e1 + f1) + d1) + c1);
         b1 *= (d1 *(e1 + f1)+ e1 * f1) + c1 *((e1 + f1) + d1);