Fortran Subroutine For Simpson's Method of Numerical Integration

Abstract

The study discusses Simpson's method for numerical integration using a Fortran subroutine. The routine aims to compute the numerical integral of the function \( f(x) \) from point \( a \) to point \( b \) using Simpson's one-third rule. The function is evaluated at \( n+1 \) points, dividing the interval into smaller segments. Simpson's rule states that the integral can be computed using a combination of function values at various points. The routine starts with an initial interval \( n \) and repeats the computation by doubling the interval to achieve a more accurate estimate of the integral. The truncation error is compared with a user-defined tolerance, and if the error exceeds this value, the process is repeated by doubling the intervals for improved results.

يتناول البحث طريقة سيمبسون للتكامل العددي باستخدام روتين فرتران. يهدف الروتين إلى حساب التكامل العددي للدالة \( f(x) \) من النقطة \( a \) إلى النقطة \( b \) باستخدام قاعدة سيمبسون للثلث. يتم حساب الدالة عند \( n+1 \) نقطة، حيث يتم تقسيم الفترة إلى فترات صغيرة. تنص قاعدة سيمبسون على أن التكامل يمكن حسابه من خلال تركيبة من قيم الدالة عند النقاط المختلفة. يبدأ الروتين بتحديد الفترة الأولية \( n \) ثم يقوم بتكرار الحساب عن طريق مضاعفة الفترة للحصول على تقدير أكثر دقة للتكامل. يتم مقارنة خطأ التقريب بالقيمة المعطاة من قبل المستخدم، وإذا تجاوز الخطأ هذه القيمة، يتم تكرار العملية مع مضاعفة الفترة للحصول على نتائج أفضل.