还是有问题TypeError: Axes3D.plot_surface() takes 4 positional arguments but 5 were given

非常抱歉，我之前的回答有误。这是因为我引用了错误的函数名称。在绘制曲面时，应该使用`plot_surface`函数而不是`contour3D`函数。

请尝试以下修改后的代码：

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# 创建网格点
x = np.linspace(-4, 4, 100)
y = np.linspace(1, 3, 100)
z = np.linspace(-4, 4, 100)
x, y, z = np.meshgrid(x, y, z)

# 第一个曲面方程
eq1 = x**2 + z**2 + y**2 - 4*y

# 第二个曲面方程
eq2 = x**2 + y**2 + z**2 - 4

# 设置绘图
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

# 绘制第一个曲面
ax.plot_surface(x, y, z, cmap='viridis', alpha=0.8)

# 绘制第二个曲面
ax.plot_surface(x, y, z, cmap='plasma', alpha=0.8)

# 设置坐标轴标签
ax.set_xlabel('X')
ax.set_ylabel('Y')
ax.set_zlabel('Z')

# 显示图形
plt.show()

这样修改后，你应该能够成功绘制出三维曲面图。非常抱歉之前的回答给你带来了困扰。