function Z2 = Z2_min(AR, D, g1, f1, r1)

[R1, S1, S2, d1, g2, a1l, a1r] = common(AR, D, g1, f1);
fx = fx_max(AR, D, g1, f1, r1);
da1 = atan2(R1, S1 + S2) - atan2(R1 - fx, S1 + S2);
Da1 = 2 * a1r;
Z2 = Da1 / da1;

return;