load nfft_times_frodo_m_2_flops.m
load nfft_times_frodo_flops.m

A=nfft_times_frodo_flops;
B=nfft_times_frodo_m_2_flops;
n1=A(1:20,1);
r1=A(1:20,2);
s1=A(1:20,4);
t1=B(1:20,4);

n2=A(21:29,1);
r2=A(21:29,2);
s2=A(21:29,4);
t2=B(21:29,4);

n3=A(30:34,1);
r3=A(30:34,2);
s3=A(30:34,4);
t3=B(30:34,4);

h=plot(n1,5*2.^n1.*n1 ./ r1*10^-6, 'k-', ...
     n1,5*2.^n1.*n1 ./ r1*10^-6, 'ko', ...
     n1,(2.^n1 + 5*2.^(n1+1).*(n1+1) + (4*9+0)*2.^n1) ./ s1*10^-6, 'k-', ...
     n1,(2.^n1 + 5*2.^(n1+1).*(n1+1) + (4*9+0)*2.^n1) ./ s1*10^-6, 'kx', ...
     n1,(2.^n1 + 5*2.^(n1+1).*(n1+1) + (4*5+0)*2.^n1) ./ t1*10^-6, 'k-', ...
     n1,(2.^n1 + 5*2.^(n1+1).*(n1+1) + (4*5+0)*2.^n1) ./ t1*10^-6, 'k+', ...
     n2,5*2.^n2.*n2 ./ r2*10^-6, 'k:', ...
     n2,5*2.^n2.*n2 ./ r2*10^-6, 'ko', ...
     n2,(2.^n2 + 5*2.^(n2+2).*(n2+2) + (4*9^2+0)*2.^n2) ./ s2*10^-6, 'k:', ...
     n2,(2.^n2 + 5*2.^(n2+2).*(n2+2) + (4*9^2+0)*2.^n2) ./ s2*10^-6, 'kx', ...
     n2,(2.^n2 + 5*2.^(n2+2).*(n2+2) + (4*5^2+0)*2.^n2) ./ t2*10^-6, 'k:', ...
     n2,(2.^n2 + 5*2.^(n2+2).*(n2+2) + (4*5^2+0)*2.^n2) ./ t2*10^-6, 'k+', ...
     n3,5*2.^n3.*n3 ./ r3*10^-6, 'k--', ...
     n3,5*2.^n3.*n3 ./ r3*10^-6, 'ko', ...
     n3,(2.^n3 + 5*2.^(n3+3).*(n3+3) + (4*9^3+0)*2.^n3) ./ s3*10^-6, 'k--', ...
     n3,(2.^n3 + 5*2.^(n3+3).*(n3+3) + (4*9^3+0)*2.^n3) ./ s3*10^-6, 'kx', ...
     n3,(2.^n3 + 5*2.^(n3+3).*(n3+3) + (4*5^3+0)*2.^n3) ./ t3*10^-6, 'k--', ...
     n3,(2.^n3 + 5*2.^(n3+3).*(n3+3) + (4*5^3+0)*2.^n3) ./ t3*10^-6, 'k+');

axis([3,22,0,2400]);
set(gcf,'Position',[100,100,1000,333]);
set(h,'LineWidth',1.8/3); set(h,'Markersize',10/3); set(gca,'FontSize',20/3); 
set(gca,'DataAspectRatio',[0.25*9.5 1200 1]);
%set(gca,'XTickLabel',[-n,-n+N/2,-N/2,0,N/2-1,n-N/2-1,n]);
%set(gca,'YTick',[10^-15,10^-10,10^-5,1]);
print nfft_times_flops.eps -deps
!epstopdf nfft_times_flops.eps
