- pl = pl1 * (1 - pe) + pl * (1 - pe1) - pl * pl1;
- pe = pe + (1 - pe) * pe1;
- n2 = hv_iterinit(res[n]);
+ pl = pl1 * (1 - (pd + pe)) + pl * (1 - (pd1 + pe1)) - pl * pl1;
+ pd2 = pd1 * (1 - pe) + pd - pd * pd1;
+ pe = pe1 * (1 - pd) + pe - pe * pe1;
+ pd = pd2;
+ n2 = hv_iterinit(res[n]);