#include <stdio.h>
#include <math.h>

int main() {
        float V1, V2, V3, TV1, TV2, PI, A, B, R;

	// case 1
	V1 = -7.10521984100341796875000000000000000f;
	V2 = -0.40112295746803283691406250000000000f;
	V3 = 2.0f;
	PI = 3.141592f;

	TV1 = -tanf(V1*(PI/180.0f));
	TV2 = tan(V2);

	A = acosf(TV1*TV2);
	B = 1.0f+0.033f*cosf(2.0f*PI*V3/365.0f);

	R = 15.392f*B*(A*sinf(V1*(PI/180.0f))*sinf(V2)+cosf(V1*(PI/180.0f))*cosf(V2)*sinf(A));

	printf("Case 1:\n");
        printf("V1:  %40.35f\n", V1);
        printf("V2:  %40.35f\n", V2);
        printf("TV1: %40.35f\n", TV1);
        printf("TV2: %40.35f\n", TV2);
        printf("A:   %40.35f\n", A);
        printf("B:   %40.35f\n", B);
        printf("R:   %40.35f\n", R);

	// case 2
	V1 = -10.81274986267089843750000000000000000f;
	V2 = -0.40112295746803283691406250000000000f;
	V3 = 2.0f;
	PI = 3.141592f;

	TV1 = -tanf(V1*(PI/180.0f));
	TV2 = tan(V2);

	A = acosf(TV1*TV2);
	B = 1.0f+0.033f*cosf(2.0f*PI*V3/365.0f);

	R = 15.392f*B*(A*sinf(V1*(PI/180.0f))*sinf(V2)+cosf(V1*(PI/180.0f))*cosf(V2)*sinf(A));

	printf("Case 2:\n");
        printf("V1:  %40.35f\n", V1);
        printf("V2:  %40.35f\n", V2);
        printf("TV1: %40.35f\n", TV1);
        printf("TV2: %40.35f\n", TV2);
        printf("A:   %40.35f\n", A);
        printf("B:   %40.35f\n", B);
        printf("R:   %40.35f\n", R);

        return(0);
}