کد متلب حل معادله دیفرانسیل معمولی با روش عددی رانگ کوتای مرتبه چهار با توضیحات خط به خط در ورد
برنامه ابتدا در خطوط 7 و 8 به ترتیب، مقادیر x0 و xn یا همان نقطهی ابتدا و انتهای بازه را از کاربر دریافت میکند. سپس در خط 9 مقدار h که طول گام است را از کاربر گرفته و در خط 10 مقدار تابع در نقطه شروع یعنی y(x0) را از کاربر دریافت میکند.
سپس در خط 11 تابع را از کاربر سوال میکند و آن را در متغیر fun ذخیره میکند.
در خط 13 برنامه رشتهای که از x0 شروع و تا xn ادامه دارد با طول گام h میسازد. در خط 14 این رشته را در ستون اول ماتریس y که تمام جواب های نهایی در آن قرار دارد، میریزد. در خط 15 نیز مقدار اولیه را در ماتریس نهایی جواب قرار میدهد.
در بلوک for مربوط به خطوط 16 تا 21 برنامه مقدار تابع را به دست میآورد به این شکل که ابتدا مقادیر k1, k2, k3, k4 را حساب کرده و سپس با افزودن مقدار (k1+2*k2+2*k3+k4)/6، به نقطهی قبلی مقدار تابع در نقطه جدید را محاسبه میکند.
در بلوک مربوط به خطوط 26 تا 30 ضمن محاسبهی مشتق تابع در تمام نقاط بازه مقادیر x,y,y’ را که به ترتیب در ستون اول تا سوم ماتریس y قرار دارند را نمایش میدهد. در خط 29 هم با 0.1 ثانیه تاخیر به سراغ نقطه بعدی میرود.
در خطوط 31 تا 36 هم مقادیر x, y , y’ رسم میشوند.