In this paper, the numerical solution of the improper fuzzy integral ~I=ò+¥~a ~¦(~x)~dx is proposed where, ~¦(~x) is a bounded and closed fuzzy valued function defined on the closed fuzzy real number system and a is ~a fuzzy number with triangular or bell-shape membership function. For this purpose, the a -level sets of the fuzzy number ~I, -Ia and -Ia, 0£a£1, are evaluated whose end points are the crisp improper integrals. Then, a reliable method is introduced to estimate -Ia and -Ia by using the stochastic arithmetic. In this case, an algorithm is presented which evaluates these integrals by using the Simpson rule. By using the CESTAC method, we find the optimal natural numbers -ma and -ma such that these improper integrals are estimated by definite integrals. Also, for a given r, we evaluate the approximate value of the membership function m~I (r) and determine the accuracy of the results. At last, two fuzzy integrals are computed using the given algorithm to show the results of the research.