如何使用 Scipy 整合柱坐标？

Question

我知道如何将函数与 Scipy 集成。我是这样做的：

    from scipy import *
    from scipy import integrate
    integral = integrate.simps(y,x)

通过这种方式，我将函数 y(x) 与 Simpsons 规则积分，但如果我想将此函数积分到柱坐标中怎么办？ 我的意思是，而不是积分 \int y(x)dx 我想解决积分 \int y(x)*2*pi*x*dx。

我试过了

    integral = integrate.simps(y,x**2.*pi)

但它似乎没有返回正确的结果。

我正在集成数据点，而不是函数！

Answer 1

你应该试试

integral = integrate.simps(y*x*2*pi,x)

其中y 应该是位于x 位置的函数值数组。请注意，2*pi*x 是添加到被积函数的必要因素（正如您自己指出的那样）。