Đây là mô phỏng của con lắc hình cầu, hay con lắc hai chiều, có Lagrangian được cho bởi
$$
\mathcal{L} = \dfrac{1}{2}mL^2\dot{\theta}^2 + \dfrac{1}{2}mL^2\dot{\phi}^2\sin^2\ theta +
$$
$$
+ mgL\cos\theta,
$$
trong đó $\theta$ và $\phi$ là tọa độ hình cầu truyền thống.
Các phương trình chuyển động được tích hợp bằng phương pháp Runge-Kutta bậc 4 với bước thời gian $h$.
Đối với các góc ban đầu lớn (lớn hơn 120°),
cần phải giảm $h$ để đảm bảo sự hội tụ của nghiệm.
Gia tốc trọng trường được định nghĩa là $g = 9,8$ m/s².
Đọc thêm: