Период отдыха для простого маятника с использованием Рунге Кутты – python ⇐ Python

[Anonymous]

Я пытаюсь смоделировать простой маятник, используя методы Рунге-Кутты в Python, и период должен составлять 2 секунды. Мало того, мои значения угловой скорости и углового смещения также неверны. Я не уверен, в чем моя ошибка в моем коде.
Это код, который я пробовал: v1 — это шаг скорости, а t1 — шаг углового смещения.

Код: Выделить всё

g = 9.81
pi = math.pi
l = g/pi

t0 = 0
t = t0
t_list = [t0]

v0 = 0
v = v0
v_list = [v0]

theta0 = 0.1
theta = theta0
theta_list = [theta0]

dt = 0.02
tfinal = 20

def theta2(t, v, theta):
return -(g/l)*math.sin(theta)

def theta1(t, v, theta):
return v

while t 

Подробнее здесь: [url]https://stackoverflow.com/questions/79008357/period-off-for-simple-pendulum-using-runge-kutta-python[/url]