function fRGB = freq_RGB(AR, D, g1, f1, r1)

[R1, S1, S2, d1, g2, a1l, a1r, w1] = common(AR, D, g1, f1);
Z1 = Z1_min(AR, D, g1, f1, r1);
Z2 = Z2_min(AR, D, g1, f1, r1);
fRGB = pi * Z1 * Z2 * w1 / (atan2(R1, S1 + S2) * (g2 - g1));

return;