time-changed Brownian motion

inverse subordinator E 如何在0.9-stable下数值模拟以及相应的时变布朗运动B(E(t))如何模拟