e = 10^-8;
actual = [0, 0, 6370]';
errorlist = [];
for i = -1:1
	for j = -1:1
		for k = -1:1
			for l = -1:1
				guess = [0, 0, 0, 0]';
				for m = 1:10
					guess = timeerror4(guess(1), guess(2), guess(3),...
						  	guess(4), i*e, j*e, k*e, l*e);
				end
				[i*e, j*e, k*e, l*e];
				spaceguess = [guess(1), guess(2), guess(3)]';
				guesserror = norm((actual - spaceguess), Inf);
				errorlist(end + 1) = guesserror;
			end
		end
	end
end
maxerror = norm(errorlist, Inf)

maxerror =

    0.0199

Published with MATLAB® R2013a