A new algorithm is developed to solve stochastic problems with random fields of non-finite variances. Developing this algorithm motives from an attempt of representing random fields following the Lévy distribution. The first step of current algorithm is deriving moving least square reproducing kernel (MLSRK) approximations of random fields. These...