[zip city state lat lon zone dst] = textread('zipcode.txt','%s%s%s%f%f%d%d','headerlines',1);

n  = length(zip);	
figure(1);
for i = 1:100:(n-100) 
    c = [max((-n+2*i-1.0)/n,0) min((2*i-2.0)/n,(2*n-2*i+2.0)/n) max((n-2*i+1.0)/n,0)];
%    plot(lon(i:i+101),lat(i:i+101),'Color',c, 'LineStyle','-');
    plot(lon(i:i+101),lat(i:i+101),'Marker','o','MarkerFaceColor',c, ...
		    'MarkerEdgeColor','none','MarkerSize',2,'LineStyle','none');
    hold on;
end
hold off;

figure(2);
for i = 1:100:(n-100) 
    xx = sin((90-lat)*pi/180).*cos((lon)*pi/180);
    yy = sin((90-lat)*pi/180).*sin((lon)*pi/180);
    zz = cos((90-lat)*pi/180);
    c = [(i-1.0)/n 0.0 (n-i+1.0)/n];
%    [X,Y,Z] = sphere(100);
%    surf(X,Y,Z); hold on;
    sphere(100);
%    plot3(xx(i:i+101),yy(i:i+101),zz(i:i+101),'-','Color',c);
    plot3(xx(i:i+101),yy(i:i+101),zz(i:i+101), 'Marker','o','MarkerFaceColor',c, ...
		    'MarkerEdgeColor','none','MarkerSize',2,'LineStyle','none');
    hold on;
end
hold off;
