Given constants $r,\sigma, K>0$, considering the Black-Scholes PDE of a Europoean call: \begin{cases} \frac{\partial c}{\partial t}+rs\frac{\partial c}{\partial s}+\frac{1}{2}\sigma^2s^2 \frac{\partial^2 c}{\partial s^2}-rc = 0\\ c(T,s) =(s-K)^+\\ c(t,0) = 0 \end{cases} I learned that we can do the substitution like this: \begin{cases} \tau=T-t\\ u=ce^{r\tau}\\ x=\ln\frac{s}{K}+(r-\frac{1}{2}\sigma^2)(T-t) \end{cases} Set $u=u(\tau,x)$, then $c(t,s)=e^{-r\tau}u(\tau,x)$. Then we have: \begin{gather*} \frac{\partial c}{\partial t} = re^{-r\tau}u+e^{-r\tau}(\frac{\partial u}{\partial \tau} \frac{d \tau}{dt} + \frac{\partial u}{\partial x} \frac{\partial x}{\partial t}) = e^{-r\tau}[ru-\frac{\partial u}{\partial \tau}- (r-\frac{1}{2}\sigma^2)\frac{\partial u}{\partial x}]\\ \frac{\partial c}{\partial s} = e^{-r\tau}\frac{\partial u}{\partial x} \frac{\partial x}{\partial s} = e^{-r\tau}\frac{1}{s}\frac{\partial u}{\partial x}\\ \frac{\partial^2 c}{\partial s^2} = e^{-r\tau}[-\frac{1}{s^2}\frac{\partial u}{\partial x}+\frac{1}{s^2}\frac{\partial^2 u}{\partial x^2}]=e^{-r\tau}\frac{1}{s^2}(\frac{\partial^2 u}{\partial x^2}-\frac{\partial u}{\partial x}) \end{gather*} Thus: \begin{align*} \frac{\partial c}{\partial t}+rs\frac{\partial c}{\partial s}+\frac{1}{2}\sigma^2s^2 \frac{\partial^2 c}{\partial s^2}-rc &= e^{-r\tau}[-\frac{\partial u}{\partial \tau}- (r-\frac{1}{2}\sigma^2)\frac{\partial u}{\partial x} + r\frac{\partial u}{\partial x} + \frac{1}{2}\sigma^2(\frac{\partial^2 u}{\partial x^2}-\frac{\partial u}{\partial x})]\\ &= e^{-r\tau}[-\frac{\partial u}{\partial \tau} + \frac{1}{2}\sigma^2\frac{\partial^2 u}{\partial x^2}]\\ &=0 \end{align*} Namely we get the heat equation with boundary condition: \begin{cases} \frac{\partial u}{\partial \tau} = \frac{\sigma^2}{2}\frac{\partial^2 u}{\partial x^2}\\ u(0,x)=K(e^x - 1)^+\\ u(\tau, -\infty) = 0 \end{cases} I am stuck here. Since the boundary condition is not integrable, I cannot use Fourier transformation with regard to $x$ here. I know there is a probabilistic way of doing it, by computing the conditional expectation. I want to solve the PDE directly. After getting the heat equation here, what should I do next to solve it?
Besides, is there any way to start directly from the PDE of $c$, without change of variable? And, how do we know that change of variable work here? It seems to me very coincidentally the PDE becomes a heat equation in the end.
Thank you so much!
You can apply the Fourier transform with complex wavenumber, $k = k_r + ik_i$.
Since $u(0,x) = \mathcal{O}(e^x)$ as $x \to +\infty$ and $u(0,x) = 0$ for $x \leqslant 0$ the Fourier integral (applied to the initial condition) converges when $\Im{(k)} > 1$,
$$\hat{u}(0,k) = \int_{-\infty}^\infty u(0,x) e ^{ikx} \, dx = \int_{-\infty}^\infty \underbrace{u(0,x) e^{-k_i x}}_{\in L^1(\mathbb{R})}e ^{ik_rx} \, dx \\ = \int_{0}^\infty K(e^x-1)e ^{ikx}\,dx = -\frac{K}{1 + ik} + \frac{K}{ik} \\ = \frac{K}{ik-k^2}$$
Apply the transform to the PDE in the usual way and obtain an ODE for the transform $\hat{u}(\tau,k)$ of the form
$$\frac{\partial \hat{u}}{\partial \tau} = - \frac{\sigma^2k^2}{2} \hat{u},$$
with the solution
$$\hat{u}(\tau,k) = \hat{u}(0,k) e^{-\sigma^2k^2\tau/2} = \frac{Ke^{-\sigma^2k^2\tau/2}}{ik - k^2}$$
The inverse transform takes the form of a contour integral in the complex plane
$$u(\tau,x) = \frac{1}{2\pi}\int_{i\beta-\infty}^{i\beta+\infty} \hat{u}(\tau,x) e ^{-ikx} \, dk, $$
for any $\beta > 1$, leading ultimately to the Black-Scholes solution.