A PDE method for estimation of implied volatility

In this paper it is proved that the Black–Scholes implied volatility satisfies a second order non-linear partial differential equation. The obtained PDE is then used to construct an algorithm for fast and accurate polynomial approximation for Black–Scholes implied volatility that improves on the existing numerical schemes from literature, both in speed and parallelizability. We also show that the method is applicable to other problems, such as approximation of implied Bachelier volatility.