Kula działki z Julią i PyPlot

Question

Ostatnio próbowałem wykreślić kulę za pomocą PyPlot/Julia i niestety było to trudniejsze, niż myślałem. Prawdopodobnie jest coś nie tak z generowaniem punktów, ale nie mogę zrozumieć, dlaczego moja implementacja nie zadziałała. Chociaż wszystko jest w porządku z oryginalnym kodem Pythona.

Starałem się dostosować demo2 z matplotlib powierzchni działki doc jak MWE:

using PyPlot 
u = linspace(0,2*π,100); 
v = linspace(0,π,100); 

x = cos(u).*sin(v); 
y = sin(u).*sin(v); 
z = cos(v); 

surf(x,y,z) 

A ja dostaję zamiast

Tak, to jest dokładnie to źle w moim realizacji Julia?

Answer 1

x, y i z powinny być macierzami, nie wektorami - w przeciwnym razie na powierzchni sfery zamiast samej powierzchni rysuje się tylko krzywa.

using PyPlot 
n = 100 
u = linspace(0,2*π,n); 
v = linspace(0,π,n); 

x = cos(u) * sin(v)'; 
y = sin(u) * sin(v)'; 
z = ones(n) * cos(v)'; 

# The rstride and cstride arguments default to 10 
surf(x,y,z, rstride=4, cstride=4) 

Narysowana początkowo krzywa odpowiada przekątnej tych macierzy.

plot(diag(x), diag(y), diag(z), color="yellow", linewidth=3)