Fractional Brownian motion stationary and self-similar increments. It has stationary increments: For "H" > ½ the process exhibits long-range dependence, Sample-paths are almost nowhere differentiable. However, almost-all trajectories are locally Hölder continuous of any order strictly less than "H": for each such trajectory, for every "T" > 0 and for every "ε" > 0 there exists a (random) constant "c" such that for 0 < "s","t" < "T". With probability 1, the graph of "B"("t") has both Hausdorff dimension and box dimension of 2−"H". As for regular Brownian motion, one can define stochastic integrals with respect to fractional Brownian motion, usually called "fractional
Fractional Brownian motion Suppose we want to simulate the values of the "fBM" at times formula_12 using the Cholesky decomposition method. In order to compute formula_15, we can use for instance the Cholesky decomposition method. An alternative method uses the eigenvalues of formula_16: Note that the result is real-valued because formula_36. Note that since the eigenvectors are linearly independent, the matrix formula_39 is invertible. It is also known that where "B" is a standard Brownian motion and Where formula_47 is the Euler hypergeometric integral. Say we want to simulate an "fBm" at points formula_48. The integral may be efficiently computed by Gaussian quadrature.