[x(t_0) = x_0] ——————— import time from math import prod print("--- Khởi động Mô phỏng Hệ thống Động ---") # Các Điểm thu hút Y = [104, 101, 108, 108, 111, 32, 119, 111, 114, 108, 100] # Hàm Tiến hóa Phân tích x(t): Nội suy Lagrange điều khiển động lực học f = lambda t: sum(y * prod((t-k)/(j-k) for k in range(len(Y)) if k!=j) for j,y in enumerate(Y)) # Tiến hóa Thời gian: Bước qua thời gian rời rạc t, quan sát trạng thái x(t) acc = "" for t in range(len(Y)): c=chr(int(round(f(t)))); acc+=c; print(f"Thời gian t={t:2} | x(t): '{c}' | Đường đi: {acc}"); time.sleep(0.2)