def bilateration(latA, lonA, distA, latB, lonB, distB):
    """
    @param latA: latitude of point 1
    @type latA: C{float}
    @param lonA: longitude of point 1
    @type lonA: C{float}
    @param distA: distance from point1 to the third point in miles
    @type distA: C{int}
    @param latB: latitude of point 1
    @type latB: C{float}
    @param lonB: longitude of point 1
    @type lonB: C{float}
    @param distB: distance from point2 to the third point in miles
    @type distB: C{int}
    """
    #Converting miles to nautical miles
    distA *= 0.868976
    distB *= 0.868976

    xA = math.cos(lonA)*math.cos(latA)
    yA = math.sin(lonA)*math.cos(latA)
    zA = math.sin(latA)

    xB = math.cos(lonB)*math.cos(latB)
    yB = math.sin(lonB)*math.cos(latB)
    zB = math.sin(latB)

    x1 = numpy.array([xA, yA, zA])
    x2 = numpy.array([xB, yB, zB])

    #Converting from nautical miles to radians
    r1 = distA / 60
    r2 = distB / 60

    q = numpy.dot(x2, x1)
    a = (math.cos(r1) - math.cos(r2)*q) / (1 - pow(q,2))
    b = (math.cos(r2) - math.cos(r1)*q) / (1 - pow(q,2))

    n = numpy.cross(x1,x2)

bilateration(37.673442, -90.234036, 123.709, 36.109997, -90.953669, 166.863)