function r = p6()

c = 299792.458;  %(km)
dt = 1e-8;

% Some Possible ways of arranging -1e-8 and 1e-8 for the four Ti's
dt_combinations = [-dt dt -dt dt -dt dt -dt dt -dt;
                   dt -dt dt -dt dt -dt dt -dt dt;
                   dt dt dt dt -dt dt -dt dt -dt;
                   -dt dt -dt dt -dt dt dt dt dt;
                   -dt -dt -dt dt dt dt -dt -dt dt;
                   dt -dt -dt -dt -dt dt dt dt -dt;
                   -dt -dt -dt -dt -dt dt -dt dt -dt;
                   -dt dt -dt dt -dt -dt -dt -dt -dt;
                   dt dt dt -dt -dt -dt dt dt dt;
                   -dt -dt -dt dt dt dt -dt -dt -dt;
                   -dt dt dt dt dt dt -dt dt -dt;
                   dt dt -dt dt dt dt dt dt dt;
                   dt dt dt dt -dt dt dt dt dt;
                   dt dt -dt dt dt dt dt dt -dt];
               
% Gauss-Newton Solution without changing t's
x0 = gauss_newton_p6(3, 0);

max_position_error = 0;
max_EMF = 0;

max_tdiff = 1e-8;  % Maximum change in t
BE = c*max_tdiff;  % Backward Error


for i = 1:14
   x(i, 1:4) = gauss_newton_p6(3, 1, dt_combinations(i, 1:9));   % Run dt adjusted Gauss-Newton
   xdiff(i, 1:4) = abs(x(i, 1:4)' - x0);  % Difference between new x's and original
   max_sdiff(i) = max(xdiff(i, 1:4));
   max_position_error = max(max_position_error, max_sdiff(i));
   FE(i) = max_sdiff(i);            % Forward Error
   EMF(i) = FE(i)/BE;               % Get error magnification factor
   max_EMF = max(max_EMF, EMF(i));  % Find max EMF (cond number estimate)
end

max_position_error
condition_number_estimate = max_EMF

r = condition_number_estimate;