function petaflops_LU

nV=[5 10 25 50 75 100 125 150];
sp=2*10^(17);

anno=24*60*60;

for k=1:length(nV)
    n=nV(k);
    ops(k,1)=factorial(n);
    ops(k,2)=n^3/3;
    secs=ops/sp;
end


for k=1:length(nV)/2
    if secs(k,1) > anno
        secs(k,1)=secs(k,1)/anno;
    end
    fprintf('\n \t $%3.0f$ & $%1.1e$ & $%1.1e$ & $%3.0f$ & $%1.1e$ & $%1.1e$ \\\\',...
        nV(k),secs(k,1),secs(k,2),nV(k+4),secs(k+4,1),secs(k+4,2));
end
