کدنویسی های بسیار طولانی و وقت گیری داره ولی شکر همه اطلاعات مورد نیاز رو بصورت کامل توضیح داده توی خود مقاله
متد دقیقا براساس مقاله پیاده سازی شده
تنها قسمتی که ممکن هست بخواهید ویرایش کنید تابع Lambda برای مقادیر ویژه است که من از یک تقریب سریع استفاده کردم که مقداری سرعت کد بره بالا بازم اگر برای معادله مشخصه رابطه دقیقتری برای محاسبه مقادیر ویژه برای تیر اویلر برنولی دارید برام قرار بدید کمتر از چند دقیقا براتون اعمال میکنم

همه توابعی که نوشتم برای استفاده کلی هستند یعنی میشه باهاش هر مثال دلخواهی رو حل کرد.
تا حد امکان سعی کردم دقیقا مشابه خود مقاله کاراکترهای توابع و نمادها رو نامگذاری و استفاده کنم.
اسکریپت اصلی برای اجرا Main_code هست و بقیه توابع هستند که توی برنامه فراخوانی میشن
توی مقاله از مقدار n=m=5 استفاده شده ولی هیچ محدودیتی برای افزایش یا کاهش این مقادیر و حتی توابع شرایط مرزی f,g وجود نداره و می تونید هر تابعی تعریف کنید. یه مقداری تابع بسط دهنده باعث کندی تحویل کار شد.
هر جا براتون سوال بود یا براتون مبهم بود برام قرار بدید.
مقاله رو خونده باشید همه توابع و روابط بین آنها رو کامل متوجه خواهید شد.