%sup de velocidad 0 en 3D
close all
C=165;

z=[0:1e3:2*GEO(1)];

for k=1:length(z)
    Surff3(C,GEO,mu,wt,z(k));
end
[xE, yE, zE] = ellipsoid(0,0,0,GEO(1),GEO(2),GEO(3),30);
hold on
surf(xE, yE, zE)