clear all

results=[];
data=[];

imax=400;           % maximum set size
for ii=1:10,
densn=.05/ii;
%rng('default');        % reset random number generator to default at start
rng('shuffle');        % reset random number generator using current time
k=randi([65 70],1);     % k is number of sets
%[Aprojeigbnd,upperAbnd,n,Ecard,cputimeprojeigbnd,cputimeupper]=...
%                        testsparseprojeigbnd(imax,k,densn);
[Aprojeigbnd,upperAbnd,n,Ecard,cputimeprojeigbnd,cputimeupper]=...
                        newV0sparsetest(imax,k,densn);
data=[data
     n k Ecard Ecard/( n*(n-1)/2 ) ]; % density fixed by TK, Jan 13, 2014
   disp('ii ceil(Aprojeigbnd)  floor(upperAbnd) cputimeprojeigbnd cputimeupper')
   tempd=[ii ceil(Aprojeigbnd) floor(upperAbnd) cputimeprojeigbnd cputimeupper];
   disp([num2str(tempd)]);
results=[results
     ceil(Aprojeigbnd)  floor(upperAbnd) cputimeprojeigbnd cputimeupper];
end
relerror=(floor(results(:,2))-ceil(results(:,1)))./(floor(results(:,2))+ceil(results(:,1)));

textdata=latextable([data(:,1:4) ceil(results(:,1))  floor(results(:,2)) relerror results(:,[3 4])
],'%i','%i','%i','%3.2e','%i','%i','%5.2f','%5.1f','%5.1f'); % display changed by TK, Jan 13, 2014

save resultV0 data results relerror textdata
