function plot_freq_2(AR, D, g1, f1)

x = 0.001*(100:-1:1)';
f2 = x;
fc = x;
fL1 = x;
fL2 = x;

for i = 1:size(x, 1),
  f2(i) = freq_2(AR, D, g1, f1, x(i));
  fc(i) = freq_RGB(AR, D, g1, f1, x(i));
  fL1(i) = L1(AR, D, g1, f1, x(i));
  fL2(i) = L2(AR, D, g1, f1, x(i));
end

write('plot_freq_2.xy', [x, f2]);
write('plot_freq_RGB.xy', [x, fc]);
write('plot_L1.xy', [x, fL1]);
write('plot_L2.xy', [x, fL2]);

return;