In this paper, by using the CESTAC method and the CADNA library a procedure is proposed to solve a fuzzy initial value problem based on RBF-meshless methods under generalized H-differentiability. So a reliable approach is presented to determine optimal shape parameter and number of points for RBF-meshless methods. The results reveal that the proposed method is very effective and simple. Also, the numerical accuracy of the method is shown in the tables and figures, and algorithms are given based on the stochastic arithmetic. The examples illustrate the efficiency and importance of using the stochastic arithmetic in place of the floating-point arithmetic.