SymPy定积分中的ZeroDivisionError；可变积分界
SymPy定积分中的ZeroDivisionError；可变积分界
提问于 2020-04-17 06:01:30
我使用Julia并调用SymPy对其两个变量a和n进行二变量正态分布函数的积分。我希望在集成边界中包含第三个变量p，但当我尝试此变量时，SymPy会抛出ZeroDivisionError。
运行以下命令，
using SymPy
σ = 0.2
ρ = -0.9
f(a,n) = (1/(2*π*(σ^2)*sqrt(1-ρ^2)))*exp(-(1/(2*(1-ρ^2)))*((a^2)/σ^2 + (n^2)/σ^2 -((2*ρ*a*n)/σ^2)))
@vars a n p
integrate(a*f(a,n), (n, -oo, p*a), (a, -oo, oo))收益率
ERROR: PyError ($(Expr(:escape, :(ccall(#= C:\Users\Anshu\.julia\packages\PyCall\zqDXB\src\pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'ZeroDivisionError'>奇怪的是，以下尝试运行没有任何问题：
只是对f进行积分，而不是a*f：
integrate(f(a,n), (n, -oo, p*a), (a, -oo, oo))将p替换为常量1：
integrate(a*f(a,n), (n, -oo, a), (a, -oo, oo))回答 1
Stack Overflow用户
发布于 2020-04-17 08:00:50
在定义sigma和rho之后，我得到：
julia> integrate(a*f(a,n), (n, -oo, p*a), (a, -oo, oo)) |> string
"0.159154943091895*Integral(a*exp(-a^2/(-2*ρ^2*σ^2 + 2*σ^2))*Integral(exp(-n^2/(-2*ρ^2*σ^2 + 2*σ^2))*exp(2*a*n*ρ/(-2*ρ^2*σ^2 + 2*σ^2)), (n, -oo, a*p)), (a, -oo, oo))/(σ^2*sqrt(-(ρ - 1)*(ρ + 1)))"也许这是你的症状版本，我有
julia> sympy.__version__
"1.5.1"页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/61260761
复制相关文章
相似问题
腾讯云开发者
