Suponemos que el proceso del precio de las acciones $\{S_t,\,t>0\}$ satisface, bajo la medida de probabilidad del mundo real $P$ una SDE de la forma \begin{align*} dS_t=S_t\big((\mu-q)dt+\sigma dW_t\big), \end{align*} donde $\{W_t, \, t >0\}$ es un movimiento browniano estándar. En este caso, tenemos que considerar el activo de retorno total $e^{qt}S_t$ es decir, el activo con los pagos de dividendos invertidos en la misma acción subyacente. Consideramos una cartera de autofinanciación localmente libre de riesgo de la forma \begin{align*} \pi_t =\Delta_t^1 \big(e^{qt}S_t\big) + \Delta_t^2 V_t, \end{align*} donde $V_t$ es el precio de la opción. Entonces, \begin{align*} d\pi_t &= \Delta_t^1 d\big(e^{qt}S_t\big) + \Delta_t^2 dV_t\\ &= \Delta_t^1 e^{qt}\big(q S_t dt + dS_t \big) + \Delta_t^2\left(\frac{\partial V}{\partial t}dt + \frac{\partial V}{\partial S}dS_t + \frac{1}{2}\frac{\partial^2 V}{\partial S^2} \sigma^2S_t^2 dt\right)\\ &=\left[\mu\Delta_t^1 e^{qt} S_t + \Delta_t^2\left(\frac{\partial V}{\partial t} + (\mu-q) S_t \frac{\partial V}{\partial S} + \frac{1}{2}\frac{\partial^2 V}{\partial S^2} \sigma^2S_t^2 \right)\right]dt \\ &\qquad\qquad\qquad\qquad\qquad\quad + \left(\sigma\Delta_t^1 e^{qt}S_t + \sigma \Delta_t^2 S_t \frac{\partial V}{\partial S}\right)dW_t. \end{align*} Desde $\pi_t$ está localmente libre de riesgo, suponemos que $\pi_t$ gana el tipo de interés sin riesgo $r$ Es decir, \begin{align*} d\pi_t = r \pi_t dt, \end{align*} Entonces, \begin{align*} &\left[\mu \Delta_t^1 e^{qt} S_t + \Delta_t^2\left(\frac{\partial V}{\partial t} + (\mu-q) S_t \frac{\partial V}{\partial S} + \frac{1}{2}\frac{\partial^2 V}{\partial S^2} \sigma^2S_t^2 \right)\right]dt \\ &\qquad\qquad\qquad\qquad\qquad + \left(\sigma\Delta_t^1 e^{qt} S_t + \sigma \Delta_t^2 S_t \frac{\partial V}{\partial S}\right)dW_t= r \pi_t dt. \end{align*} En consecuencia, \begin{align*} \sigma\Delta_t^1 e^{qt}S_t + \sigma \Delta_t^2 S_t \frac{\partial V}{\partial S}=0, \tag{1} \end{align*} y \begin{align*} \mu e^{qt} \Delta_t^1 S_t + \Delta_t^2\left(\frac{\partial V}{\partial t} + (\mu-q) S_t \frac{\partial V}{\partial S} + \frac{1}{2}\frac{\partial^2 V}{\partial S^2} \sigma^2S_t^2 \right) = r(\Delta_t^1 e^{qt}S_t + \Delta_t^2 V_t). \end{align*} Desde $(1)$ , \begin{align*} \Delta_t^1 = -e^{-qt} \Delta_t^2 \frac{\partial V}{\partial S}. \end{align*} Entonces, \begin{align*} -\mu \Delta_t^2 S_t \frac{\partial V}{\partial S}+ \Delta_t^2\left(\frac{\partial V}{\partial t} + (\mu-q) S_t \frac{\partial V}{\partial S} + \frac{1}{2}\frac{\partial^2 V}{\partial S^2} \sigma^2S_t^2 \right) = r\Big(-\Delta_t^2 S_t\frac{\partial V}{\partial S} + \Delta_t^2 V_t\Big), \end{align*} o \begin{align*} \Delta_t^2\left(\frac{\partial V}{\partial t} -q S_t \frac{\partial V}{\partial S} + \frac{1}{2}\frac{\partial^2 V}{\partial S^2} \sigma^2S_t^2\right) &= r\Delta_t^2\Big(-\frac{\partial V}{\partial S} S_t + V_t\Big). \tag{2} \end{align*} Anulación del plazo $\Delta_t^2$ de ambos lados de $(2)$ obtenemos la ecuación de Black-Scholes de la forma \begin{align*} \frac{\partial V}{\partial t} + (r-q) S_t \frac{\partial V}{\partial S} + \frac{1}{2}\frac{\partial^2 V}{\partial S^2} \sigma^2S_t^2 -rV = 0. \end{align*}