function d = objectiveFunc(vars, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3)
    X = vars(1);
    Y = vars(2);
    Z = vars(3);
    
    % Compute squared distances from (X, Y, Z) to each circle
    d1 = min((X - X1).^2 + (Y - Y1).^2 + (Z - Z1).^2);
    d2 = min((X - X2).^2 + (Y - Y2).^2 + (Z - Z2).^2);
    d3 = min((X - X3).^2 + (Y - Y3).^2 + (Z - Z3).^2);
    
    % Sum of squared distances
    d = d1 + d2 + d3;
end