Trato de resolver este sistema de EDOs y me aparece el error cannot unpack non-iterable float object y no sé como resolverlo
Esta es la función a integrar:
def sist(z, t, k, m):
phi1, phi2, phi3, w1, w2, w3 = z
w = k/m
phi1_dot = w1
phi2_dot = w2
phi2_dot = w3
w1_dot = -w*(phi1 - phi2) - w*(phi1 - phi3)
w2_dot = -w*(phi2 - phi3) - w*(phi2 - phi1)
w3_dot = -w*(phi3 - phi1) - w*(phi3 - phi2)
return w1, w2, w3, w1_dot, w2_dot, w3_dot
Este es el método:
# MÉTODO DE RK4:
t0 = 0.
tf = 50.
k = 1.
m = 1.
z = solve_ivp(sist, t_span=(t0,tf), y0=[1., 1., 1., np.pi/40, np.pi/20, np.pi/10], method='RK45',
args=(k, m), vectorized=False)
sist()
comoz
es unfloat
y cuando intentas hacerphi1, phi2, phi3, w1, w2, w3 = z
no puede desempaquetarlo.phi1 = phi2... = z
. Esta respuesta te puede ayudar a entender el unpacking