New Algorithm for Stochastic Problems with Random Fields of Non-finite Variances
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 MLSRK approximations are derived over local support domains in the probability space. Thus, equating such approximations is still possible, even if the variance of random fields to be studied is infinite. The stochastic problem is next solved with respect to these MLSRK approximations. Testing the succeeding algorithm finds that it doesn’t require many samples and any empirical coefficient to represent accurately random fields following such as Lévy, Cauchy, and multivariate Cauchy distributions. It also provides accurate computation of means and variances of the option price with the stochastic volatility following two empirical Pareto- Lévy and non-stable Lévy distributions. Except for MLSRK approximations of the option price and stochastic volatility, such computation is tested with a deterministic meshless collocation formulation of the Black-Scholes equation.