1:
2:
3:
4:
5:
6:
7:

i1 = -3:0.1:3;
f1 = -3:0.1:3;
[I1, F1] = meshgrid(i1, f1);
x1 = cos(I1).*sin(F1);
y1 = sin(I1).*sin(F1);
z1 = cos(F1);
surf(x1,y1,z1) 

1:
2:
3:
4:
5:

a = 5;
x2=cos(I1).*sin(F1)+a*ones(length(I1));
surf(x1,y1,z1)
surf(x2,y1,z1)

1:
2:
3:
4:

r = 0.2;
n = 100
[A,B,C] = cylinder(r,n);
mesh(C,A,B)

1:

mesh(0.9:2.7,A,B)

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:

r = 0.2;
n = 100
[X,Y,Z] = cylinder(r,n);
mesh(Z,X,Y)


hold on
%---
i1 = -3:0.1:3;
f1 = -3:0.1:3;
[I1, F1] = meshgrid(i1, f1);
x1 = cos(I1).*sin(F1);
y1 = sin(I1).*sin(F1);
z1 = cos(F1);
surfl(x1,y1,z1) % erste Kugel

hold on
xlabel('x')
ylabel('y')
zlabel('z')
a =3; % Abstand zwischen den Kugeln
x2=cos(I1).*sin(F1)+a*ones(length(I1));
surfl(x2,y1,z1) % zweite Kugel
cm = hsv(1000);
cm(:,3)= 0 * cm(:,3);
cm = cm(500:1000,:);
colormap(cm);
shading interp;

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:

a = input('Geben sie den Abstand der Kugeln an: '); % Abstand zwischen den Kugeln
i1 = -3:0.1:3;
f1 = -3:0.1:3;
[I1, F1] = meshgrid(i1, f1);
x1 = cos(I1).*sin(F1);
y1 = sin(I1).*sin(F1);
z1 = cos(F1);

surf(x1,y1,z1) % erste Kugel

hold on
xlabel('x')
ylabel('y')
zlabel('z')

x2=cos(I1).*sin(F1)+a*ones(length(I1));
surf(x2,y1,z1) % zweite Kugel

% der Radius der Kugeln ist ja max 1 da mit sin und cos gerechnet wird
xa = 0.9; % Gewählter Wert für den Beginn des Zylinders da nah an 1
xe = 0.9*(1+a); % berechneter Endwert des Zylinders
x = xa:0.01:xe; % Bereich in dem der Zylinder geplottet werden soll
t=linspace(0,2*pi,length(x));
[X,T]=meshgrid(x,t);
Y=0.2*cos(T); % 0.2 als Radius des Zylinders 20% des Kugel Radius
Z=0.2*sin(T);
mesh(X,Y,Z)
hold on
set(gcf,'color','white')
hold on
farbe = gray(1000);
colormap(farbe);
shading interp;