Решить дифференциальное уравнение пружинного маятника:
a = — k/m x
с помощью электронно-вычислительной машины.
Python-программа:
m = 1.5
x = 5.e-2
k = 2.
t = 0.
v = 0.
dt = 0.001
while t < 2.*3.14*(m/k)**0.5:
print(t, x, v)
t = t + dt
x = x + v*dt
v = v - k/m*x*dt
график положения грузика:
график скорости грузика:
