%Hosszusagi kor szogek (rekta aszcenczio):
rac = 0:pi/64:2*pi;
%Szelessegi kor szogek (deklanacio):
dec = 0:pi/64:pi;
%!Figyelem ez csak a felso terfelt rajzolja meg!
%!Ha a teljesre szukseg van, akkor dec = -pi:pi/64:pi;!

%Kezd
X = [];
Y = [];
Z = [];
k = 1;
%Szamolas
for i = 1:length(rac),
for j = 1:length(dec),
%Ha pl. r = 1 a fuggveny, akkor egy gombot rajzol
%Ez egy pelda fuggveny: r = 2*sin(rac)^2 + cos(dec)^3
r = 2*sin(rac(i))^2 + cos(dec(j))^3;
%konvertalas
X(k) = r * cos(rac(i))*cos(dec(j));
Y(k) = r * sin(rac(i))*cos(dec(j));
Z(k) = r * sin(dec(j));
%lepes
k = k+1;
end
end
%Rajzolas
plot3(X,Y,Z,'.');