function R = R( a, e2, fi )
% Strední polomer krivosti
% R = R( a, e2,fi )

R = sqrt(M(a, e2, fi) * N(a, e2, fi));