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

int main() {
    float V1, V2, V3, TV1, TV2, SV1, SV2, CV1, CV2, SA;
    float PI, A, B, R;

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

    TV1 = -tan(V1*4.0*atan(1.0)/180.0f); 
    TV2 = tan(V2);
    SV1 = sinf(V1*2*asin(1.0)/180.0f);
    SV2 = sinf(V2);
    CV1 = cosf(V1*acos(-1.0)/180.0f);
    CV2 = cosf(V2);

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

    R = 15.392f*B*(A*SV1*SV2+CV1*CV2*SA);

    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("SV1: %40.35f\n", SV1);
    printf("SV2: %40.35f\n", SV2);
    printf("CV1: %40.35f\n", CV1);
    printf("CV2: %40.35f\n", CV2);
    printf("SA : %40.35f\n", SA);
    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 = -tan(V1*4.0*atan(1.0)/180.0f);
    TV2 = tan(V2);
    SV1 = sinf(V1*2*asin(1.0)/180.0f);
    SV2 = sinf(V2);
    CV1 = cosf(V1*acos(-1.0)/180.0f);
    CV2 = cosf(V2);

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

    R = 15.392f*B*(A*SV1*SV2+CV1*CV2*SA);
 
    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("SV1: %40.35f\n", SV1);
    printf("SV2: %40.35f\n", SV2);
    printf("CV1: %40.35f\n", CV1);
    printf("CV2: %40.35f\n", CV2);
    printf("SA : %40.35f\n", SA);
    printf("A:   %40.35f\n", A);
    printf("B:   %40.35f\n", B);
    printf("R:   %40.35f\n", R);

    // case 3
    V1 = -7.23411989212036132812500000000000000f;
    V2 = -0.40112295746803283691406250000000000f;
    V3 = 2.0f;
    PI = 3.141592f;

    TV1 = -tan(V1*4.0*atan(1.0)/180.0f); 
    TV2 = tan(V2);
    SV1 = sinf(V1*2*asin(1.0)/180.0f);
    SV2 = sinf(V2);
    CV1 = cosf(V1*acos(-1.0)/180.0f);
    CV2 = cosf(V2);

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

    R = 15.392f*B*(A*SV1*SV2+CV1*CV2*SA);

    printf("Case 3:\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("SV1: %40.35f\n", SV1);
    printf("SV2: %40.35f\n", SV2);
    printf("CV1: %40.35f\n", CV1);
    printf("CV2: %40.35f\n", CV2);
    printf("SA : %40.35f\n", SA);
    printf("A:   %40.35f\n", A);
    printf("B:   %40.35f\n", B);
    printf("R:   %40.35f\n", R);

    return(0);
}