function fy = fy_max(AR, D, g1, f1, r1)

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

fy = 2 * r1 / sin(2 * g1);

return;