function [ a, b, e2 ] = bessel()
% parametry Besselova elipsoidu
% [ a, b, e2 ] = bessel()

a = 6377397.15508;
b = 6356078.96290;
e2 = ( (a^2-b^2)/a^2 );