1. #include<stdio.h>  
    2.  
         
  2. #include<math.h>  
    3.  
         
  3. #define PI  3.1415926  
    4.  
         
  4. #define R   6371.137   
    5.  
         
  5. #define RAD PI/180.0  
    6.  
         
  6.  
    7.  
         
  7. int main()           
    8.  
         
  8. {                   
    9.  
         
  9.         float   a1,b1,a2,b2;  
    10.  
         
  10.     float   x1,y1,A1,B1;   
    11.  
         
  11.     float s;  
    12.  
         
  12.         printf("please input first longitude and latitude\n");  
    13.  
         
  13.      scanf("%f%f",&a1,&b1);  
    14.  
         
  14.      printf("plsease input second logitude and latitude\n");  
    15.  
         
  15.      scanf("%f%f",&a2,&b2);  
    16.  
         
  16.  
    17.  
         
  17.      A1=b1 * RAD;  
    18.  
         
  18.      B1=b2 * RAD;  
    19.  
         
  19.  
    20.  
         
  20.      x1 = A1 - B1;  
    21.  
         
  21.      y1 = (a1-a2)*RAD;   
    22.  
         
  22.  
    23.  
         
  23.      s=2*asin(sqrt(pow(sin(x1/2),2) +   
    24.  
         
  24.             cos(b1)*cos(b2)*pow(sin(y1/2),2)));  
    25.  
         
  25.      s *= R;  
    26.  
         
  26.      printf("The distance is %f\n",s);  
    27.  
         
  27. } 
    28.