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 = np.meshgrid(x, y) 第一个曲面方程 eq1 = x2 + z2 + y**2 - 4*y 第二个曲面方程 eq2 = x2 + y2 + z**2 - 4 设置绘图 fig = plt.figure() ax = fig.add_subplot(111, projection='3d') 绘制第一个曲面 ax.plot_surface(x, y,z,eq1, cmap='viridis') 绘制第二个曲面 ax.plot_surface(x, y,z,eq2,cmap='plasma') 设置坐标轴,标签 ax.set_xlabel('X') ax.set_ylabel('Y') ax.set_zlabel('Z') 显示图形 plt.show() 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 = np.meshgrid(x, y) 第一个曲面方程 eq1 = x2 + z2 + y**2 - 4*y 第二个曲面方程 eq2 = x2 + y2 + z**2 - 4 设置绘图 fig = plt.figure() ax = fig.add_subplot(111, projection='3d') 绘制第一个曲面 ax.plot_surface(x, y,z,eq1, cmap='viridis') 绘制第二个曲面 ax.plot_surface(x, y,z,eq2,cmap='plasma') 设置坐标轴,标签 ax.set_xlabel('X') ax.set_ylabel('Y') ax.set_zlabel('Z') 显示图形 plt.show() TypeError: Axes3D.plot_surface() takes 4 positional arguments but 5 were given

这个错误是由于你在绘制曲面时传递了多余的参数导致的。`Axes3D.plot_surface()` 方法只接受三个位置参数，分别是 `X`、`Y` 和 `Z`，而你在传递参数时多传递了一个 `eq1` 或 `eq2`。

要解决这个问题，你可以将 `eq1` 和 `eq2` 移除，并将其替换为对应的函数表达式。下面是修改后的代码：

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 = np.meshgrid(x, y)

# 第一个曲面方程
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, eq1, cmap='viridis')

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

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

# 显示图形
plt.show()

请注意，在计算 `eq1` 和 `eq2` 时，我将乘法操作符 `*` 用于对应的变量，以获得正确的计算结果。