我试图解决一个物理问题的微分方程:高尔夫球的坠落。我正在使用欧拉解析方法,我有这个代码:
def F(Z, t):
res = [Z[1], (-K/m)*Z[1]*sqrt((Z[1]**2) + (Z[1]**2)), Z[3], (-K/m)*Z[3]*sqrt((Z[3]**2) + (Z[3]**2))]
return res
def reso_z_euler(liste_t, F, CI, K):
n = len(liste_t)
Z = CI
pas = (liste_t[-1] - liste_t[0])/(n-1)
listeZ = [Z]
for i in range(1, n):
Z = F(Z, temps[i-1])
listeZ.append(Z)
return listeZ
listeZ = reso_z_euler(temps, F, Z0, K)
当我运行此代码时,终端抛出一个错误:OverflowError:(34,“结果太大”)我从未见过此错误,我不知道该怎么办。你的帮助会很好..谢谢
波斯汪
相关分类