h1=plot(cubicdelay5thttp0500ms(1:180,3),cubicdelay5thttp0500ms(1:180,2)/1000,'--');

hold on

h2=plot(cubicdelay5thttp0500ms(1:180,3),theorapiddelay5thttp0500ms(1:180,2)/1000);

h4=plot(cubicdelay5thttp0500ms(1:180,3),gap200tau10delay5thttp0500ms(1:180,2)/1000,'o-');

[h1,h2,h4]=legend('Cubic','Rapid w/0 BASS','Rapid with BASS');

plot(cubicdelay5thttp0500ms(1:180,3),gap200tau10delay5thttp0500ms(1:180,2)/1000);

h3=plot(cubicdelay5thttp0500ms(1:180,3),gap200tau10delay5thttp0500ms(1:180,2)/1000,'o');

xdata=get(h3,'XData');

ydata=get(h3,'YData');

set(h3,'XData',xdata(1:15:end));

set(h3,'YData',ydata(1:15:end));

xlabel('Time (s)','FontSize',16)

ylabel('Throughput (Gbps)','FontSize',16);

set(gca,'FontSize',16)

set(gcf,'PaperUnits','points')

set(gcf,'PaperSize',[680,480])

set(gcf,'PaperPositionMode','manual')

set(gcf,'PaperPosition',[0 0 680 480])

print(gcf,'-dpdf','protocol_bct_rtt5ms.pdf')