Объединение субплотов в один сюжет в Python

Question

Итак, у меня есть текущая функция, которая дает мне график с кругом и точкой в ​​центре данного круга:

def epicentral_distance(station, differential_time_variable, Vp = 5.5, Vs = 3.2):
    ri = differential_time_variable*(1/Vs - 1/Vp)**(-1)
    print( np.round(ri,3) )
    return np.round(ri,3)


def station_earthquake_GraphFunction(station,differential_time_variable, Vp = 5.5, Vs = 3.2):
    radius = epicentral_distance(station, differential_time_variable, Vp = 5.5, Vs = 3.2)
    fig, ax = plt.subplots(figsize = (5,5))

    graph_station = plt.plot(station[0],station[1],'^r')
    graph_circle = plt.Circle(station, radius, color = 'r',fill = False)

    ax.add_artist(graph_circle)
    plt.xticks(np.arange(-100,100,20))
    plt.yticks(np.arange(-100,100,20))
    return graph_station, graph_circle

Я пытаюсь запустить функцию со следующими значениями:

station_earthquake_GraphFunction((10,10), 3.922)
station_earthquake_GraphFunction((15,45),3.272)
station_earthquake_GraphFunction((60,30),3.874)

Два вопроса здесь:

1. каждый раз, когда я запускаю функцию, я получаю отдельный график. Я хочу запустить функцию несколько раз и получить результаты на одном графике. Как мне это сделать?

2. После трехкратного выполнения этой функции у меня будет 3 круга, пересекающихся в данной точке. Как я могу определить значения x и y пересечения?