function y = EvaluateSatelliteEquations(satellites, position)
speed_of_light = 299792.458;
speed_of_light_squared = speed_of_light^2;

[rows, columns] = size(satellites);
y = zeros(rows, 1);

for i = 1 : rows
    left_side = 0;
    for j = 1 : 3
        left_side = left_side + (position(j) - satellites(i, j))^2;
    end
    right_side = (speed_of_light*(satellites(i, 4) - position(4)))^2;
    y(i) = left_side - right_side;
end

