در این مقاله، حل عددی انتگرال ناسره فازی ~I=ò+¥~a~f (~X)d~x ارایه می گردد که ~f (~X) تابعی با مقدار فازی، بسته و کراندار تعریف شده روی دستگاه اعداد حقیقی فازی بسته است و ~a عددی فازی با تابع عضویت مثلثی با زنگی شکل است. به این منظور، مجموعه های -a برش عدد فازی Ia, ~I و 0£a£1,`Ia که نقاط انتهایی آنها انتگرالهای ناسره قطعی هستند محاسبه می شوند. سپس، یک روش قابل اعتماد برای تخمین Ia و `Iaبا استفاده از حساب تصادفی معرفی می شود. در این حالت، الگوریتمی ارایه می شود که این انتگرالها را با استفاده از قاعده سیمپسون محاسبه می کند. با استفاده از روش CESTAC، اعداد طبیعی بهینه ma و `maرا می یابیم به طوری که این انتگرالهای ناسره توسط انتگرالهای معین تقریب زده می شوند.همچنین، برای یک r مفروض، مقدار تقریبی تابع عضویت m~I® را تقریب زده و دقت نتایج را مشخص می کنیم. سرانجام، دو انتگرال فازی با استفاده از الگوریتم ارایه شده جهت نشان دادن نتایج این تحقیق محاسبه می گردند.