# Stochastic Simulation vs percentile-to-percentile map

I was wondering why someone would go to the trouble to generate random variables in scenarios that are not path dependent. Let me provide a simple (although somewhat contrived) example. Lets say that we want the terminal distribution of a vanilla European call option. We could generate terminal values for the underlying, then plug it into $\mathrm{max}(0;S_T - K)$, and there we have our terminal distribution of the call...

So here's the question. Why not simple do a percentile-to-percentile map from the underlying distribution to the derivatives distribution using equally spaced points?

I do not know if the following is correct, but someone I know told me that we do this because things occasionally fall apart when we use equally spaced points in a multivariate context... I however do not understand why this would happen?

I was hoping someone could shed some light on this issue for me.

What you describe is a very simple quasi monte carlo, where the 'random' points are equally spaced in probability space. Like numerical integration.

Sometimes you can use it, but in general you will need the cumulative distribution to do percentile mapping. This very frequently is not known in closed form, and can be very expensive to compute numerically.

In fact, if you have the cdf in closed form then you will also have the option price in closed form and there is no need to evaluate numerically (Baksi and Madan have shown this amongst others).

If you have the characteristic function in closed form, then it is much more efficient to numerically invert for the option price directly rather than going through probabilities.

Your approach is sensible for a single variate case if the cdf is available, but does (as your friend said) break down for more variates.

One issue with multivariate case is the "curse of dimensionality" - as the number of variates increases your number of samples will get infeasibly large very quickly. To address this, one can use a low discrepancy sequence (eg Sobol) rather than a uniform grid.

Another issue is that even if the terminal marginal distributions are all known, the terminal joint distribution may not be. Of course, one can pick a copula, but this may not be sufficiently accurate.