quadメソッドのpythonを使った統合

Question

erfメソッドを使用せずに、scipyを使ってエラー関数を再現しようとしています。

これは私のコードです：

#!/usr/bin/env python 
import numpy as np 
import math 
from scipy.integrate import quad 

def error_func(y): 
    return 2/math.sqrt(np.pi)*quad(np.exp(-y**2), 1, np.inf, args=(y))[0] 

g = [error_func(x) for x in np.arange(-1,1,0.2)] 

print g 

このコードは、次のエラーメッセージが返されます：

File "./test.py", line 9, in <module> 
    g = [error_func(x) for x in np.arange(-1,1,0.2)] 
    File "./test.py", line 7, in error_func 
    return 2/math.sqrt(np.pi)*quad(np.exp(-y**2), 1, np.inf, args=(y))[0] 
    File "/usr/local/lib/python2.7/site-packages/scipy/integrate/quadpack.py", line 316, in quad 
    points) 
    File "/usr/local/lib/python2.7/site-packages/scipy/integrate/quadpack.py", line 383, in _quad 
    return _quadpack._qagie(func,bound,infbounds,args,full_output,epsabs,epsrel,limit) 
quadpack.error: quad: first argument is not callable 

を私が正しく理解していれば、クワッドの最初の引数は関数でなければなりません。私はそれを正しく伝えます。 私のコードで何が問題になっていますか？

Answer 1

あなたの最初の問題は、

np.exp(-y**2) 

は関数ではないということです。それは数字です。具体的には、e ^（ - y^2）の値です。しかし、あなたが他の問題を抱えている

lambda y: np.exp(-y**2) 

：あなたはnp.exp(-y**2)にyをマップする関数を定義したい場合は、最も簡単な方法は、lambda構文である

(y)

は args=(y)でタプルではありません。グループを小括弧で囲んだのはちょうど yです。 1要素タプルは (y,)です。

また、args引数をquadに渡すべきではありません。なぜなら、あなたの関数は、統合している引数以外の追加引数を必要としないからです。

最後に、統合の境界が間違っています。 0からエラー関数の引数に統合する必要があります。