plot two tangent circles in matlab

main.m
% Define the centers and radii of the two circles
center1 = [0 0];
radius1 = 1;
center2 = [2 0];
radius2 = 1;

% Calculate the distance between the centers
d = norm(center2 - center1);

% Calculate the intersection points of the circles
a = (radius1^2 - radius2^2 + d^2) / (2*d);
h = sqrt(radius1^2 - a^2);
p = center1 + a * (center2 - center1) / d;
intersect1 = p + h * [0 1] * (center2 - center1) / d;
intersect2 = p - h * [0 1] * (center2 - center1) / d;

% Plot the circles and intersection points
figure;
viscircles(center1, radius1, 'Color', 'b');
hold on;
viscircles(center2, radius2, 'Color', 'r');
plot(intersect1(1), intersect1(2), 'ro');
plot(intersect2(1), intersect2(2), 'ro');
axis equal;
703 chars
25 lines