Определить величину ускорения тела в каждый момент времени

import numpy as np
import math
import matplotlib.pyplot as plt
from scipy.integrate import simps
 
# Задаем скорости по осям x и y
vx = np.array([21.65, 15.27, 11.88, 9.78, 8.35 ,7.31 ,6.51, 5.87, 5.34, 4.86, 4.44, 4.04, 3.68, 3.33, 3.01, 2.71, 2.43, 2.17, 1.94, 1.72, 1.53])
vy = np.array([12.50, 8.01, 5.37, 3.53, 2.11, 0.93, -0.10, -1.02, -1.86, -2.63, -3.34, -3.98, -4.55, -5.05, -5.49, -5.87, -6.19, -6.46, -6.68, -6.87, -7.02  ])
t = np.array([0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0])
 
 
# Находим ускорения по осям x и y
ax = np.diff(vx) / t
ay = np.diff(vy) / t
 
# Находим полное ускорение
a = np.sqrt(ax**2 + ay**2)
 
print("Ускорение по x:")
print(ax)
print("Ускорение по y:")
print(ay)
print("Полное ускорение:")
print(a)
 
def f(x):
    return x**2
    
def trapz(f, a, b, n):
    h = (b - a) / n
    result = 0.5 * (f(a) + f(b))
    for i in range(1, n):
        result += f(a + i * h)
    result *= h
    return result
 
def simps(f, a, b, n):
    h = (b - a) / n
    result = f(a) + f(b)
    for i in range(1, n, 2):
        result += 4 * f(a + i * h)
    for i in range(2, n-1, 2):
        result += 2 * f(a + i * h)
    result *= h / 3
    return result
 
# График изменения ускорения тела во времени
time = np.arange(len(ax))
plt.plot(1,2,2)
plt.plot(t, ax, label='Ускорение по x', color='blue')
plt.plot(t, ay, label='Ускорение по y', color='red')
plt.title('Изменение ускорения тела во времени')
plt.xlabel('Время')
plt.ylabel('Ускорение')
plt.legend()
plt.grid(True)