python - 将最后一个迭代值分配给python中的变量

Question

我有迭代的最后一组值:y1[i],y2[i],y3[i](通过耦合 ODE 积分获得)现在我需要将这些值分配给另一个变量:Y1、Y2、Y3，然后在函数 f(y,t) 中使用这些变量(必须再次积分的 ODE)。这是我的代码的一部分:

#constants used
H = 2.27e-18
l = 1.5
G = 6.637*(10**(-11))
k = (8*3.14*G)**0.5
om_de = 0.75
omega_matter = 1 - om_de
w0 = -0.85
rho_c = 3*(H**2)/(k**2)
c = ((om_de*(1 + w0)*rho_c)/(H**2))**0.5
v0 = (om_de*(1 - w0)*rho_c)/(2*np.exp(-l*k))

#for iteration
p = 0.0
q = 1.0
n = 10000.0

def H(Y1,Y2,Y3):
    return -3*Y3*(H**2)*(omega_matter/(Y1**3) + (H**2)*(Y3**2)/(2.0*rho_c) + (v0*np.exp(-l*Y2*k))/(rho_c))**0.5 + (l*k*v0*np.exp(-l*Y2*k))/(H**2)

dH = (q-p)/n
Y1 = [y1[j]]
Y2 = [y2[j]]
Y3 = [y3[j]]
But I am not able to do any further operations on them.

TypeError: unsupported operand type(s) for ** or pow(): 'function' and 'int'

不断出现。完整代码的链接在这里: https://pastebin.com/mdQE29N9 (可能有助于理解我想要实现的目标并且代码太长而无法在此处发布)

我认为问题在于我分配值的方式，但我不太确定。任何有关导致此问题的原因以及我应该如何解决此问题的帮助都会很棒。

Answer 1

第 57 行有一个错误:

def F(Y1,Y2,Y3):
    return Y1*(omega_matter/(Y1**3) + (H**2)*(Y3**2)/(2.0*rho_c) + (v0*np.exp(-l*Y2*k))/(rho_c))**(1.0/2.0)

在表达式 H**2 中，其中 H 在第 62 行被定义为一个函数:

def H(Y1,Y2,Y3):
    ...

我猜你真的想说:H(Y1,Y2,Y3)**2。

但这还不是全部。第 52 行:

y1 == y1.append(y1[j] + (dh/6.0)*(k1 + 2.0*k2 + 2.0*k3 + k4))

无需测试 y1 和 None 之间的等价性(这就是 append() 返回的内容)。我不确定您在这里尝试实现的目标，也许是这个？

y1.append(y1[-1] + (dh/6.0)*(k1 + 2.0*k2 + 2.0*k3 + k4))

此外，要访问 y1 的最后一个元素(第 66 行)，您可以使用负数索引:

Y1 = [y1[-1]]   # instead of: Y1 = [y1[j]]