In this problem we were expected to add on additional satellites to our equations.
Normally this would require a combined least squares/multivariate Newton's method.
However, Matlab, when given more equations than variables, already expects this and does the least squares method for us.
To accomplish this we modified our original f and dF matrix to be constructed within a for loop allowing us to dynamically establish the number of satellites used.
For our data we added 3 more satellites for a total of 7 unbunched satellites.
We received a position error of 2.5 meters and an error magnification factor of .835
Code - Gauss Function
Code - Polar/Init
Output