clc
clear

Ro=180/pi;   %radian
R=6380;


u2 = [10/Ro; 20/Ro; 30/Ro; 40/Ro];
v2 = [10/Ro; 15/Ro; 18/Ro; 20/Ro];




%GRAFIKA-----------------------------------------------------------------

%kresli rovnobezky
for lat = -90/Ro:10/Ro:90/Ro %beru rovnobezky po 10 stupnich
    
    %pro kazdou rovnobezku hledam pruseciky s poledniky po 10 stupnich
    v = -pi : 10/Ro : pi; %vektor 37 souradnic zem. delek pruseciku

    %velikost vektoru lon (1,37)
    [m,n] = size(v)
    
    %vektor jednicek (1,37) vynasobeny sirkou rovnobezky
    u = ones(m,n)*lat;
    
    %spocita souradnice vsech 37 pruseciku v danem zobrazeni 
    x = R*sin(v).*cos(u)
    y = R*u
    
    %nakresli graf (spoji 37 pruseciku, modre)
    plot (x, y, 'b');
    
    %graf zustane otevren pro dalsi kresleni
    hold on   
end


%kresli krivku

x2 = R*sin(v2).*cos(u2)
y2 = R*u2
    
%nakresli graf (spoji pruseciky, cervene)
plot (x2, y2, 'r');
    
%graf zustane otevren pro dalsi kresleni
hold on   

    
    

    
    

