% earth.m  
% mapping latitude and longitude to 3D

longitude_range = linspace(0,2*pi,50);
latitude_range = linspace(-pi/2,pi/2,50);
[longitude, latitude] = meshgrid(longitude_range,latitude_range);
X = cos(latitude) .* cos(longitude);
Y = cos(latitude) .* sin(longitude);
Z = sin(latitude);
surfl(X,Y,Z)
axis equal