%Radius wird frei gewählt mit R = 1
R = 3;

[PHI, THETA] = meshgrid(linspace(0, 2*pi, 50), linspace(0, pi, 50));

X = R * sin(THETA) .* cos(PHI);
Y = R * sin(THETA) .* sin(PHI);
Z = R * cos(THETA);

%Darstellung der oberen Hälfte
%surf(X, Y, Z, EdgeColor="none");
fsurf('x^2')
hold on;
%Darstellung der unteren Hälfte
%surf(X, Y, Z, EdgeColor="none");

axis equal;%Gleichmäßige Achsenverhältnisse
xlabel("X");
ylabel("Y");
zlabel("Z");
title("Kugeloberfläche");

annotation("textarrow", [0.2, 0.5], [0.6, 0.5], String="Mitte des Bildes");