This paper presents an extension of a general computational scheme for asymptotic expansions proposed in earlier works by the authors and coworkers. In the earlier works, a new method was developed for the computation of an arbitrary-order expansion with a normal benchmark distribution in a multidimensional diffusion setting. In particular, a new algorithm was proposed for calculating coefficients in an expansion by solving a system of ordinary differential equations. In the present note, by a change of variable technique, and by various ways of setting the perturbation parameters in the expansion, we provide the flexibility of setting the benchmark distribution around which the expansion is made and an automatic way for computation up to any order in the expansion. For instance, we introduce new expansions, called the lognormal expansion and the CEV expansion. We also show some concrete examples with numerical experiments, which imply that a high-order CEV expansion will produce more a precise and stable approximation for option pricing under the SABR model than other approximation methods such as the log-normal expansion and the well-known normal expansion. In the application of the asymptotic expansion based on Watanabe theory, we need to calculate certain conditional expectations which appear in the expansions and play a key role in computation. In the first place [52,53] have developed the formulas necessary for the second-order expansion. Subsequently, [33,34,[38][39][40] have derived new formulas up to the third order. (Also, multi-dimensional formulas were provided in [33,34].) In many applications, these formulas give sufficiently accurate approximation, but in some cases such as in the cases with long maturities or/and with highly volatile underlying variables, the approximation up to the third order may not provide satisfactory accuracies. Thus, the formulas for the higher order computation are desirable.Recently, [41,42,48] have proposed two alternative computational schemes for any order expansions in an automatic manner. In fact, one of their new methods does not rely on direct evaluation of the conditional expectations, but on solving a certain system of ordinal differential equations with grading structure. Independently, [19] has developed a new computational method for the conditional expectations necessary for high order expansions. As a consequence, their approximations generally showed sufficient accuracy with computation of high order expansions, which was confirmed by numerical experiments.Furthermore, in terms of approximation it is important to set the limiting or benchmark distribution around which an expansion is made. Typically a normal distribution is chosen, which enables us to compute higher order correction terms due to the nice and well-known Gaussian properties. For (shifted) log-normal local volatility cases in (jump-diffusion) stochastic volatility models, the same technique is applied. (For instance, see [15,29,31,44,48].) This paper introduces a change of variable technique i...