function tuple = metatest(numbersats,iterations,phi,theta)
maxerror = 0;
maxsats = [];
X = [];
Y = [];
for i = 1:iterations
    M = [ones(numbersats,1)*26570 (pi*phi/2).*rand(numbersats,1)+(pi/2)*rand (2*pi*theta).*rand(numbersats,1)+(2*pi)*rand];
    error = testsats(M);
    X = [X distance(M)];
    Y = [Y error];
    if error > maxerror
        maxsats = M;
        maxerror = error;
    end
end
%semilogy(X,Y,'o')
tuple = Y;