--- /dev/null
+%load('pink.mat');
+%figure;
+
+s = 7400;
+max = t(length(t));
+
+ax = [0 max -5 5];
+ytick = -4:2:4;
+
+subplot(2,1,1);
+hold off;
+plot(t, pink2(s:(s+length(t)-1)), 'k');
+hold on;
+plot([0 max], [0 0], 'k:');
+axis(ax)
+set(gca,'XTick', 0:5:max)
+set(gca,'XTickLabel', [])
+set(gca,'YTick', ytick)
+ylabel('2 poles');
+
+%subplot(3,1,2);
+%hold off;
+%plot(t, pink3(s:(s+length(t)-1)), 'k');
+%hold on;
+%plot([0 max], [0 0], 'k:');
+%axis(ax)
+%set(gca,'XTick', 0:5:max)
+%set(gca,'XTickLabel', [])
+%set(gca,'YTick', ytick)
+%ylabel('5 poles');
+
+subplot(2,1,2);
+hold off;
+plot(t, pink10(s:(s+length(t)-1)), 'k');
+hold on;
+plot([0 max], [0 0], 'k:');
+axis(ax)
+set(gca,'XTick', 0:5:max)
+set(gca,'YTick', ytick)
+ylabel('10 poles');
+
+xlabel('Time / s');
+print('pinktimesub.eps')