%% ‘æ‚T�}

sys1 = rperf0_clsys(Np0+1:end,Nm0+1:end);
sys2 = rperf_clsys(Np0+1:end,Nm0+1:end);

figure;
bodemag(sys1, 'k:', sys2, 'k-');
set(gcf,'PaperUnits','inches','PaperPosition',[0 0 6 8])
title('Bode diagram');
xlabel('frequency');
ylabel('magnitude');
print -depsc2 -r200 bodegains.eps 


