0

我试图构建一个图形，显示非线性周期（相对于线性周期，T0）作为初始位置的函数，θ0 和两种方法之间的差异。

def f(theta, theta_0, L, g):
    return 4 * np.sqrt(L/(2*g)) * 1 / np.sqrt(np.cos(theta) - np.cos(theta_0))

theta_0 = 0.1
L = 1 
g = 1
quad(f, 0, theta_0, args=(theta_0, L, g))

(6.287114549357374, 5.232019262280119e-10)

和特殊功能：

def f(theta, theta_0, L, g):
return 4 * np.sqrt(L/g) * ellipkinc(np.pi/2, np.sin(theta_0/2) )

theta_0 = 0.1
L = 1 
g = 1
quad(f, 0, theta_0, args=(theta_0, L, g))

(0.6363979153194195, 7.065436184110701e-15)

如何生成它们之间的误差差异图？在此处输入图像描述

4

0 回答 0