关于scipy.integrate.quad嵌套积分函数的实现优化问询

首先，你的代码存在核心错误：scipy.integrate.quad要求第一个参数是可调用的被积函数，但你传入的myreturn1、myreturn2是基于输入x计算出的固定数值，这会直接导致quad抛出类型错误。

结合你需要嵌套积分、返回两个结果的需求，以下是优化实现：

基础优化版本（修正错误+减少冗余） #

import numpy as np
from scipy.integrate import quad

def myfunction(x):
    # 计算基础表达式
    function_0 = 1.25 * x ** 2
    function_1 = function_0 * 9.81 + 2 * x
    
    # 定义内层被积函数（这里假设内层积分变量为t，可根据实际需求替换）
    def integrand(t):
        # 这里填写实际的内层被积逻辑，示例中假设function_1是被积项的系数
        return function_1
    
    # 仅计算一次基础积分，避免重复调用quad浪费资源
    base_integral, _ = quad(integrand, 0, np.inf)
    
    # 基于基础积分生成两个返回值
    val1 = base_integral * np.pi
    val2 = base_integral * np.pi ** 2
    
    return val1, val2

# 调用示例：传入具体x值
x = 1.5
val1, val2 = myfunction(x)

优化点说明 #

• 修正quad参数错误：传入可调用的integrand函数，符合quad的参数要求。
• 减少重复计算：只执行一次内层积分，再通过系数生成两个结果，提升运行效率。
• 结构清晰：将基础计算、被积函数定义、积分逻辑拆分，可读性更强。

若需外层积分（二重嵌套场景） #

如果你的"两层嵌套积分"是指外层对x积分，内层对其他变量积分，可以进一步封装：

import numpy as np
from scipy.integrate import quad

def calculate_inner_integrals(x):
    function_0 = 1.25 * x ** 2
    function_1 = function_0 * 9.81 + 2 * x
    
    base_integral, _ = quad(lambda t: function_1, 0, np.inf)
    return base_integral * np.pi, base_integral * np.pi ** 2

# 外层对x从0到10积分，得到两个总积分结果
total_val1, _ = quad(lambda x: calculate_inner_integrals(x)[0], 0, 10)
total_val2, _ = quad(lambda x: calculate_inner_integrals(x)[1], 0, 10)

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