0

计算此积分的任何帮助，F函数都使用f涉及第一次积分的函数定义，最后是积分F。

from scipy.integrate import quad

f = lambda x,a : a**2*x

def F(s,a):
  return quad(f,0,s,args=(a,))
quad(F,0,5,args=(4,))

得到错误：

      2 def F(s,a):
      3   return quad(f,0,s,args=(a,))
----> 4 quad(F,0,5,args=(4,))
      5 
    446     if points is None:
    447         if infbounds == 0:
--> 448             return _quadpack._qagse(func,a,b,args,full_output,epsabs,epsrel,limit)
    449         else:
    450             return _quadpack._qagie(func,bound,infbounds,args,full_output,epsabs,epsrel,limit)

TypeError: must be real number, not tuple
4

2 回答 2

1

看看返回值scipy.integrate.quad：

回报：

y: floatfunc 从 a 到 b 的积分。

abserr:float结果中绝对误差的估计。

...

所以有多个返回值（一个元组），这就是你收到TypeError: must be real number, not tuple消息的原因。

我想，你只是对积分值感兴趣，quad(...)[0]所以这就是你F应该返回的：

from scipy.integrate import quad

f = lambda x, a: a**2 * x
F = lambda x, a: quad(f, 0, x, args=(a,))[0]

I = quad(F, 0, 5, args=(4,))
print(I)

哪个打印：

(333.33333333333337, 3.700743415417189e-12)
于 2019-07-25T16:47:11.253 回答
0

另一种看待问题的方法是意识到您正在将函数积分到由点、和a**2 * y跨越的三角形上。然后，您可以使用quadpy（我的项目之一）中的三角形求积来计算值：[0, 0][5, 0][5, 5]

import quadpy

a = 4.0


def g(x):
    return a ** 2 * x[1]


scheme = quadpy.triangle.dunavant_05()
val = scheme.integrate(g, [[0.0, 0.0], [5.0, 0.0], [5.0, 5.0]])
print(val)

与嵌套四边形方法相比，这应该需要更少的函数评估。

于 2019-10-30T20:35:38.853 回答