h=R(1-sqrt(1+cos(2*arcsin(1-sqrt(c/(2*R))))/2));R=6371,c=50