Let's consider a Vanilla European Call with a Strike K and Maturity T. Its payoff function is as below:
Assume the underlying asset follows GBM process:
The Fokker Planck Equation with respect to this process is:
Integrate both side of with respect to :
where is the density function of at time
Take the derivative of with respect to
replace using :
We can use Integration by part on rhs separately. For the first part of rhs:
For the second part of , we also do integration by parts. However, note that there is a second order differential term, so we have to do integration by parts twice. Here I just skip the calculation:
With and , can be written as :
Note that can be written as:
Also, we take the derivative of with respect to K:
Stochastic Volatility Modeling by Bergomi, Lorenzo
2. The Volatility Surface A Practitioners Guide (Wiley Finance) by Jim Gatheral, Nassim Nicholas Taleb
3. Local Volatility Model: Dupire PDE and Valuation/Pricing PDE Derivations and Comparisons - YouTube