function b1 = beta1(AR, D, g1, f1, r1)

[R1, S1, S2, d1, g2, a1l, a1r] = common(AR, D, g1, f1);

b1 = atan2(d1 * (S1 + S2), R1 * S2);

return;