def METHOD_NAME(delta_rho, a, z, stationSpacing):
    respEW = 0
    gravity_change = []
    xmax = 10.0
    npts = int(1 / stationSpacing)
    x = np.linspace(-xmax, xmax, num=npts)
    y = x.copy()
    X, Y = np.meshgrid(x, y)
    for each in range(npts):
        gravity_change.append(calculategravity(delta_rho, a, z, x[each]))
    return respEW, gravity_change, X, Y