اصلاح کد نویسی اندرکنش خاک و سازه

سلام و عرض ادب. موضوع بحث امواج و پراکنش امواج در محیط نیمه بی نهایت هست.یک موج p به سیستم پی و سازه برخورد میکنه در مدل ما که یک تیر برشی هست جابجایی هایی به وجود میاره.که عبارت اند از فی دلتا وی.در کتاب تمام معادلات لازم جهت به دست آوردنش هست.حالا باید طبق منبع ما یعنی کتاب ماریا این معادلات بازکدنویسی بشن و به پلات هایی که اونا رسیدن ما هم برسیم. الان من کد رو تا اینجا رسوندم.یکی به ازای کا بتا آ برابر 1 به من خروجی میده ولی وقتی همون رو برای 4 زاویه و در بازه کا بتا آ صفر تا 6 ران می کنم.ران نمیشه. در وهله اول این مشکل فایل پلات حل بشه و در وهله دوم 450 خط کد بازبینی و با معادلات ماریا تطابق انجام بشه.به خصوص توابع هنکل و و بسل.تا در نهایت پلات هایی مشابه اون گرفته شود.مشکلی هم که هست اعداد بسیار بزرگ خروجی میده و دلیلش رو نمیدونم.من فایل کد به ازای عدد و فایل پلات رو میفرستم. در یک کلام به پلات های ماریا ما هم دست پیدا کنیم.پیام بدید تا توضیحات بیشتر و فایل کتاب رو بفرستم.