如何在Scipy中使用柱坐标进行积分？

Question

我知道怎么用Scipy来进行函数的积分。我是这样做的：

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

这样我就用辛普森法则对函数y(x)进行了积分。但是，如果我想用圆柱坐标来积分这个函数呢？也就是说，我想计算的不是普通的积分 \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 是你在计算时需要加到积分函数里的一个必要因素（就像你自己提到的那样）。