La respuesta de Gordon es agradable (+1). Quiero añadir la otra solución.
Observación 1
Dejemos que $X=(X_1,X_2,...,X_n)$ donde el componente $X_i$ tiene una diferencial estocástica de la forma $$dX_i(t)=\mu_i(t)dt+\sum_{j=1}^{d}\sigma_{ij}(t)dW_j(t)$$ donde $dW_k(t)dW_j(t)=\rho_{kj}dt$ para todos $k,j\in\{1,2,...,d\}$ . Sea $f:\mathbb{R}^+\times\mathbb{R}^n\to\mathbb{R}\in\mathbb{C}^{1,2}.$ Por aplicación del lema de Ito, tenemos $$df(t,X_1,...,X_n)=\frac{\partial f}{\partial t}dt+\sum_{i=1}^{n}\frac{\partial f}{\partial x_i}dX_i+\frac 12\sum_{i=1}^{n}\sum_{l=1}^{n}\frac{\partial ^2f}{\partial x_i\partial x_l}dX_idX_l\tag 1$$
Observación 2
Supongamos que el tipo corto $r_t$ sigue el proceso Ito descrito por la siguiente ecuación diferencial estocástica $$dr_t=\mu(t,r_t)dt+\sigma(t,r_t)dW_t$$ y $P(t,T)$ denota el precio del bono de cupón cero con vencimiento $T$ . Podemos mostrar $$\frac{\partial P}{\partial t}+\mu(t,r_t)\frac{\partial P}{\partial r}+\frac{1}{2}\sigma^2(t,r_t)\frac{\partial^2 P}{\partial r^2}-r_tP=0\tag 2$$
Observación 3
Dejemos que $$\qquad dS_t=rS_t+\sigma S_tdW_1(t)\\ \quad\qquad\quad dr_t=\kappa(\theta-r_t)dt+\Sigma dW_2(t)\\ dW_1(t)dW_2(t)=0\tag 3$$
Ahora formamos una cartera compuesta por una opción $V=V(t,S,r,T,K)$ (Posición corta), $\Delta_1$ unidades de la acción (posición larga) y $\Delta_2$ unidades del $T-$ precio del bono de cupón cero (posición larga). La cartera tiene valor $$\Pi=\Delta_1S_t+\Delta_2 P(t,T)-V(t,S,r,T,K)$$ por lo tanto $$d\Pi=\Delta_1dS_t+\Delta_2dP(t,T)-dV(t,S,r,T,K)\tag 4$$ Por aplicación del lema de Ito, tenemos $$dV_t=\frac{\partial V}{\partial t}dt+\frac{\partial V}{\partial s}dS+\frac{\partial V}{\partial r}dr+\frac 12\left(\sigma^2\frac{\partial^2 V}{\partial s^2}+\Sigma^2\frac{\partial^2 V}{\partial r^2}\right)dt\tag 5$$ y $$dP=\frac{\partial P}{\partial t}dt+\frac{\partial P}{\partial r}dr+\frac{1}{2}\Sigma^2\frac{\partial^2 P}{\partial r^2}dt\tag 6$$ $(4)\,,\,(5)\,,(6)$ y $(3)$ $$d\Pi =-\left( \frac{\partial V}{\partial t}+\frac{1}{2}{{\sigma }^{2}}{{S}^{2}}\frac{{{\partial }^{2}}V}{\partial {{S}^{2}}} \right)dt+\left( \Delta _1-\frac{\partial V}{\partial S} \right)dS\\+\left( \Delta _2\frac{\partial P}{\partial r}\,-\frac{\partial V}{\partial r} \right)dr+\Delta _2\,\left( \frac{\partial P}{\partial t}+\frac{1}{2}{{\Sigma }^{2}}\frac{{{\partial }^{2}}P}{\partial {{r}^{2}}} \right)dt\,$$ A continuación, encontramos los valores de $\Delta_1$ y $\Delta_2$ que hace que la cartera no tenga riesgo. En efecto, fijamos $$ \Delta_1=\frac{\partial V}{\partial S}\\ \Delta_2=\frac{\frac{\partial V}{\partial r}}{\frac{\partial P}{\partial r}}$$ así $$d\Pi =-\left( \frac{\partial V}{\partial t}+\frac{1}{2}{{\sigma }^{2}}{{S}_{t}}^{2}\frac{{{\partial }^{2}}V}{\partial {{S}^{2}}}+\frac{1}{2}{{\Sigma }^{2}}\frac{{{\partial }^{2}}V}{\partial {{r}^{2}}} \right)dt+\Delta {{}_{2}}\,\left( \frac{\partial P}{\partial t}+\frac{1}{2}{{\Sigma }^{2}}\frac{{{\partial }^{2}}P}{\partial {{r}^{2}}} \right)dt\tag 7$$ $(2)$ y $(7)$ $$d\Pi =-\left( \frac{\partial V}{\partial t}+\frac{1}{2}{{\sigma }^{2}}{{S}^{2}}\frac{{{\partial }^{2}}V}{\partial {{S}^{2}}}+\frac{1}{2}{{\Sigma }^{2}}\frac{{{\partial }^{2}}V}{\partial {{r}^{2}}} \right)dt+\Delta {{}_{2}}\,\left( {{r}_{t}}\,P-\kappa \left( \theta -r \right)\,\frac{\partial P}{\partial r} \right)dt\tag 8$$ La condición de que la cartera gane la tasa libre de riesgo, $r$ implica que la variación del valor de la cartera es $$d\Pi=r\Pi dt$$ En otras palabras $$d\Pi =(\Delta _1r_t\,S_t+\Delta_2r_tP-r_tV)dt\tag 9$$ $(8)$ y $(9)$ $$\frac{\partial V}{\partial t}+{{r}_{t}}\,{{S}_{t}}\frac{\partial V}{\partial S}+\kappa \left( \theta -r \right)\frac{\partial V}{\partial t}+\frac{1}{2}{{\sigma }^{2}}{{S}^{2}}\frac{{{\partial }^{2}}V}{\partial {{S}^{2}}}+\frac{1}{2}{{\Sigma }^{2}}\frac{{{\partial }^{2}}V}{\partial {{r}^{2}}}-{{r}_{t}}V=0$$