এটি গোলাকার পেন্ডুলামের একটি সিমুলেশন, বা দুটি মাত্রায় পেন্ডুলাম, যার ল্যাগ্রাঞ্জিয়ান দ্বারা দেওয়া হয়েছে
$
\mathcal{L} = \dfrac{1}{2}mL^2\dot{\theta}^2 + \dfrac{1}{2}mL^2\dot{\phi}^2\sin^2\theta +
$
$
+ mgL\cos\theta,
$
যেখানে $\theta$ এবং $\phi$ হল প্রথাগত গোলাকার স্থানাঙ্ক।
গতির সমীকরণগুলি চতুর্থ ক্রম রুঞ্জ-কুট্টা পদ্ধতি ব্যবহার করে সময় ধাপ $h$ সহ একত্রিত করা হয়েছে।
বড় প্রাথমিক কোণের জন্য (120º এর বেশি),
সমাধানের কনভারজেন্স নিশ্চিত করতে $h$ কমাতে হবে।
অভিকর্ষের ত্বরণ $g = 9.8$ m/s² হিসাবে সংজ্ঞায়িত করা হয়।
আরও পড়ুন: