محاكاة مونت كارلوا في برنامج R

في هذا المثال ، سنستخدم محاكاة مونت كارلو لتقدير قيمة تكامل الدالة f (x) = x ^ 2 خلال الفترة [0 ، 1].في هذا المثال ، نستخدم الدالة (runif)  لإنشاء عدد n من النقاط عشوائيًا بين 0 و 1 ، وهو ما يمثل إحداثيات x للنقاط ضمن فترة التكامل. نقوم أيضًا بإنشاء عدد n من النقاط إحداثيات y عشوائية بين 0 والقيمة القصوى للدالة خلال فترة التكامل. 

نحسب بعد ذلك عدد النقاط التي تقع تحت الدالة f (x) = x ^ 2 ونقرب المساحة الواقعة أسفل المنحنى كنسبة النقاط أسفل الدالة إلى العدد الإجمالي للنقاط. أخيرًا ، نقرب تكامل الوظيفة على أنه المساحة الواقعة أسفل المنحنى مضروبة في القيمة القصوى للدالة ، ونرسم النقاط العشوائية والدالة باستخدام ألوان مختلفة للتمييز بين النقاط الموجودة أسفل الدالة وفوقها. 

لاحظ أن دقة التقريب تعتمد على عدد النقاط المتولدة. ستؤدي زيادة عدد النقاط بشكل عام إلى تقريب أكثر دقة ، ولكنها ستستغرق أيضًا وقتًا أطول للحساب.استخدم الكود في برنامج ار قم بنسخ ولصق الاكواد ثم اعطيه امر Run

 قم بتحميل ملف الاكواد من هنا

بعد تشغيل الاكواد سوف تحصل على الصورة الاتية