


import math

lat1 = 35.660386
long1 = -78.975496
lat2 = 35.699844
long2 =-78.942238

def calc_bearing(lat1, long1, lat2, long2):
    # Convert latitude and longitude to radians
    lat1 = math.radians(lat1)
    long1 = math.radians(long1)
    lat2 = math.radians(lat2)
    long2 = math.radians(long2)

    ### Calculate the bearing
    bearing = math.atan2(
        math.sin(long2 - long1) * math.cos(lat2),
        math.cos(lat1) * math.sin(lat2) - 
        math.sin(lat1) * math.cos(lat2) * 
        math.cos(long2 - long1)
    )###
    # Convert the bearing to degrees
    bearing = math.degrees(bearing)
    
print (lat1)
print (long1)

print
print(lat2)
print(long2)
print
print (bearing)
