
RTT=0.05;
b=1;

%To=0.1;
%L=1;
%d=2;

x=0.001:0.001:1;

%for n=1:5:101
%for ind=0:1
p=zeros(1000,1);
fid=fopen('data-fig4.txt','w');
for p=1:1:1000
  y1(p,1)=eqtcp(RTT,b,p/1000,1,1,2,0);
  y2(p,1)=eqtcp(RTT,b,p/1000,1,1,2,1);
  y3(p,1)=eqtcp(RTT,b,p/1000,1,4,2,2);
  y4(p,1)=eqtcp(RTT,b,p/1000,1,1,1.25,2);
  y5(p,1)=eqtcp(RTT,b,p/1000,0.1,1,2,2);
  y6(p,1)=eqtcp(RTT,b,p/1000,0.1,4,1.25,2);    
  fprintf(fid,'%f %f %f %f %f %f %f \n', p/1000, y1(p,1), y2(p,1), y3(p,1), y4(p,1), y5(p,1), y6(p,1));  
 end
fclose (fid);

