この質問は、 GSL 統合の動作がおかしいという以前の質問のフォローアップのようなものです
スケーリングx -> (1-t)/tが無限積分法で望ましくない結果を引き起こすため、gsl_integration_qagi私は現在、有限サポート上の積分を使用することに頼っています。私が今していることはこれです：
私は（離散）一連Sの実数を持っています。これを（それぞれ）特定の指数関数で畳み込みます
exp(-t/T) ... T = decay constant
(min(S) - 10*T, max(S) + 10*T)関数からの「重要な」貢献のほとんどをカバーできるように、統合のサポートを選択します。
gsl_integration_qagi` gsl_integration_qag(無限サポート) は数ミリ秒しかかかりませんが、間違った結果を生成します。誰も理由を知っていますか？
畳み込みが指数ではなくガウスである場合、同じgsl_integration_qagことがうまく機能します。
前もって感謝します、ニキル