看起来你的尝试很接近。
诀窍是使用f(x)作为半径创建磁盘。所以每个磁盘的面积是pi * f(x)**2 ...
我相信这就足够了。

import numpy as np
from scipy import integrate
import matplotlib.pyplot as plt

def f(x):
    ''' some function (to be roatated around the x-axis)'''
    return x**2

# make some x and y values
x = np.linspace(0, 3)
y = f(x) 

def g(x):
    ''' volume of the curve by integrating the circle slices '''
    radius = f(x)
    return np.pi * radius**2

volume, error = integrate.quad(g, 0, 3)

print(volume)

字符串
在该（示例）情况下的结果是152.68…
注意，要围绕y-axis旋转，您需要切换x和y。