function M = M(a, e2, fi)
% Meridianovy polomer krivosti
% M = M(a, e2, fi)

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

M = a*(1-e2) / W^3;