format long
a = [16.641038888888888];     % 1 component 
b = [16 38.462333333333333];  % 2 components  
c = [16 38 27.74];            % 3 components 
d = [16 38 27 74];            % 4 components
radians = [deg2rad(a) deg2rad(b) deg2rad(c) deg2rad(d)]


f = deg2rad([49 3 37 2]); 
l = deg2rad([16 38 27 74]); 
BRNO = gg2gc(f,l)           


format long                              
f1 = deg2rad([49 45 18]);                % Praha 
l1 = deg2rad([14 7 58]);                 %   position
f2 = deg2rad([49 3 37 2]);               % Brno
l2 = deg2rad([16 38 27 74]);             %   position
A1praha = [200; 140; -5.790];            % 1. A. seen from Praha 
A2praha = [157; 397; -9.870];            % 2. A. seen from Praha 
A1brno = transf(f1, l1, A1praha, f2, l2) % 1. A. seen from Brno

a1praha = transf(f2, l2, A1brno, f1, l1) % inv. transformation

A2brno = transf(f1, l1, A2praha, f2, l2) % 2. A. seen from Brno 

a2praha = transf(f2, l2, A2brno, f1, l1) % inv. transformation


