change in function eta, it accepts a third optional argument. Change in i2sm1m2_c, the products x1*x2 and (1-x1)*(1-x2) are computed explicitly