Considere la función de recompensa $$ V_T = max(S_T^1 - K, 0) 1_{\{L<S_T^2<U\}} = (S_T^1 - K)1_{\{S_T^1 > K\}}1_{\{L<S_T^2<U\}}$$
donde $S_T^1$ y $S_T^2$ son dos valores distribuidos por GBM con correlación, $\rho$ . ¿Cómo encontrarías $V_t$ sin utilizar la simulación MC (pista: utilizar una cópula gaussiana)?
Mi intento: $$V_t = e^{-r(T-t)}E[(S_T^1 - K)1_{\{S_T^1 > K\}}1_{\{L<S_T^2<U\}}]$$ $$V_t = e^{-r(T-t)}\int_{L}^{U}\int_{K}^{\infty}(x-K)f_{x,y}(x,y)dxdy$$
donde $f_{x,y}$ es el pdf conjunto procedente de la cópula guasiana. ¿Cómo puedo ir más allá ahora? Si no puedo, ¿cómo implementaría la integral doble en un ordenador?