在Python中将子绘图合并到单个绘图中

所以我有一个当前函数，它给出了一个图，图中有一个圆和一个位于给定圆中心的点：

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)

这里有两个问题：

每次运行该函数时，我都会得到一个单独的绘图。我想运行该函数几次，并将结果绘制在单个图形上。我该怎么做

运行此函数3次后，我将有3个圆在给定点相交。如何确定交点的x和y值

---

每次调用plt.subplot都会生成一个新的数字。建议您在函数外部调用它，并将ax作为parameter@JodyKlymak因此，将plt.subplot移到外面肯定给了我想要的，但是，我不知道如何将ax作为参数传递。你能解释一下那个部分吗？

台站地震（（10，10），3.92，ax）

？显然，您还必须更改函数