The paper is devoted to the numerical solutions of fractional PDEs based on its probabilistic interpretation, that is, we construct approximate via certain Monte Carlo simulations. main results represent upper bound errors between exact solution and approximation, estimate fluctuation appropriate central limit theorem (CLT) construction confidence intervals. Moreover, provide rates convergence ...