find the intersection of 3 xyz functions in matlab

main.m
% Define the functions
f1 = @(x, y, z) x.^2 + y.^2 + z.^2 - 25;
f2 = @(x, y, z) (x - 5).^2 + (y - 5).^2 + z.^2 - 25;
f3 = @(x, y, z) (x + 5).^2 + (y + 5).^2 + z.^2 - 25;

% Create a grid of values for x, y, and z
[x, y, z] = meshgrid(-10:0.1:10, -10:0.1:10, -10:0.1:10);

% Evaluate the functions on the grid
v1 = f1(x, y, z);
v2 = f2(x, y, z);
v3 = f3(x, y, z);

% Plot the surfaces
figure;
isosurface(x, y, z, v1, 0, 'r');
hold on;
isosurface(x, y, z, v2, 0, 'g');
isosurface(x, y, z, v3, 0, 'b');
xlabel('X');
ylabel('Y');
zlabel('Z');
title('Intersection of Three XYZ Functions');
grid on;

% Find the points of intersection
[xi, yi, zi] = spherefit(x(v1<0), y(v1<0), z(v1<0), 2);
plot3(xi, yi, zi, 'ro', 'MarkerSize', 10);
[xi, yi, zi] = spherefit(x(v2<0), y(v2<0), z(v2<0), 2);
plot3(xi, yi, zi, 'go', 'MarkerSize', 10);
[xi, yi, zi] = spherefit(x(v3<0), y(v3<0), z(v3<0), 2);
plot3(xi, yi, zi, 'bo', 'MarkerSize', 10);
legend('f1', 'f2', 'f3', 'Intersection Point');
974 chars
34 lines