fixfigure(1);
plot(f, 20*log10(snonoise+eps), 'b');  hold on;
plot(f, 20*log10(sn1+eps), 'r');  hold on;
plot(f, 20*log10(sn2+eps), 'm');
axis([0 50e3 -150 50]);
xlabel('Frequency  [Hz]');
ylabel('Output Spectrum [dBWN]  /  Int. Noise [dBV]');
legend('No noise', '1st Integrator', '2nd Integrator');
grid on;  hold off;

fixfigure(2);
plot(f, 20*log10(snonoise+eps), 'b');  hold on;
plot(f, 20*log10(sn1+eps), 'r');  hold on;
plot(f, 20*log10(sn2+eps), 'm');
axis([0.49*fs 0.5*fs -150 50]);
xlabel('Frequency  [Hz]');
ylabel('Output Spectrum [dBWN]  /  Int. Noise [dBV]');
legend('No noise', '1st Integrator', '2nd Integrator', 3);
grid on;  hold off;
