| 36 | } |
| 37 | } |
| 38 | static void do_arc(void) { |
| 39 | double az; |
| 40 | |
| 41 | printLL(phi2, lam2); |
| 42 | putchar('\n'); |
| 43 | for (az = al12; n_alpha--;) { |
| 44 | al12 = az = adjlon(az + del_alpha); |
| 45 | geod_pre(); |
| 46 | geod_for(); |
| 47 | printLL(phi2, lam2); |
| 48 | putchar('\n'); |
| 49 | } |
| 50 | } |
| 51 | static void /* generate intermediate geodesic coordinates */ |
| 52 | do_geod(void) { |
| 53 | double phil, laml, del_S; |