def METHOD_NAME(delta_rho, a, z, x) -> float:
    partial = (4000 / 3 * np.pi * G) * delta_rho * a * a * a * z
    g = partial / float(np.power(x * x + z * z, 1.5))
    response = g * pow(10, 5)
    return response