DEFINT I-O
'               program to calculate the sides of the pyramid
'        INPUT "HEIGHT OF PYRAMID, LATITUDE, ALPHA =  ", H, ALAT, ALPHA
'               THREE POINTS ON THE SIDE TRIANGLE :
'               A - SOUTH POINT
'               B - APEX (0,0,H)
'               D - BASE POINT OF NORTH FACE
'
'               EACH POINT HAS X,Y,Z
        PI = 3.14159
'               height
        h = 10
'               latitude
        ALAT = 34
        CLS
        PRINT "    ALPHA    A - B    A - D    B - D   2(C - D)"
        FOR ALPHA = 30 TO 60 STEP 2
        ALATRADIANS = ALAT * PI / 180
        ALPHARADIANS = ALPHA * PI / 180
        TANLAT = TAN(ALATRADIANS)
        TANALPHA = TAN(ALPHARADIANS)
        COSLAT = COS(ALATRADIANS)
        AX = -h / TANLAT
        AY = 0
        AZ = 0
        BX = 0
        BY = 0
        BZ = h
        DX = h * TANLAT
        DY = h / COSLAT / TANALPHA
        DZ = 0
'               CALCULATE SIDES   
        AB = SQR((AX - BX) ^ 2 + (AY - BY) ^ 2 + (AZ - BZ) ^ 2)
        AD = SQR((AX - DX) ^ 2 + (AY - DY) ^ 2 + (AZ - DZ) ^ 2)
        BD = SQR((BX - DX) ^ 2 + (BY - DY) ^ 2 + (BZ - DZ) ^ 2)
        D2 = 2 * DY
        PRINT USING "######.##"; ALPHA; AB; AD; BD; D2
        NEXT ALPHA
        END