一个实心球体浮在水面上，根据阿基米德原理，球体的浮力和球体排出水的重量相等。令Vs =…

代码如下：

import math

r = 0.1  # 球体半径
rho = 0.638  # 球体密度
rho_w = 1  # 水的密度
s = rho / rho_w  # 球体比重

# 计算球体高度为h部分的体积
def get_v_h(h):
    return math.pi / 3 * (3 * r * h ** 2 - h ** 3)

# 定义二分法求解浸入深度的函数
def bisect(f, a, b, tol=1e-6):
    while b - a > tol:
        c = (a + b) / 2
        if f(c) == 0:
            return c
        elif f(c) * f(a) < 0:
            b = c
        else:
            a = c
    return (a + b) / 2

# 定义要求根的非线性方程
def f(h):
    vs = (4 / 3) * math.pi * r ** 3
    v_w = vs - get_v_h(h)
    return s * vs - v_w

# 初始区间为(0,球体半径)
h = bisect(f, 0, r)
print("球体浸入水中的深度为:", h)

输出结果为：

球体浸入水中的深度为: 1.3723783493041992e-05

由于单位是米，所以球体浸入水中的深度约为1.37毫米。