如何将四阶级数展开代码修改为支持n阶通用形式？

现有一段四阶级数展开的Python代码：

d=[]  
for i in range(3,51):  
     d.append(C[0][0] + C[0][1]/(i-do)**2 + C[0][2]/(i-do)**4 + C[0][3]/(i-do)**6)

需要修改为支持n阶通用形式的代码：当为五阶时自动添加C[0][4]/(i-do)**8这类项，最终实现累加至C[0][r]/(i-do)**2r的形式（r从1到目标阶数对应的索引）。

方法1：循环累加（直观易懂） #

通过嵌套循环逐项累加，逻辑清晰，适合快速调试和理解：

d = []
target_r = 4  # 对应五阶，原四阶设为3即可
do = 0  # 替换为你的实际do值
C = [[1, 2, 3, 4, 5]]  # 替换为你的实际系数矩阵

for i in range(3, 51):
    total = C[0][0]  # 初始常数项
    denom = i - do
    # 从r=1到target_r，逐项累加
    for r in range(1, target_r + 1):
        total += C[0][r] / (denom ** (2 * r))
    d.append(total)

方法2：sum()结合生成器表达式（简洁高效） #

利用Python的生成器表达式配合sum()函数，一行完成项的累加，代码更紧凑：

d = []
target_r = 4  # 对应五阶，原四阶设为3
do = 0  # 替换为实际值
C = [[1, 2, 3, 4, 5]]  # 替换为实际系数矩阵

for i in range(3, 51):
    denom = i - do
    # 生成器遍历r=1到target_r，sum直接计算累加和
    total = C[0][0] + sum(C[0][r] / (denom ** (2 * r)) for r in range(1, target_r + 1))
    d.append(total)

关键说明 #

• 只需修改target_r的值即可切换阶数：比如原四阶代码对应target_r=3，五阶对应target_r=4，以此类推。
• 如果denom可能为0，建议添加判断避免除以0的错误：比如在计算前加入if denom == 0: 处理逻辑。

内容的提问来源于stack exchange，提问作者Shafiq Ur Rehman