function Result=SatelliteConditionNumber(NumberOfSatellites)
Theta=linspace(0,2*pi,NumberOfSatellites+2);
Phi=linspace(0,pi/2,NumberOfSatellites+2);
Rho=26570.*ones(1,NumberOfSatellites);

Theta=Theta(2:NumberOfSatellites+1);
Phi=Phi(2:NumberOfSatellites+1);

A=Rho.*cos(Phi).*cos(Theta);
B=Rho.*cos(Phi).*sin(Theta);
C=Rho.*sin(Phi);

%c is light speed
c=299792.458;

InitialGuess=[0,0,6370,0.0001];

R=sqrt(A.^2+B.^2+(C-6370).^2)
T=InitialGuess(4)-R./c
Result=MVNewtonsInputs(InitialGuess,A,B,C,T);

