Julia双重积分问题求助：同类函数一个积分正常另一个出现类型转换错误

核心问题分析 #

你遇到的MethodError本质是**hcubature要求被积分函数返回单个Float64标量，但你的cop2019函数在接收标量输入时，返回的是长度为1的Vector{Float64}**，导致类型转换失败。

对比能正常运行的joedensity：当传入标量u、v时，它直接返回单个浮点数，完全符合hcubature的要求；而cop2019的逻辑是为批量向量输入设计的——即使你传入标量x[1]、x[2]，函数内部会生成2×1的矩阵y，通过eachcol(y)迭代后，map会输出长度为1的向量，这和hcubature期望的标量返回值冲突。

修复方案（两种可选） #

方案1：修改cop2019函数，适配标量/向量双输入

调整函数逻辑，自动判断输入类型，标量输入返回标量，向量输入返回向量：

function cop2019(u,v,delta,thetaW)
    # 检测输入是否为标量
    is_scalar_input = isa(u, Number) && isa(v, Number)
    # 根据输入类型初始化y矩阵
    y_len = is_scalar_input ? 1 : length(u)
    y = zeros(2, y_len)
    y[1,:] .= quanty.(u,delta)
    y[2,:] .= quanty.(v,delta)
    
    # 计算核心结果
    joint_vals = map(x -> jointdist(x,delta,thetaW), eachcol(y))
    marg_prod = marg.(y[1,:],delta) .* marg.(y[2,:],delta)
    res = joint_vals ./ marg_prod
    
    # 标量输入返回单个值，否则返回向量
    return is_scalar_input ? first(res) : res
end

修改后直接调用原积分代码即可：

hcubature(x -> cop2019(x[1],x[2],delta,thetaW),[0,0],[1,1],abstol=0)

方案2：在积分调用时直接提取标量

如果不想修改cop2019，可以在传给hcubature的匿名函数里，直接取返回向量的第一个元素：

hcubature(x -> first(cop2019(x[1],x[2],delta,thetaW)),[0,0],[1,1],abstol=0)

验证结果 #

修复后运行积分代码，你会得到接近1的结果，比如类似：

(0.9999999999712456, 1.087234519876e-8)

符合你预期的积分结果接近1的要求。

内容的提问来源于stack exchange，提问作者maria96