I'm working in double exponential Jump-diffusion model (we know as kou model) with following form , under the physical probability measure $P$: \begin{equation} \frac{dS(t)}{S(t-)}=\mu dt+\sigma dW(t)+d(\sum_{i=1}^{N(t)}(V_i-1)) \end{equation} where $W(t)$ is a standard Brownian motion, $N(t)$ is a Poisson process with rate $\lambda$ , and $\{V_i\}$ is a sequence of independent identically distributed (i.i.d.) non negative random variables such that $Y = log(V)$ has an asymmetric double exponential distribution with the density \begin{equation} f_Y(y)=p.\eta_1 e^{-\eta_{1}y}\upharpoonleft_{y\geq 0}+q.\eta_2 e^{\eta_2 y} \upharpoonleft_{y<0},\eta_{1}>1,\eta_{2}>0 \end{equation} where $p, q \ge 0$, $p+q = 1$, represent the probabilities of upward and downward jumps.
$$ $$ Solving the stochastic differential equation gives the dynamics of the asset price: \begin{equation} S(t)=S(0)\exp\{(\mu- \frac{1}{2}\sigma^2)t+\sigma W(t)\} \prod_{i=1}^{N(t)}V_i \end{equation} and also The stock price process, $(S_t)_{t \geq 0} $, driven by these model, is given by: \begin{equation} S_{t}=S_{0}e^{L_t} \end{equation} where $S_0$ is the stock price at time zero and $L_t$ is defined by: \begin{equation} L_t:=\gamma_{c}t+\sigma W_t+\sum_{i=1}^{N_i}Y_i \end{equation} here,$\gamma_{c}$ is a drift term , $\sigma$ is a volatility, $W_t$ is a Brownian motion, $N_t$ is a Possion process with intensity $ \lambda $, $ Y_i $ is an i.i.d. sequence of random variables.Since $\sigma>0$ in up equation, there exists a risk-neutral probability measure $Q$ such that the discounted process $\{e^{-(r-q)} S_t\}_{t \geq 0}$ becomes a martingale,where $ r $is the interest rate and $ q $ is the dividend rate.Then under this new measure $ Q $, the risk-neutral Levy triplet of $L_t$ can be described as follows: \begin{equation*} (\gamma_{c},\sigma,\nu) \end{equation*} where \begin{align*} \gamma_{c} & = r-q-\frac{1}{2}\sigma^2+ \int_{\mathbb{R}} (e^x-1) \nu(dx) \\ & = r-q-\frac{1}{2}\sigma^2+ \lambda \eta \end{align*} Here we focus on the case where the Levy measure is associated to the pure-jump component and hence the Levy measure$ \nu(dx) $ can be written as $ \lambda f(x) dx $, where the weight function $ f(x) $ can take the following form: \begin{equation} f(x):=p.\eta_1 e^{-\eta_{1}x}\upharpoonleft_{x\geq 0}+(1-p).\eta_2 e^{\eta_2 x} \upharpoonleft_{x<0},\eta_1>1,\eta_2>0 \end{equation}
Also in $ \eta = \int_{\mathbb{R}}(e^x-1)f(x) dx $ represents the expected relative price change due to a jump. Since we have defined the Levy density function $ f(x) $ for double exponential Jump-diffusion model, $ \eta $ can be computed as: \begin{equation} \eta= \frac{p \alpha_1}{\alpha_{1}-1}+\frac{(1-p)\alpha_2}{\alpha_2+1}-1 \end{equation}
This is found by integrating $ e^x $ over the real line by setting $ \alpha_1 >1 $ and $ \alpha_{2}>0 $.
We let $\tau=T-t$, the time-to-maturity, where $T$ is the maturity of the financial option under consideration and we introduce $x = log S_t$, the underlying asset's log-price. If $u(x; \tau )$ denotes the values of some (American and European) contingent claim on $S_t$ when $log St = x$ and $\tau = T - t$, then it is well-known, see for example, (Cont and Tankov, 2004) that $u$ satisfies the following $PIDE$ in the non-exercise region: \begin{align*} \partial_\tau\, u(x,\tau) & = \frac{1}{2}\sigma^2 \partial_{x}^2 u +(r-q-\frac{1}{2}\sigma^2-\lambda \eta)\partial_x u-(r+\lambda)u \\ &+ \lambda \int _{\mathbb{R}} u(x+y,\tau) f(y) dy \end{align*} with initial value \begin{equation} u(x,0)=g(x):=G(e^x)= \begin{cases} max\{e^x-k,0\}, & \text{call option} \\ max\{k-e^x,0\}, & \text{put option} \end{cases} \end{equation}
my question is how we can derive the above $PIDE$ I've searched a lot of article but most of them only mention the $PIDE$ and we said you can find in Cont and Tankov Book and also I've searched in this book but I could not find the Exactly above $PIDE$.
thanks for help.
