function position = NewtonsMethodGPS(satellites, initial_guess, steps)
position = initial_guess;
for i = 1 : steps
    J = SatellitesJacobianMatrix(satellites, position);
    R = EvaluateSatelliteEquations(satellites, position);
    Jt = J';
    J = Jt*J;
    R = Jt*R;
    v = J\R;
    for j = 1 : 4
        position(j) = position(j) - v(j);
    end
end

