function N = N(a, e2, fi)
% Pricny polomer krivosti
% N = N(a, e2, fi)

W = sqrt( 1 - e2*( sin(fi) )^2 );

N = a / W;